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ABSTRACT 


This  report  develops  a dynamia  map  projection  especially  suited 
for  processing  and  display  of  satellite  electro-optical  remote  sensing 
of  the  earth's  surface.  The  new  map  projection  (the  Space  Oblique 
Mercator)  projects  the  satellite  ground-track  from  the  ellipsoid  into 
the  map  plane,  free  of  length  distortion  and  free  of  normal  view 
curvature  distortion.  The  length  and  curvature  distortions  in  the 
finite  sensed  region  are  negligible  for  most  applications.  The  report 
details  the  formulation,  provides  numerical  examples  for  the  LANDSAT 
multi-spectral  scanner,  and  includes  FORTRAN  IV  software  as  an  appendix. 
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1,0  SUMMARY 


We  report  herein  upon  a research  and  development  effort  which  has 
led  to  the  following  results: 

• A rigorous  mathematical  development  of  a Space  Oblique  Mercator 
map  projection.  This  projection  is  especially  designed  for 
processing  and  displaying  data  from  LANDSAT-type  satellite- 
borne  multi-spectral  scanners  and  is  characterized  by  the 
following  desirable  features: 

• Zero  length  and  angle  distortion  along  the  satellite 
groundtrack,  and  only  a few  parts  in  10,000  length 
distortion  as  far  as  200  km  off  typical  groundtracks . 

• Rigorously  valid  for  an  arbitrary  continuous  satellite 
orbit;  the  formulation  can  be  routinely  interfaced  to 
state-of-the-art  orbit  integration  programs,  or  can  use 
simplified  Kepler ian  or  circular  orbit  approximations, 
depending  upon  the  needs  of  particular  applications. 

• Rigorously  valid  for  an  arbitrary  ellipsoid  reference 
figure  for  the  earth. 

• Computationally  efficient;  the  most  expensive  calculations 
are  orbit-dependent  integrals  which  need  only  be  determined 
once  (for  each  specific  nominal  orbit). 

• Prototype  software  (FORTRAN  IV)  has  been  developed,  checked  out, 
and  is  demonstrated  herein.  The  software  has  been  developed 
and  all  calculations  performed  on  the  University  of  Virginia's 
CYBER  172  computer  system. 
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2.0  PREFACE 


This  report  constitutes  the  final  report  of  Phase  III  of  contract 
no.  DAAG-53-76-C-0067  performed  by  the  University  of  Virginia  for  the 
U.S.  Army  Engineer  Topographic  Laboratories,  Fort  Belvoir,  Virginia, 
under  the  sponsorship  of  the  U.S.  Geological  Survey,  (USGS) , Reston, 
Virginia . 

The  authors  acknowledge  the  competent  guidance  and  technical 
liasion  of  Mr.  L.A.  Gambino  (technical  mointor,  USAETL)  and  Dr.  A.P. 
Colvocoresses  (EROS  Program  Cartography  Coordinator,  USGS). 
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3.0  INTRODUCTION 


In  Ref,  1,  Colvocoresses  conceived  a dynamic  map  projection  concept 
especially  suited  to  satellite  mapping.  Unlike  classical  stat^a  map 
projections,  the  line  which  is  projected  distortion-free  is  not  restricted 
to  be  an  equator,  a meridian,  a parallel,  or  an  oblique  great  circle; 
rather,  the  distortion-free  line  is  the  satellite  sub-point  path  (ground- 
track)  on  a reference  ellipsoid.  Colvocoresses  developed  a geometrical 
analog  involving  an  oscillating  cylinder  to  which  projections  are  made 
from  the  reference  ellipsoid,  the  cylinder  oscillation  is  such  that  the 
cylinder  instantaneously  osculates  with  the  normal  sub-point  on  the 
ellipsoid.  A small  region  near  the  sub-point,  when  projected  from  the 
ellipsoid  onto  the  oscillating  cylinder  and  developed  onto  a plane,  is 
projected  with  negligible  length  and  angle  distortions. 

Unfortunately,  Colvocoresses'  elegant  geometric  analog  was  not 
supported  by  a mathematical  formulation  of  the  map  projection  equations; 
rather,  he  Issued  a challenge"! or  the  cartographic  community  to  undertake 
a considerable  and  dedicated  effort  to  develop  the  mathematical  model 
and  associated  computer  programs"  to  implement  this  map  projection 
concept.  The  present  report  documents  our  response  to  the  above 
challenge.  While  motivated  (and  occasionally  perplexed!)  by  the 
oscillating  cylinder  analog,  we  have  not  used  this  concept  in  our 
formulations.  We  elected  instead  to  set  doiim  the  mathematical  constraints 
underlying  Colvocoresses'  objectives.  These  lead  immediately  to 
differential  equations  which  are  the  key  to  a rigorous  solution  to  the 
problem. 
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We  have  successfully  developed  a most  general  version  of  this  space 
oblique  mercator  (SOM)  map  projection  formulation  (really,  an  infinite 
family  of  map  projections,  depending  upon  definition  of  the  nominal 
orbit).  The  derivation,  computational  summary,  and  software  for  the 
SOM  projection  are  documented  herein. 

The  present  report  is  arranged  in  a modular  fashion.  The  computational 
summary  concentrates  on  the  structure  of  the  solution  and  relegates 
details  to  appendixes,  Exclusive  of  Section  6,  the  report  is  directed 
primarily  toward  readers  seeking  to  understand  the  essence  of  and  learn 
to  use  the  map  projection.  Section  6 is  intended  to  document  the 
geometrical,  mathematical,  and  intuitive  details  and  to  develop  the 
equations  in  a fashion  which  parallels  their  invention. 
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4,0  SOM  C'"»PUTATIONAL  SUMMARY 

4.1  Forward  Transformation 

The  map  coordinates  (x.y)  are  related  to  the  corresponding 

ellipsoidal  coordinates  (()),X)  by  the  formulae 
t* 

X = f Vcosf  dt  + R ln[tan(y  + ■^)]cos  Y /cos9  (la) 

' C 4Z  S S 

O 

t* 

y = / Vsin  f dt  + R^ln[tan (-^  + -I-) ] sin  Yg/cosB^  + y (0)  (lb) 

where 


t = time  since  some  selected  initial  point  in  the  orbit 
t*=  the  "scan  instant"  for  which  the  scan  vector  passes  through  ((?,a). 
V = satellite  sub-point's  instantaneous  speed  relative  to  the  earth 
p = the  instantaneous  radius  of  curvature  of  the  sub-point  path,  as 

projected  into  the  plane  tangent  to  the  ellipsoid  at  the  sub-point 


f S I — dt : T = dummy  (time)  integration  variable. 

R^=  the  local  radius  of  curvature  of  the  ellipsoid  in  a plane 


(2) 


whose  normal  is  the  earth-fixed  velocity  vector  of  the  sub-point 

(a,Y  ,6  ) angles  defined  in  §6.3. 
s s 

Methods  for  evaluating  the  right-hand  side  of  equations  (1)  are  now 
discussed . 

t t t y 

4.1.1  Calculation  of  the  Integrals  / Vcos  f dt , / Vsin  f dr , / — dr 

0 O 0 ^ 

As  will  be  seen  in  §6.1,  the  integral  terms  of  equations  (la)  and 

(lb)  are  correctly  interpreted  as  (x  ,y  ),  the  map  coordinates  of  the 

S 8 

groundtrack  or  sub-point  path.  We  can  write  Instead  of  the  three  integrals 
the  three  differential  equations 
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(3a) 


dx 

X = ^ = Vcos  f 

g dt 


dy 

y = — ^ = Vsin  f 
■’g  dt 

(3b) 

■ df  _ V 

^ ' dt  p 

(3c) 

Given  the  appropriate  formulas  (Appendixes  A and  B)  for  calculation  of 


V and  p,  the  given  initial  conditions  {x^(o),  y^(o) , f(o)),  then 

equations  (3)  can  be  Integrated  numerically  to  evaluate  the  integrals 

in  equations  (1)  and  (2).  However,  it  is  expensive  to  carry  out  these 

integrations  many  times,  and  since  these  integrals  have  been  found  to  be 

very  smooth  functions  of  time,  they  can  be  conveniently  and  accurately 

replaced  by  their  harmonic  series  as 

00  t 

X (t)  = X t + y {S  sin(-^^^)  + C cos  (■^^^) } = / V cosf  dT(')a) 
g g ^ „ x_n  P X n P { 

O 


n=0  g 


y (t)  = y C + y {S  + C cos  / V sinf  dxC^b) 

2 2 P yn  P ^ 


n=0 


n=0 


f(t)  = f t ^ I sin  C-^)  + cos  (^)}=  ! - dr 


where 


(4c) 


P = orbital  period 

_ dx  P dx 

X = the  average  value  of  — rf-  = rr  j 3— ^(  r)  di , x e-y 
g dtP^'^dT  gg 


(5) 


The  symbol  (interpreted  as  "replaced  by")  allows  the  three  equations 

, P dx 


X = 


f = 


0 dt 


P dy 


1 f .df,  . 
p ! 'dt'-*' 


(footnote  continued  on  next  page) 
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(6) 


gn  o ^ S 


T , X -*-y  ,f 
g g 


*^x  = 7 / - X t:]cos(^)  dx,  ^ ^y  , f (7) 

gn  o ® ^ & & 

For  certain  approximate  orbits  and  choice  of  the  initial  point  (position  in 
the  orbit  for  which  t=0) , many  of  the  above  coefficients  are  either  zero  or 
negligible.  For  example,  for  circular  LANDSAT  orbits,  (99°  inclination,  103 
min.  period)  equations  (4)  have  been  found  to  reduce  to  simply 


X (t)  = X t + A sin  (-^)  + A sin(^^)+  . 
g g 1 P 2 P 


, . . f'9TTt  ^ 

. + AgSin(— ) 


yg(t)  = + B^cos(^)  + B2Sin('^^)+  . . . . + B^sin  (•^^) 

f(t)  = C^sin(^)  + C2sin(-^~)+  . . . + C^sin  (•^^) 

where  the  specialized  version  of  equations  (5) , (6) , and  (7)  are 
_ 1 P dx 

X = — I S.  dx 

g p J a-r 


^ = P- 


/^[x  (x)  - X x]sin(^)di 


^ p J cos(^^)dt  (9c) 

^n  = I ^0  sin(5^)dx  (9d) 

and  with  t=0  the  instant  for  which  the  satellite  is  at  its  northernmost 
lattitude.  In  eqs.  8,  all  but  two  coefficients  arc  near-negligible 


in  each  of  the  three  series  (see  §5.1). 

The  solution  procedure  (restricting  the  discussion  to  the  LANDSAT  case) 

for  the  coefficients  (9)  is  as  follows: 

• Using  the  Runge-Kutta  algorithm  of  Appendix  B and  the  calculations 

leading  to  instantaneous  values  for  V and  f (established  in  Appendix 

B) , calculate  x by  numerical  Integration  of  the  differential 


footnote  cont'd.:  to  be  written  as  only  the  first,  eq.  (5),  since  replacing 
By  y yields  the  second  and  replacing  x by  f yields  the  third.  This 

shorthand  notations  is  used  in  all  subsequent  equations  for  compactness. 
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equation 


^ (Xg)  = (i)  Vcos  f (10) 

using  zero  as  the  initial  condition. 

■ The  integrals  (9)  are  evaluated  by  simultaneous  Runge-Kutta 

solution  (Appendix  C)  of  the  following  system  of  (3+3n)  differential 
equations 


dt  p 
dx 

— = Vcos  f 


(lla) 

(llb) 


dy 

— j-j  = Vsin  f (11c) 


dA 
n 

dt 

dB 
n 

dt 

dC 
n 

dt 


2 

= p •^]sin 

= p Yg  (t)  cos  (— -) 

= - f (t)  sin  0-^) 


(lid) 

(lie) 


(Ilf) 


using  zeros  for  initial  values. 

Analogous  integrations  establish  the  coefficients  for  the  general 
case  of  equations  (5),  (6),  and  (7). 

4.1.2  Scan  Time  Determination 

For  given  ellipsoidal  coordinates  (‘t’,X),  the  corresponding  scan  time 
t*  is  defined  as  the  Instant  that  the  plane  established  by  the  unit  vector 
_n  (normal  to  the  ellipsoid)  and  the  scan  vector  w (tangent  plane 
projection  of  the  orbit  normal)  contain  the  vector  [from  the  sub-point 
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(<(>  ^ ) to  the  point  of  interest  on  the  ellipsoid];  this  geometrical 

g.  g 

condition  requires  that  these  three  vectors  satisfy  the  constraint 
F(t*)  = [(w  X n)  • = 0 

A/s  AA  AA  AA  AA  AA 

= [AR  (w  n - w n ) + AR  (w  n - w n ) + AR  (w  n - w n )] 

x yz  zv  vzx  xz  zxv  yxt- 


z y 


y z X 


y X'  ■t=t* 
(12) 


This  constraint  neglects  the  finite  scanner  sweep  time;  upon  finding 
a t*  satisfying  (12) , an  additive  correction  (13b)  is  introduced  to 
account  for  the  scanner  sweep  time. 

In  Appendix  B,  the  explicit  earth-fixed  components  of  w and  n^  are 
given,  and  the  components  of  AR  are  developed  in  §6.3;  these  expressions  (6, 
and  their  derivatives  are  required  to  calculate  F(t)  and  . The 

time  t*  for  which  (12)  vanishes  is  determined  via  a Newton's  successive 


approximation  algorithm  as 

.(k+1) 


(k)  F[t^] 


dt' 


(13a) 


t=t 


(k) 


(o) 


with  the  approximate  estimates  t determined  by  calculating  the  angle 
A9  between  the  initial  sub-point  position  vector  and  the  position  vector 
to  (♦>,X),  then  dividing  by  2tt/P:  t^°^  = . with  appropriate 


quadrant  checks. 


The  condition  (12)  and  the  resulting  converged  t*  from  (13a)  must 

be  corrected  to  account  for  the  finite  scan  time  At  for  thi  scanner 

to  sweep  from  the  center  of  scan  [t(ij)  ,X  )]out  to  the  sensed  point  at 

g S 

'i^j.X).  This  correction,  added  to  (13a),  has  the  form 


At  = T (f~) 
max 


(13b) 


13) 
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where 


T = % the  scan  period 

= 18.355  milliseconds  (for  the  LANDSAT  scanner) 

c = the  maximum  scanner  beam  deflection  angle  away  from  local  vertical 

= 5.78  degrees  (for  the  LANDSAT  scanner) 

£ = instantaneous  scanner  beam  deflection  angle  away  from  the  local 

vertical 

-1 

= cos"  ( — g^)  • sign  (AR-£)  (13c) 

where  from  figures  6.4C,  6.4D,  the  displacement  of  the  sensed  point 
from  the  satellite  is, 

S = R - £ 

Equation  (13b)  assumes  a constant  linear  scan  rate  and  that  center-of-scan 
is  exactly  on  the  sub-point  path;  these  idealizations  should  be 
replaced  if  more  precise  vehicle  altitude  and  scanner  dynamics  are 
available . 

4.1.3  Calculation  of  SOM  Map  Coordinates 

Upon  convergence  (usually  4 or  5 iterations)  of  eqn.  (13a),  with 

t*  the  converged  value,  the  integral  terms  of  the  map  projection  equations 

(1)  can  be  evaluated  from  the  series  (8a)  and  (8b) . The  second  terms 

of  equations  (1)  are  then  calculated  immediately,  given  the  values  of 

a,  Y > 3nd  6 from  Appendix  B for  the  instant  t=t*. 
s s 

4.2  Inverse  SOM  Transformation 

The  ellipsoidal  coordinates  C}i  ,X)  corresponding  to  map  coordinates 
(x,y)  are  determined  as  follows: 
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First  the  instant  t*  that  the  scan  vector  w is  colinear  with  the 

displacement  vector  6R  from  the  groundtrack  point  of 

interest  (x,y)  is  determined  via  the  iteration 
Jk+1) 


(k)  G[^] 


d£| 

dt' 


(14) 


t=t 


(k) 


where 


G(t)  = (w  X n) 


6R 


6R  = [x  - X (t)]  t + [y  - y (t)]j 

o o 

Explicit  equations  are  obtained  for  calculation  of  G(t)  and 


(15) 

(16) 


dG(t) 

dt 


by  substituting  n and  w from  Appendix  B.  The  starting  approximation 

for  (14)  is  taken  as  t^°\l  x/x  . The  converged  value  [usually  requiring 

5 

3 or  4 iterations  of  (14)]  is  t*.  The  unit  vectors  i and  J are  parallel 
to  the  X and  y axis  of  the  mapping  plane.  The  correction  for  finite  scan 
time,  eqn.  (13b),  should  be  added  to  the  converged  t*  from  (14).  The 
associated  partials  assume  an  infinite  scan  speed. 

Knowing  t*,  it  is  a simple  matter  to  determine  (4>,^)  via  the  2 


variable  successive  approximation 


9x  1 

(k) 

3x  1 

(k) 

3(^'' 

3y  1 

(k) 

(k) 

34)' 

at' 

-1 


X - X 


(k) 


(k) 


(17) 

where  and  are  substituted  into  eqns.  (1)  to  determine  x , y 

and  the  partial  derivative  matrix  is  determined  via  the  formulas 

summarized  in  §6.4.  The  iteration  (17)  has  been  found  very  well  behaved 

(usually  converging  in  4 to  5 iterations) . The  starting  estimates  are 

taken  as  4)^°^  = ip  (t*),  = X (t*),  calculated  from  formulas  of 

g g 

Appendix  A. 


11 


5.0  COMPUTATIONAL  TESTS  FOR  THE  ERTS-1  (LANDSAT)  ORBIT 
We  adopted  a nominal  orbit  with  an  inclination  of  99°,  a period  of 
103.267  minutes,  and  zero  eccentricity.  This  orbit  approximates  that 
of  the  Earth  Resources  Technology  Satellite  (ERTS-1,  also  known  as 
LANDSAT).  For  the  reference  ellipsoid,  we  adopted  the  values 
a = equatorial  radius  = 6378.165  km 
b = polar  radius  = 6356.783  km 

These  correspond  to  a flattening  of  1/298.3  or  an  eccentricity  of 
0.0818130.  In  earth-fixed  coordinates,  this  satellite  generates  the 
sub-point  trace  shown  in  the  oblique  views  of  figures  5.1a  and  5.1b. 

For  the  special  case  of  the  LANDSAT  orbit,  we  summarize  below  numerical 
results  of  the  major  options  of  the  software  documented  in  Appendix 
D. 

5.1  LANDSAT  Groundtrack  Projection 

Following  the  approach  of  4.1.1,  the  software  of  Appendix  D was 

used  to  evaluate  the  groundtrack' s projection  into  the  map  plane  as 

9 

X (t)  = f^Vcos  fdT  = At  + y A sin  (-^)  (5.1a) 

g ■'  o ^ , n P 


t 9 

y (t)  =f  Vsin  fdT  + y(t)=y  B cos  (^-^) 
^g  J ^g  O n P 


where  the  coefficients  were  found  to  be 


_ ^ P dx 

A = X = zr  f dt  = 6.504961  km/sec 

o g P ^ dt 

p 


(5.1b) 


(5.2a) 


(5.2b) 


(5,2c) 
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n 

A 

n 

B 

0 

6.50496 

1 

+.00000 

-.00000 

2 

-.00000 

916.45110 

3 

-.00000 

.00000 

4 

9.73648 

-.00000 

5 

.00000 

.00000 

6 

.00000 

-.09648 

7 

-.00000 

.00000 

8 

.00232 

-.00000 

9 

-.00000 

.00000 

As  is  evident  from  the  coefficients,  only  two  harmonic  terms  each  are 
actually  required.  This  fact  is  evident  only  after  the  coefficients 
are  determined  and  in  general  is  a function  of  the  particular  orbit. 


Appendix  D provides  Table  D1  listing  (f)  , X , x , y , for  100  equally 

8 S 8 S 

spaced  time  increments  spanning  the  orbit.  The  map  plane  projection 

(x  ,y  ) of  the  LANDSAT  groundtrack  is  the  bold  line  of  figure  5.2. 

8 8 

5.2  Example  Forward  and  Inverse  Transformations 

Using  the  LANDSAT  orbit  data  and  ref.  ellipsoid  of  §5.0  and  the 
calculation  sequence  summarized  in  §4.1,  the  following  ellipsoid  point 
(iJ>,X)  = (-0°. 16549  , -7°. 02310) 

resulted  [eqns.(l),  using  software  of  Appendix  d]  in  the  following  map 
plane  coordinates 

(x,y)  = (10099,66  km,  - 58.587  km). 


On  page  109. 
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When  substituted  into  the  Inverse  transformation  (as  discussed  in  §4.2 
and  implemented  in  Appendix  D) , the  above  (x,y)  values  result  in 
recovery  of  the  original  (i}>,X)  values  (to  within  10  . The  printouts 

of  the  "Forward  and  Inverse  Test  Cases"  of  Appendix  D supply  values  of 
intermediate  quantities  en-route  to  these  end  results. 

5.3  The  LANDSAT  Graticule 

By  simply  selecting  a sequence  of  closely  spaced  (i{i,X)  values 
along  lines  of  constant  (}>  and  X,  the  forward  transformation  equations 
of  §4.1  (as  implemented  in  Appendix  D)  result  in  the  projections  of 
the  meridians  and  parallels  into  the  map  plane.  The  resulting 
graticule  is  displayed  in  Figure  5.2.  Except  for  (j)  = + 85°,  the 
parallels  are  plotted  at  a 10°  Interval  (-80°,  -70°,  -60°,  70°,  80°) 
and  the  meridians  are  plotted  at  a 5°  interval.  The  region  displayed 
is  for  roughly  + 20°  (central  angle)  of  the  groundtrack. 

Observe,  qualitatively,  that  angles  are  well-preserved,  even  20° 
off  the  groundtrack  (meridians  and  parallels  intersect  at  right  angles) . 
However,  observe  that  slight  shape  distortions  are  evident  in  the 
departure  of  the  80°  and  85°  parallels  from  circles.  As  is  evident  in 
the  formulation  herein  (and  as  is  clear  in  the  numerical  error  analyses 
of  §5.4),  rigorous  satisfaction  of  constant  scale  and  conformality  are 
achieved  along  the  groundtrack.  These  distortions  off  the  groundtrack 
are  entirely  satisfactory  within  + 100  km  of  the  groundtrack  (the 
approximate  length  of  the  LANDSAT  scan-lines). 
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5.4  Length  Distortion  Analysis 
Adopting  the  notation 


s = arc  length  measured  along  some  line  on  the  ellipsoid 
s'  = arc  length  measured  along  the  corresponding  line  in  the  map 
plane 


then  the  basic  equations  for  length  distortion  analysis  are 

2 2 T /?  2 

/3s’.  ^ (1-e  sin  ()))  ' r .ax.  2 .3y.  -|-*j 

dS  A 


(5.3a) 


a(l-e2) 

local  length  distortion  along  a meridian  (line  of  X=constant) 

2 


(5.3b) 


= local  length  distortion  along  a parallel  (line  of  <i)=constant)  . 
The  partial  derivatives 


8x 

ax~ 

3 <() 

ax 

Iz 

Iz 

ax_ 

are  developed  in  §6.4. 

Evaluation  of  equations  (5.3)  at  various  points  along  and  near  the 
groundtrack  resulted  in  the  table  of  test  case  4 (Appendix  D) . As  is 
evident,  the  absolute  length  distortions  are  zero  along  the  groundtrack, 
worst  case  errors  of  less  than  3 parts  per  10,000  occur  at  the  outer 
fringe  of  the  sensed  region  (+  1°) . 
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6.0  DETAILED  FORMULATION  OF  THE  MAP  PROJECTION 
6.1  Preliminary  Remarks 

In  this  section,  we  present  the  derivation  of  the  map  projection 
roughly  in  the  manner  it  was  developed  originally,  (sans  the  unproductive 
blind  alleys'.)  attempting  to  expose  major  features  of  the  logical 
process  underlying  the  derivation.  This  logical  process  (and  the 
resulting  mathematics)  partitions  naturally  into  two  major  steps. 

With  reference  to  Figure  6.1,  the  first  step  is  to  project  the 
satellite  sub-point  path  (groundtrack)  from  the  reference  ellipsoid 
onto  the  mapping  plane.  This  transformation  should  be  such  that  the 
groundtrack  length  is  not  subject  to  local  length  distortions  (zero 
scale  distortion)  and  the  "shape"  of  the  groundtrack  should  be  preserved 
(curvature  constraint)  at  every  point.  These  two  objectives  lead 
Immediately  to  two  corresponding  differential  equations  which  are 
developed  and  solved  in  §6.2. 

With  reference  to  Figure  6.2,  the  second  step  is  the  projection  of 
all  sensed  points  in  a finite  region  on  the  ref.  ellipsoid,  centered  on  the 
groundtrack  (i.e.,  the  "sensed  ribbon"  on  the  earth's  surface)  into  the 
mapping  plane.  This  problem  is  approached  using  the  intuitively 
clear  notion  that  small  displacements  are  made  on  the  ellipsoid  from  a 
"locally  nearly  straight  line"  (the  groundtrack;  the  radius  of  curvature 
of  the  groundtrack ' s tangent  plane  projection  for  LAMDSAT  orbits  varies 
from  over  three  earth  radii  to  Infinity,  for  example) . This  logic  led 
us  to  consider  the  idealization  that  displacements  near  the  groundtrack 
to  nearby  points  might  be  well-approximated  by  displacements  near  the 
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equator  of  an  oblique  mercator  projection  (where  the  oblique  equator  is 

locally  tangent  to  the  groundtrack) . The  resulting  map  projection 

(based  upon  this  idealization)  is  developed  such  that  rigorous  conformality 

and  length  preservation  are  satisfied  only  along  the  groundtrack,  but 

the  approximation  is  excellent  within  several  hundred  km  of  typical 

satellite  groundtracks . Since  this  intuitively  appealing  approximation 

has  worked  out  so  well,  we  have  not  pursued  the  possibility  that  a more 

rigorous  approximation  to  conformality  and  length  preservation  may  be 

feasible  (we  did  devote  sufficient  attention  to  this  question,  however 

to  rule  out  the  existence  of  exact  conformal  solutions  for  any  cases 

other  than  trivial  cases  such  as  equatorial  orbits  or  inclined  orbits  about  a 

non-rotating  spherical  earth'.).  We  develop  the  geometry  and  equations 
for  the  local  mercator  approximations  in  §6.3.  These,  in  conjunction 
with  the  groundtrack  projection  equations  of  §6.2,  are  the  essential 
formula  of  the  desired  map  projection. 

In  §6.4,  we  develop  the  equations  necessary  to  rigorously  compute 
the  partial  derivatives 


which  are  necessary  to  do  distortion  analysis  (§5.2)  and  which  are 
required  in  the  inverse  transformation  (§4.2)  from  given  (x,y)  in  the  map 
plane  to  the  corresponding  ((}),X)  on  the  ellipsoid. 

6.2  Length  and  S hape  Preserving  Projection  of  the  Satellite  Groundtrack 
First  we  state  two  desired  constraints  which  lead  directly  to  the 
projection  of  the  sub-point  path  G from  the  ellipsoid  to  the  map  plane 
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(see  Fig.  6.1).  Observe  that  a transformation  is  desired  from  (<j>  ,X  ,t) 

g g 

along  a line  on  the  ellipsoid  to  the  corresponding  (x  ,y  ,t)  along  a 

g g 

line  in  the  map  plane.  The  line  is  generated  as  t varies  from  zero  to  an 
orbital  period.  Since  a unique  "two-to-two"  mapping  is  desired,  it  is 
clear  that  two  independent  constraints  are  necessary  and  sufficient 
to  establish  the  desired  transformation. 

Constant  Scale  Constraint 


s = arc  length  measured  along  the  satellite  sub-point  trace 
(groundtrack) 


= / V(t)  dt  (6.2) 

0 

s*^  = arc  length  measured  along  the  projection  of  the  groundtrack  onto 
the  map  plane 


t / dx  „ dy  „ 

= I dr. 


We  require,  for  no  scale  distortion  along  the  groundtrack,  that  s=s'  for 
all  values  of  t.  By  inspection  of  equations  (6.2)  and  (6.3),  this  is 


possible  only  if 
« dx  „ 


(6.4) 


is  satisfied  at  every  point  (x  ,y  ,t)  along  the  map  plane  projection  of 

g g 

the  groundtrack. 

Curvature  Constraint 

To  preserve  the  groundtrack  "shape",  we  require  that  the  radius  of 
curvature  of  the  groundtrack  projection  (in  the  map  plane)  equal  at 
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every  point  the  instantaneous  radius  of  curvature  (p)  o£  the  projection 


of  the  groundtrack  into  the 

sub-point  (see  Figure  6.3). 

2 2 
d X „ d y „ T 
(_^)2  + 

ds  ds  P 


S olution  of  eqns.  (6.4)  and 


the  independent  variable  of 

(6.2)  that 


plane  tangent  to  the  ellipsoid  at  the 
This  requires  that 


(6.5) 

(6.5)  for  the  Groundtrack  Projection.  To  change 
eqn.(6.4)  from  t to  s,  observe  from  eqn. 


ds 

^ V 

dt 

dx 

dx 

= —3. 

dt 

ds 

dy 

dy 

= —3 

dt 

ds 

ds 


dx 


= V 


ds 


V 


ds 


so  that  the  length  constraint  can  be  written  as  simply 

dx  „ dy  „ 

(-^)^  + (^)  = 1 

ds  ds 


(6.6a) 

(6.6b) 

(6.6c) 

(6.6d) 


Thus  we  require  a simultaneous  solution  for  x and  y satisfying  equations 

8 8 

(6.5)  and  (6.6d).  By  inspection  of  (6.6d)  it  is  clear  that  the  solution 
for  the  derivatives  can  be  taken  as  simply 


dx 

= cos  f 
ds 


dy 


ds 


= sin  f 


(6.7a)* 


(6.7b)* 


Strictly  speaking  + signs  should  accompany  these  equations;  numerical 
experiments  support  the  conclusion  that  the  positive  signs  adonted  are 
correct;  note  that  the  sign  of  the  right  hand  side  of  eqn.  (6.7) 
ran  change  with  the  variation  of  p (see  App>-ndix  B)  . 
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where  f is  an  appropriate  function  to  guarantee  satisfaction  of  the 
curvature  constraints  eqn  (6.5).  Substitution  of  (6.7)  into  (6,5)  leads 
immediately  to  the  conclusion  that  the  first  derivative  of  f must  satisfy 


= i 

ds  p 

Making  use  of  eqn.  (6.6a)  eqns.  (6.7)  can  be  written  as  the 


(6.7c)* 


differential  equations 


dx 

— ^ = Vcos  f 
dt 

(6.8a) 

dy 

— ^ = Vsin  f 
dt 

(6.8b) 

ii  = I 

dt  p 

(6.8c) 

The  solution  of  which  is  indicated  formally  as 


f(t)  =/  dx 

" pm) 

0 

(6.9a) 

t 

X (t)  = 0 + / V(x)  cos  f (t)  dx 
^ » t 

(6.9b) 

y (t)  = y (o)  + / V(t)  sin  f (x)  dx 
g g ' 

(6.9c) 

where  V(t)  and  p(t)  are  calculated  as  established  in  Appendices  A and  B. 

It  is  obvious  that  the  scale  and  curvature  constraints,  in  the  form 
of  equations  (6.4)  and  (6.5)  are  satisfied  by  equation  (6.8)  and  therefore 
equations  (6.9) 

As  is  discussed  in  §4.4.1,  the  evaluation  of  the  integrals  in 
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eqn.  (6.9)  (if  required  for  many  t values)  is  greatly  facilitated  by 
replacing  equations  (6-9)  by  their  fourier  series.  The  integration  of 
(6.9)  can  be  done  sequential ty  or  simultaneously.  Since  the  right  hand 
sides  of  (6.9b)  and  (6.9c)  contain  f(t),  f(t)  can  be  integrated  first 
from  (6.9a)  and  then  the  (6,9b)  and  (6.9c)  integrations  can  be  carried 
out.  Alternatively,  we  have  found  it  much  more  convenient  to  do  the 
integrations  by  using  the  Runge-Kutta  algorithm  (Appendix  C)  to  integrate 
all  three  of  equations  (6.8)  simultaneously,  this  is  the  method 
recommended,  [if  a large  number  of  ((Ji,^)  points  need  to  be 
transformed  to  the  corresponding  (x,y),  the  calculation 
of  the  Fourier  series,  discussed  in  §4.1.1,  is  strongly  recommended]. 

The  software  of  Appendix  D implements  the  Runge-Kutta  integration 
procedure  to  determine  the  coefficients  of  the  Fourier  series 

expansions  for  f(t),  x (t) , and  y (t). 

8 S 

6.3  Local  Oblique  Mercator  Projection  of  the  Sensed  Region  from  the 
Ellipsoid  to  the  Map  Plane 

Since  the  scanned  (sensed)  points  represent  small  displacements  off 
the  satellite  groundtrack  (which  is  typically  relatively  straight  on  a 
local  scale,  the  LANDSAT  orbit’s  radius  of  curvature  component  in  the 
tangent  plane  varies  from  several  earth  radii  to  Infinity),  one  is 
motivated  to  consider  approximations  to  account  for  the  small  displacements 
from  the  vigorously  projected  groundtrack.  We  will  now  discuss  local 
approximations  which,  together  with  §6.2  complete  the  map  projection. 

6.3.1  Scan  Time  Determination 

Peculiar  to  the  dynamic  map  projection  under  consideration  is  the 
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necessity  to  associate  a particular  time  with  the  projection  of  (ijj,!) 
into  the  corresponding  (x,y).  In  particular,  this  time  is  denoted  t*, 
it  is  the  instant  that  the  satellite  scan  vector  passed  over  the  earth- 
fixed  point  (see  Fig.  6.4d).  Assuming  infinite  scan  rate,  this 

instant  is  characterized  by  the  fact  that  the  point  ((p,\)  and  the 


satellite  sub-point  (iji  ) must  lie  in  the  plane  determined  by  the  vector 

8 8 

^ A 

n (normal  to  the  ellipsoid)  and  the  vector  w (the  scan  vector,  nomially 

normal  to  the  orbit  as  seen  in  inertial  space) , This  condition  leads 
immediately  to  the  constraint: 

F(t*)  = [(w  X = 0 (6.10) 


where  the  time  varying  vectors 


n = n i + n i+n  k 
— X — y ^ z — 


w = w i + w i+w  k 
— X — y ^ z — 


(6.11) 


are  computed  according  to  formulas  given  in  Appendices  A and  B and 

AR=(R  -R  )i+(R  -R  )j+(R  -R  )k 
— X xg  - y y8  z 28“ 


(6.13) 


and 


= the  vector  from  (tj)  .X  ) to  (ifi.X) 

8 8 


R = X - H , R = y - H , R = z - H 
xg  x^  yg  y zg  z 

(x,y,z)  earth-fixed  components  of  the  satellite  position  vectors, 
calculated  from  eqn.  (A16) . 

(H  ,H  ,H  ) earth-fixed  components  of  the  satellite  height  vector. 
X y z 

calculated  from  eqns.  (B.2). 

^ 2,  2 2^  ^ , 2 . 2^  , -is  ^ - 

= a (a  cos  + b sin  ) cost})  cosA 


2 2 2 2^ 

Ry  = a (a  cos(^  + b sin  (Ji  ) "^  cosiji  slnX 


R = b (a  cos  i}i  + b sin  iji  ) sln(Ji 
z 


(6.1A) 
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Figure  6.4a  Ellipsoidal  Georaetry 


K 


Satellite  Groundtrack 


Figure  6. A B Ellipsoid  Geometry:  Tangent  Plane  Triads  {s  e,  n 


track 

track 


Groundtrack 
on  the  Ellipsoid 


R = C s Radius  of  Curvature 
c 

a 2 

- (x  - H ) c.i^i  c\  + ('/-//  JciJsX  + (-)('?.-  H)s<ti 
X R R ■ >1  R R h Z 


a = tan  ( — 

c 


Figure  6.4C  Ellipsoidal  Geonetry;  Ellipsoid  Curvature  in 
n.  c Plane 


When  earth-fixed  components  of  the  three  vectors  are  substituted  into 


(6.15a) 


eqn.  (6.10),  the  function  we  seek  to  find  the  zero  of  is 

F(t*)  = [(R  + //  - x)(w  n - w n ) 

X X y z z y 

+(R  + H -w)(wn  -wn) 

y y z X X z 

+ (R  + H - 3) (w  n -wn)]  j=0 

z z X y y X -'t=t* 

and  the  time  derivative  F(t)  follows  immediately  upon  taking  the 

derivative  of  (6.15a).  Newton's  method  is  used,  as  is  discussed  in 

§4.1.1  to  find  t*  such  that  (6.15a)  is  satisfied.  Usually  convergence  is 

achieved  within  5 iterations. 


Implicit  in  eqn,  (6.15a)  is  the  assumption  that  the  scan  is  instantaneous 
(all  points  under  the  scan  vector  are  imaged  simultaneously) . A 
rigorous  calculation  of  the  time  varying  location  of  the  imaged  point 
((J>,X)  requires  a rigorous  model  for  the  scanner  motion  and  the  attitude 
motion.  For  our  purposes  here,  it  is  sufficient  to  correct  linearly 
for  finite  scan  time  by  adding  the  correction 

At  = T(— ^)  (6.15b) 

E 

max 

to  the  root  t*  satisfying  eqn.  (6.15a)  , where 

T = '2  the  scan  period 

= 18.355  ms  for  the  T.ANDSAT  scanner 

e = the  maximum  scanner  beam  deflection  away  from  local  vertical 
max 

= 5.78  deg  for  the  LANDSAT  scanner 
e = instantaneous  scanner  deflection  angle  away  from  local  vertical 
(see  Fig . 6 . 4D) 

-1  -S-H 

= cos  ( — g^)  • sign  (AR-£)  (6.15c) 


where  S = R - r 
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Fif’ure  6.5  Map  Plane  Geonetry:  Local  Mercator  Approximation 
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6.3.2  Local  Oblique  Mercator  Approxlina t ion 

With  particular  reference  to  Figures  6.4d  and  6.5.  the  equations 
relating  a given  (ip,\)  on  the  ellipsoid  to  the  corresponding  (x.y)  in 
the  map  plane  are 


z-*"*  ^ r, 

X =f  Vcos  f dt  + R ln[tan(—  + — )]  ^ 

\ c ^ \ 2 cose 

0 s 

t* 

y = y„(°)  + / '^sin  f dt  + R ln[tan(7  + 


(6.16) 


siny 


s 


cos9 


where  the  integral  terms  are  the  groundtrack  coordinates  (x  ,y  ),  developed 

g g 

in  56.2  and  from  the  geometry  of  figures  6.4d  and  6.5 


0 = - sin  ^(w't)  = - sin  ^(w  t + w t ) = scan  angle 

s xxzz  “ 


7T  -1 

Y = T + 9 + tan  (-7^) 

s 2 s X 

g 


where 


w = scan  vector,  projection  into  t_,  £ plane  or 


(6.17a) 

(6.17b) 


w . t 


w = (- 


) t + ( -)c 

u — u — 


= (w'  • t)^  + (w'  • c)^ 


(6,i8a) 

(6.18h) 


and 


w = orbit  normal,  in  Inertial  frame 

= ^x  _r  = constant,  for  unperturbed  orbit(r  = spacecraft  (6.18c) 
position  vector)* 

R = ellipsoid  radius  of  curvature  in  the  _n,  c^  plane 

2 

= (x  - H )cos4>  cosX  + (u  - H jcosj)  sinX  + ^ (^  - H .'sinij)  (6.10) 
X z a ''  -j  SL  R .2  z g 


•k 

See  Appendix  A 
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, c . AR  , c (R  -R  ) + c (R  - R ) + c (R  -R  ) 

-1  , -1  , X X xg  y y yg  z z zg 

a = tan  ( — r ) = tan  ( = ^ ^ ^ 

K K 

C C 

A 

= the  angle  between  ii  and  the  vector  from  the  cross-track  center 
of  curvature  to  the  orthogonal  projection  of  point  (4>,A)  onto 
c^  (see  Figure  6.4c). 

(c  ,c  ,c  ) = earth-fixed  components  of  the  cross  track  vector  c, 

X y z — 

calculated  from  eqn.  (B.8)  of  Appendix  B. 

Notice  that  the  distance  R ln[tan(^  + ^) ] , measured  normal  to  the 

c 4 z 

groundtrack  (Figure  6.5  and  Eqns.  6.16),  is  identical  to  the  classical 
formula  for  "y"  of  the  transverse  mercator  projection  (see  References 
2 and  3).  Thus,  to  the  extent  that  the  groundtrack  approximates  a 
great  circle,  the  map  projection  approximates  an  oblique  mercator 
projection.  However,  the  degree  to  which  the  classical  oblique  mercator 
map  pr  )jection  is  approximated  is  not  terribly  important,  in  as  much 
as  the  present  projection  is  motivated  by  the  fact  that  none  of  the 
classical  projections  fulfill  the  objectives  being  pursued  here  (primarily, 
a continuous  distortion-free  mapping  of  the  entire  sensed  region  of 
typical  earth-scanning  satellites).  As  is  noted  in  §5.3,  the  length  and 


angle  distortions  are  sufficiently  small  to  compare  favorably  with  the 
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corresponding  errors  for  the  classical  transverse  mercator  projection. 
6.4  SOM  Partial  Derivatives  (for  Distortion  Analysis  and  Inverse 
Transformations) 

We  summarize  here  the  equations  in  back-substitution  form  for 
computation  of  the  partial  derivatives 

(6.20) 


3x 

3x 

3(}) 

3X 

-e- 

3X 

which  are  necessary  in  the  calculation  of  54.2,  5.2  and  5.4, 

To  compact  the  notation,  we  will  make  use  of  the  symbol 
whenever  the  identical  equation  for  X results  by  simply  replacing  <()  by 
X.  These  equations  follow  by  simply  applying  chain  partial  differentiation 
to  the  equations  of  §6.2  and  6.3.  While  references  are  frequently  made 
to  the  equations  differentiated  to  obtain  these  results,  the  present 
treatment  is  more  nearly  an  ''annotated  summary"  rather  than  a rigorous 
derivation . 


Taking  partial  derivatives  of  equations  (6.16)  yield 

(ji-'-X 


9x  ^^g  3D  ^^s  _ 

8^  = a/ ^ "°"^S  - """''s' 


3y  ^^g  3D  . ^"s 

^ = 3^  ^ ° """"s’ 


(6.21) 


b-^x 


where 


t* 

X H I Vcos  f dt 
8 o 

t* 

y =y(o)+f  Vsinfdt 
8 8 

D = ln[tan(-2  + ^)]/cos0^ 


(6.22) 


(6.23) 
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8x  9y  3y 

The  partials  (jf  , j/'  > jf 

are  determined  as  follows: 


; <|)->-X)  needed  to  calculate  eqns.  (6.21) 


X ,y  partials 
-g  S 


from  eqns.  (6.22) 

— ^ = V(t*)cos[f (t*)]  d)->X 

3(|>  3t*  34)  ^ ;cosLivc  -'J  a 

3t*  r T 3t^ 

“4  = irf  ^ = V(t-^)sin[f(t*)]  — , 


(6.24) 


D Partials 


3D  1 


3(!)  cose_'\  34.  4'^  34 


„ 2, IT  a, 

3a  ^c  ^4  + 2^ 


r,  /TT  , a, 

2 tan(-2|-  + j) 


sinG 

- R ln[tan(-?  + ■^)  ] ^ 
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‘ (cose 

s 


, <))^x 


(6.25) 


where 


DR 

3R 

34 

3R  3X  3R 

c 

= [— 
‘-34 

+ _c  _S.  + _Ji 

d<t> 

3t* 

3X  3t*  3t* 

g 

g 

3R 

2 3R  : 

-eiJ 

c4  sX 
g 

,3  X 1 

g ^2  3t*  : 

3R 

c 

= - R 

s4  cX 

- R s4  sX  + 

3()) 

g 

X 

g 

g y g g 

3R 

c 

= - R 

c4  sX  + R c4  cX 

3 X 

g 

X 

g 

g y g g 

3a  _ 

1 

3AR 

f r»  r — 

34 

■r2.p 

C 

(AR  • 

;^2  ^‘^cL34  • - 

-Y  R cifi 
2 z g 


3c 


9R 


(6.26) 

(6.26) 

(6.23) 


+ ^ (6.29) 


3AR 

?r=  ("r  ^2’  ''3) 


(6.30) 
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with 


9N 

= -5— ■ c<J’cA  - Ns4'cA 
1 3<p 

9N,, 

T_  = -r—  cificX  - N siJisX 
L d<p 

h ■ 4 <1^  »♦  + 

a 

then 

2 2 2 2 2 -^5 

N=a  (aciji  + bs(j)) 

9N  „3  (a^  - b^)  . ,,,, 

-r^  = - N sin(2(|)) 

2a^ 

9AR 
9X“  = 

with 

= - N C())sA 

= N c(t)cX 


then 


9^ 
9 (p 

9in 

"9^ 


9 

Dn 

A A 

X t ) = 

X t + n X 

dn 

9b 

Dn 

g 

— s + 



gt*-'  9b 

9t* 

9X 

8 

tv- 


dt* 


(ti->-X 


(T^,T2,T^) 


(6.31a) 

(6.31b) 

(6.31c) 

(6.32) 

(6.33) 

(6.34) 

(6.35a) 

(6.35b) 


V 


dV  9t*  . , ,, 

2 dt*  dt*^  9b  (6.36) 


(6.37) 


(6.38) 


with 


T,  = -sb  cX  (6.39a) 

1 S g 
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(6,39b) 


-s<t>  sA 
8 8 


T = c(t> 

3 8 


then 


w " (5^.82,52) 


with 


S,  = - ct])  sA 
1 8 8 


S.  = c(j)  cA 
2 g g 


S3  = 0 


then 


3v_  _ 1 r!^ 

9t*  V ^ dt 


From  eqn.  (6.17a) 

3e 


®d^R 


dt 


t=t* 


" ■ 1 - (;  . t) 

From  Eqn.  (6.18) 


8w 


2 34* 


. 3^ 

1 + 


^ 3t* 

34  3 1 3 (J> 

See  eqn.  (6.62)  for 


, 4)->-A 


3w 


3t* 


3^  dR 

_ = _i  [i  _I]=  i [ 

3(|)  34)  dt*-'  V 


®d^R 


dt* 


i IV  1^ 

V 3t*  dt*  ■‘dip  ’ 


y Partials 
s 


From  eqn.  (6.17b) 

3y  30  . 

S S , 3 

34)  ° 34'  3't> 


(jS)] 

8 


(6.39c) 

(6.40) 

(6.41a) 

(6.42b) 

(6.43c) 

(6.44) 


(6.45) 


(6.46) 


(6.47) 
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or 


9y  90  T dt  „ -1  dx  d^y  dy  d^x 

_Ji  = -^  + (— ^)^]  {—^  ^ - 

d(^  9$  ^dt  ' Mt  Mt*  , ,2  dt*  , J Stf 

d t d t * 


3t* 


(6.48) 


with 


d^x 


dt* 


2 " “d^  (f(t*))  - V(t*)  sin  (f(t*)) 


(6.49a) 


dt* 


2 ^ (f(t*))  + V(t*)  cos  (f(t*))  — 


(6.49b) 


Fron!  the  t*  derivative  of  eqn.  (6,10) 


9t* 

9<fi 


.9F  .-1  9F 
^9t*^  '9((>  ’ 


(6.50) 


when 

F(t*)  = [w  X ti) -AR] ; tp  H w X n 

= _t p • AR  ; , AR  s (AR  • _t ) ^ + (AR  . c)  c_ 

then 

F(t*)  = (t^  . t ) (aR  • t)  + (t  . £ ) (aR.  • £) 

and 

9(AR)  , , . 9(AR)  . . 

— = (— - — • t)(t  • tp)  + (-— — • c)(c  • tp)  , (j)->X  (6.51) 

3(()  3())  — — — — — 

9(AR)  . 9£  . . . 9t  . . 9tp 

• 1)  + (AR  • 37^)>  (1  • tp)  + (AR  • t){(^  • tp)  +(t  • ^)} 

9 (AR)  . 9c  „ . .31  . 9£p 

+ ((— ^ ■ S)  + (AR-9^)}(c  • tp)  + (AR-  c){(g-^  • tp)  + (£  • g-^)) 
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(6.53) 


9 1 9ri 


X w + n X 


w = r X r 


(normal  to  orbit  in  the  inertial  frame)  (6.54) 


n' , n e ' , n e ' ' , 

w=  w + m xw  , to  = 0k  = ang.  vel.  of  earth 


% = r X (=0,  for  un-perturbed  elliptical  orbits) 


(6.55) 


(6.56) 


Given  inertial  components  of  w'  and  w' , the  earth-fixed  components  are 


obtained  via  the  transformation 


{w'}  = [E-(9)]{w'} 

~ e j — n 

• • 

= [E  (0)]{%} 
e 3 n 

Te  = c0  s0  0 

'-3^'''^-’  -s0  c0  0 

0 0 1 


0 = 0 +0)  (t-t  ) to  = rotational  rate  of  the  earth 

o e o e 


9 w ^ „ 3£  3 1 

TTTC  = At  + Bc  + c - — T + D — — - 

9t*  — — 3t*  3t' 


(6.58) 


(6.59) 


(6.60) 


(6.61) 


(6.62) 


A - ^ [w'  _ t.  + w'  • 2 — {w'  • ^[w'  •_t  + w-  W'£  L£’‘£  + w'  • 


3 t w' . t 


(6.64)  B = — [w'-c  + w' -^r — j 
u — 3t 


c = 1;- 


— W C 
U — — 


^ w'  . . . , , ^ 1 . 9£ 

— r -■  ^ • 3i^^^ 

u 

(6.65) 


1 . . 
D = — w • t 


2 . ^ 2 , - 2 

u = (w'-t)  + (w’-c) 


(6.66) 


(6.67) 
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7,0  CONCLUDING  REMARKS 


This  report  documents  a fairly  general  formulation  and  implementation 

of  the  Space  Oblique  Mercator  map  projection  which  embody  the  attractive 

features  forecast  by  Colvocoresses  when  he  first  conceived  of  this 

projection.  During  the  middle  stages  of  this  development,  the  authors 

became  aware  that  a parallel  effort  was  being  carried  out  by  John 
(4) 

Snyder  . Based  upon  recent  communications,  it  is  clear  that  Snyder 
has  accomplished  an  approximately  equal  feat,  but  using  a different 
approach. * 

Snyder's  insightful  formulation  is  much  more  compact  than  the 
present  developments,  although  the  one-time  character  of  many  of  the 
calculations  would  appear  to  diminish  the  gap  in  computational 
efficiency. 

The  question  of  "which  formulation  should  te  used"  may  boil  down 
to  the  i„sue  of  "how  much  flexibility  is  required  vis-a-vis  the 
nominal  orbit".  Should  circular  orbits  be  exclusively  desired,  then 
Snyder's  formulation  may  prove  preferable.  Should  a state-of-the 
art  orbit  Integration  program  be  used  to  define  the  nominal  orbit, 

then,  withcut  question,  the  present  formulation  is  applicable  whereas 
Snyder's  is  not,  Another  important  unresolved  issue  is  how  important 

•k 

Snyder  has  not  rigorously  enforced  the  curvature  constraint  (that  the 
map  plane  projection  of  the  groundtrack  have  the  same  radius  of 
curvature  as  the  tangent  plane  projection  of  the  groundtrack),  but 
his  solution  approximately  satisfies  this  constraint  nonetheless, 

Snyder's  results  do  not  allow  as  much  generality  with  regard  to 
selection  of  the  nominal  orbit.  Whereas  the  present  formulation 
permits  routine  use  with  either  analytical  or  numerically  integrated 
nominal  orbits,  Snyder's  formulations  were  designed  for  circular 
orbits  (although  we  understand  that  he  Is  attempting  to  modify  his 
formulation  to  account  for  non-circular  orbits). 
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rigorous  enforcement  of  the  curvature  constraint  will  prove  for  orbits 
other  than  the  LANDSAT  near-polar,  near-circular,  cases  studied  thus 
far. 

In  any  event,  we  believe  both  Snyder's  formulation  and  the  present 
formulation  are  useful  contributions;  both  will  likely  be  employed 
in  future  utilizations  of  LANDSAT  and  similar  imagry. 

As  a final  remark,  we  note  that  the  considerable  logical  processes, 
algebra,  calculus,  and  computer  programming  underlying  this  work, 
coupled  with  the  usual  editorial/typographical  headaches  leave  a finite 
probability  that  significant  errors  have  escaped  the  attention  of 
several  proof-readers.  We  sincerely  hope  that  any  remaining  errors 
prove  of  no  conceptual  or  practical  consequence,  and  we  solicit  the 
readers  communication  of  all  errors  or  suspected  errors.  It  is 
anticipated  that  computational  studies  and  further  analysis  will  lead 
to  addendums  to  this  work  in  which  all  kno\vm  errors  will  be  corrected. 
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APPENDIX  A 


ORBIT  INTEGRATION  AND  TR.\NSFORMATIONS 

A.l  Comments 

Here  we  consider  two  levels  of  generality  for  orbit  calculation, 

viz  : 

(1)  The  case  of  arbitrary,  but  given,  smooth  perturbation,  (non- 
2 body  effects)  and 

(2)  Tbe  Keplerian  2-body  orbit  (including  all  possible  species  of 
ellipital  and  circular  orbits). 

The  first  class  of  orbits  clearly  Includes  the  latter  as  an 
obvious  special  case  of  zero  perturbations.  However,  the  well  known 
analytical  solution  for  the  Keplerian  case  is  so  efficient  (relative  to 
numerically  integrated  orbits)  that  we  treat  it  separately. 

A. 2 Equations  of  Motion,  Inertial  Coordinates 

We  observe  first  that  the  differential  equations  of  motion  for  the 
general  case  have  the  form 


X = -GM  X/r^ 

+ 

(X 

perturbation 

acce lerat ion) 

Y = -GM  Y/r^ 

+ 

(Y 

perturbation 

acceleration) 

7 = -GM  Z/r^ 

+ 

(Z 

perturbation 

acceleration) 

whore 

r = XI_  + Yd_  + 7X  = satellite  rectangiilar  coordinates  in  a non-rotating 
earth-centered  frame,  with  7 along  the  polar  axis. 

2 2 2 2 

r = X + Y + 7 n (A. 2) 

km 

GM  = Earth's  gravitation-mass  constant  = 398601.2  ^^^2 

• • « 

Given  initial  position  and  velocity  coordinates  (X.Y,7,X,Y,Z) 

o o o o o o 
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at  time  t^,  numerical  methods  such  as  Runge-Kutta  (Appendix  C)  can  be 
employed  to  Integrate  equations  (A.l)  to  determine  instantaneous 
position  and  velocity  (X, Y , Z ,X,Y ,Z)  at  an  arbitrary  given  time  t, 

A. 3 Satellite  Motion;  Analytical  Solution  for  Inertial  Rectangular 
Coordinates  (Keplerlan  Special  Case) 

The  following  equations,  in  order  of  solution,  determine  the  inertial 


coordinates  (X,Y ,Z ,X,Y ,Z)  at  time  t,  given 

the  same  quantities  (X  , 

o 

. 

at  time  t : 

0 

constants : 

y s CM  = 398601,2  km^/  sec^ 

2 2 2 2 

r = X + Y + Z 

O 0 O 0 

(A. 3) 

D = XX  + YY  +ZZ 

O 0 0 O O 0 o 

(A. 4) 

2 * 2 '2  ’2 

V = X + Y + Z 
o 0 0 0 

(A.  5) 

- = 2/r  - V^/y 

a 0 0 

(A. 6) 

c = 1 - r /a 

0 0 

(A. 7) 

Solve  for  the  change  in  eccentric  anomoly  E 

; (using  Newton's  method) 

from 

- D 

/y(t-t  )a  ' = E - (1-r  /a)  sinf.  + 

0 0 /y  a 

( 1-cos E) 

(A. 8) 

then 

f = 1 - a(l  - cosE)/r 

0 

(A.  9) 

n / O A A \ 

g = (t-t^)  - a (E-  sinE)/y^ 

(A. 10) 

r = a(l  - c cosF.)  - D (— )sin  E 

0 o y 

(A. 11) 

-1 

f = -(rr  ) ya  sinE 

0 

(A. 12) 

g=l-a(l  - cosK)/r 

(A. 13) 
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and  finally 


r -N 

r ->v 

C' 

x(t) 

X 

X 

0 

0 

Y(t) 

Y 

0 

0 

Z(t) 

Z 

1 

J 

r "I 

r 

r-  > 

X(t) 

X 

X 

0 

o 

\ Y(t) 

>=  f < 

Y 

) + g< 

o 

f oy 

0 

Z(t) 

Z 

Z 

L °J 

0 

^ J 

(A. 14) 


(A, 15) 


A.4  Satellite  Motion:  Transformation  to  Earth-Fixed  Rectangular  Coordinates 
Regardless  of  whether  the  satellite  orbit  is  calculated  via  the 
analytical  solution  of  5A.3  or  whether  perturbations  are  considered  and 
equations  (A.l)  are  integrated  numerically,  the  same  transformations  must 
be  applied  to  the  inertial  position,  velocity,  and  acceleration  coordinates 
{X,Y,Z:X,Y,Z;X,Y,Z} 

to  obtain  the  analogous  coordinates 

. . . •' 

{:c,y,  z -x  ,y,  s ) 

with  respect  to  earth-fixed  axes.  The  transformations  are  compactly 
written  in  matrix  form  as 


X 


cos9  sin9  0 

-sin9  cos9  0 


L 0 


0 


y ? = 


~cos9  sin9  0 
-sin9  cos9  0 

n 0 


(A. 19) 


(A. 17) 
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(A. 18) 


cos6 

sin0 

0 

> = 

-sin0 

COS0 

0 

LJ 

_ 0 

0 

1 _ 

where 

6 = 9 +0)  (t-t  ) 

o e o 


(A. 19) 


= sidereal  time  of  greenwich 

= counter  clockwise  rotation  about  the  Z=z  polar  axis 
0)^  = angular  velocity  of  the  earth 
A. 5 Satellite  Motion:  Transformation  to  Earth-Fixed  Ellipsoidal 
Coordinates 


The  earth-fixed  ellipsoidal  coordinates 
can  be  calculated  from  the  earth-fixed  rectangular  coordinates  {^;y, 


from  the  following  equations: 


Polar  Distance  and  Derivatives 


r = X X + -iv)  / r 
P " P 

..  ..  " , p . o . 2 

r = (x  X + y y + x^  + - r ) Ir 

P ■ P P 

Longitude  (X)  and  Derivatives 

X = tan  ^ f'it/x) 

X = (x  ‘ - ; x)lv^ 

P 

• 2 -22  . ' • 1 ‘ 2 3 

X = (x  I - V X + X - y )lr  - 2 ix  / - < x r /v 

P * ■ P P 


(A, 20a) 


(A. 20b) 
(A.  20c) 


(A. 21a) 
(A.2ib) 
(A. 21c) 


Lattitude  ($)  and  Height  (H)  and  Their  Derivatives 

(p  and  H are  determined  via  the  Newton  iteration 


+ D~\k)/ 


JkJ- 


r -r  (k) 
- P P 


(A. 22a) 
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(A. 22b) 


where  the 
spherical 


2 2 
■■  2r  h dK 

z + ^ [(-2)  N + H)sin(i)  - 2^  cos(f 
a a 


2 2 

2 2 
a d(}/ 


sin 


p]-2$Hcoe 


2 

r^  + ({)^[  (N+H)cos({>  + 2-^  sin(i)  - cos())]+  2'r'flsinii) 

d(j) 


(A. 22c) 

iteration  (A. 22a)  normally  converges  in  3 iterations,  using  the 
earth  starting  approximations 


(0) 


and 


D = D(!j),H)  = 


r> 

r)  ’ 

(0) 

= 

r - 

6378  km 

-• 

3a 

d<p 

dz 

3H 

2 

a 

• A 

d^ 

,2 

+ (— 2 N+H)cos<})] 
a 

[sincji] 

3r 

3r 

rdN 

J-dlj) 

3((> 

COS(J)  - 

(N+H)  sin4>] 

[cos!}>] 

(A. 23) 


D(k)  = H^^^)  = eqn(A.23)  evaluated  with  (A. 24) 


N - a (a  cos  (ji  + b sin  if>)  ^ = radius  of  curvature  in  the  prime 


vertical  plane 

dN  (a^-b^)  ..3  . 

_ = _ ^ sinzd) 

dt  2a 

^ ^ [3N^  ^ sin2({.  + 2N^  cos2({)] 

d<|,2  2a^ 


(A. 25a) 
(A. 25b) 

(A. 25c) 


n = (— 2 N + H)sin(}) 


r = (N  + H)cosij 
P 


(A. 26) 

(A. 27) 


= eqn(A.26)  = evaluated  with 
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(k)  , , (k)  (k) 

= eqn  (A, 27)  evaluated  with  ,H  , 

In  several  instances,  the  orbit  normal  unit  vector  w is  required  in 
carrying  out  calculations  needed  in  the  text  of  this  report.  The  scan 


vector  can  be  determined  by  the  cross  product 


w'=  r X f /Ir  X f|=  (h  /h)  I + (h  /h)J  + (h  /h)K 
- — X — y ~ z — 


(A. 28) 


where 


h = (YZ  ZY)  h = (ZX  - XZ)  , h = (XY  - YX) 
X y ’ z 

2 2 2 2 
h =h  +h  +h 
X y z 


(A. 29) 


(A, 30) 


The  components  (A.  29)  are  constant  for  Keplerian  orbits,  but  m.ust  be 
calculated  instantaneously  in  the  presence  of  perturbations,  The 

A 

rotating  earth-fixed  components  of  w'are  determined  via  the  transformation 


cos(Ji  sin(})  0 


-sin(J)  coscj)  0 


0 0 


(A. 31) 


and,  when  necessary,  the  tangent  (t)  crosstrack  (c)  and  normal  (n) 


components  of  w are  determined  by 


n n 

X y 


(A. 32) 


where  the  £,  £ vectors  are  defined  (and  expression  for  computing  their 


components  are  given)  in  Appendix  B(eqn  B.8). 
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APPENDIX  B 


SUB-POINT  GEOMETRY  AND  MOTION 

B.l  The  Sub-point  Vector  R and  its  Derivatives 

1 ZS 

With  reference  to  Figure  6.4b,  it  is  clear  that  the  sub-point  vector 
R and  its  earth-fixed  derivatives  are 


R = r ~ Vl  = (x  - H ) X + (y  - H + (z  - E )]s. 

_g__  y X 

^R  = = (x  - H )±  -r  (y  - H )2_  + (z  - E )k 

g__  ^ ^ a y ^ 


(B.la) 

(B.lb) 


®R  = ®r  - ®H  = (a 
-8  - - 


H )i  + (y  - H )l+  (z  - H )k 
X — 3 ^ 3 — 


(B.C) 


The  left  superscript  e denotes  that  the  derivatives  of  the  groundtrack 
(sub-point)  vector  ^ are  taken  with  respect  to  earth-fixed  axes.  The 
earth-fixed  satellite  coordinates  (,x,y ,z;x,y ,z;x,y ,z)  are  available 
frora  Appendix  A,  eqns  (A. 16),  (A. 17)  and  (A. 18).  The  11  vector  and  its 
derivatives  follow  from  the  geometry  of  Fig.  6.4a  as 


H = H cosdcosX  = HcticA 

X r _ r 

H = Hc(J)s\ 

y 

H = Hs(> 


= HccjjcX-  H(})Si{)cX  - HXciJjsX 
= HciJsX  - H^scjisX  + HXc(()cX 
H = Hs<i)  + H4ic^ 

3 

H = Hc^cX  - Hi}iS())cX  - HXs<j)sX  - H(i  +X“)c(j)cX 
X 

- 2H^s<j)cX  - 2HXcijisX  + 2Hi})Xs<}>sX 

= H c^sX  - H(|>s<j)sX  + HXc(}isX  - + X^)ci}isX 

- 2H^s<jisX  + 2HXci{)cX  - 2H^Xs!t)cX 


E = Hs({i  + H(j)C(t)  + 2HiJc(J  - Hij  s()) 

3 


(B.2a) 

(B.lb) 

(B.2c) 

(B.3a) 

(B.3b) 

(B.3c) 

(B.4a) 

(B.4b) 

(B.4c) 
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B.2  Tangent,  Crosstrack  and  Normal  Vectors 


Referring  to  figure  6.4b,  the  unit  vector  tangent  to  the  sub-point 


path  generated  in  earth-fixed  axes  is  clearly  , 

..  . X - k i)  ~ k 3-1/ 

t = V i + 1 + ^B.5) 

Referring  to  Figure  6.4a  the  unit  vector  nonnal  to  the  ellipsoid  is 
clearly 

n = (cosdi  cosX  )i  + (cosi  sinA  )j  + (sin(|)  )k  (B,6) 

- g g-  gg^  g- 


The  crosstrack,  unit  vector  is  defined  according  to  the  right  hand  rule 
as 

A A 

^=n^x^  (B,7) 

In  summary 

^ • • ■ a a ■ 


r 

t 

X 

t 

y 

t 

Z 

V 

r-  X-ff 

y-B 

r-  — ; 

V 

f ■A 
i 

r 

c 

X 

c 

y 

C 

z 

i 

) = 

(nt  -nt)(nt 
y z z y z X 

n t ) (n  t -n  t ) 
X z X y y X 

< 

1 

n 

V — > 

n 

— X 

n 

y 

n 

zLj 

k 

C(f>  cA 

L "^g  g 

C(t)  sA 

8 8 

S()) 

g 

k 

kr--2 

(B.8) 

The  ellipsoid  coordinates  and  their  derivatives  needed  in  eqns 

(B.2),  (B.3)  and  (B.4)  are  determined  in  terms  of  satellite  motion  by 

equations  (A. 20),  (A. 21)  and  (A. 22) . Thus  the  vectors  R , ^R  . and 

8 S g 

characterizing  the  sub-point's  position,  velocity  and  acceleration  can 
now  be  calculated  from  equations  (B.l). 

B.3  Radius  of  Curvature  of  the  Sub-Point  Path's  Pro-jection  in  the 
Osculating  Tangent  Plane. 

With  reference  to  Figure  6.3,  define 
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p = p(t)  = instantaneous  radius  of  curvature  of  the  groundtrack' s 
instantaneous  projection  into  the  ellipsoid's  tangent 
plane 


1 

— = c 
P - 


®d^R 



ds^ 


(B.9) 


where 


= tangent  plane  component  of  groundtrack  (sub-point)  acceleration 
with  respect  to  earth-fixed  axes 

c = unit  vector  normal  to  groundtrack,  lies  in  the  tangent  plane 

= ri  X _t  (B . 10) 

_t  = unit  vector  tangent  to  the  groundtrack 


V 


'R 


i + (y  + (2  - ^^)k] 


(B.ll) 


II  = unit  vector  normal  to  the  ellipsoid 

= (c';)cX)j^  + (c$sA)j^  + S(|i)k  (B.12) 

= speed  of  the  sub-point  with  respect  to  the  earth-fixed  axes 

(B.13) 


V = l^R 


= ^ = [ri  - ^ + (y  - H ) + (z  - H ; ]' 

-g'  dt  *■  a:  U 2 


s = arc  length  along  the  groundtrack 
t 

= / V(t)  dt 
o 


(B.14) 


“dR 

— C 
£ 

ds 


'dR  1 • 

. dt  ^ 1 ej^ 

dt  ds  V — g 


(B.15) 


®d^R 


= 


ds 


2 ds  ^ ,,2  ^ 


^R  + - 
V 


2 -g 


(B.16) 
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Substitution  of  Equations  (B,l)  and  (B,16)  and  then  equations  (B,ll) 
(B,12)  into  (B.IO)  and  (B.9)  ultimately  reduces  (B,9)  to  the  explicit 
formula 


= { (x  ~ H ) [Tb  - .7  jc(fis\  - (u  - H /)s!{>3 

+ (y  - H )lfx  - H ) - (z  - H JciPcXj 

y X z 

+ (z  - E )\_(y  - E Jc(Jc\  - (x  - E Jc(|)sX]} 

z y ^ 


(B.17) 


Equations  (B.13)  and  (B.17)  provide  the  important  equations  for  calculating 
the  sub-point  velocity  and  the  sub-point  path  radius  of  curvature. 
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APPENDIX  C 


4 CYCLE  RUNGE-KUTTA  ALGORITHM 

Given  a system  of  first  order  ordinary  differential  equations  of 


the  form 


— = f (t,x  ,x  , . . . , X ) , i=l,2  . , . ,n 
1 i z n 


(C.l) 


and  specified  initial  conditions 


{x  (t  ) , . . . X (t  ) } , (C. 

i o no 

the  following  4 cycle  Runge-Kutta  algorithm  permits  recursive,  step- 
by-step  integration  of  equations  (C.l)  to  determine  the  at  sequence 
of  times  t^.t^,  • • • 

1 k+1  ik  o il  2i  3i  4i 


i = 1,2, ... ,n 


where 


^l^i  " ^ (C.4a) 

A^x.  = At[f . (tj^  + 1^,  x^(t|^)  + , . . . , + -Y^)li=l,2,  . . ,n 

(C.4b) 

A3X.  = At[f.(t^^+^  , x^Ct^^)  , x^(t^)  +-^  )],  i=l,2...,n 


(C.4c) 


A^x.  = At[f.(t^^3,  Xj(tj^)  + A3X^,  ..,  x^(t^)  + A3X^)],  i=l,2,...n 


Sfl  = ‘^k  ^ 


(c . 4d) 
(C.4d) 


The  step  size  At  must  be  determined  empirically  to  maintain  the  desired 
number  of  significant  figures. 


54 


The  above  Runge-Kutta  algorithm  is  implemented  in  subroutine 


RIRJGE.  For  each  set  of  differential  equations  of  the  form  (C,l),  the 
functions  on  the  right  hand  side  must  be  programmed  in  subroutine 
DERIV. 

The  subroutines  are  given  in  Appendix  D,  set  up  for  the  integrations 
of  §4.1.1,  eqns.  (11). 
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APPENDIX  D 


IMPLEMENTATION  AND  DOCUMENTATION  OF  COMPUTER  PROGRAMS 
The  various  programs  developed  in  the  course  of  this  work  can  be 
divided  into  four  major  categories,  corresponding  to  the  four  major 
functions  of  this  software.  The  first  group  (D,l  through  D,12)  consists 
of  the  subroutines  necessary  to  compute  the  coefficients  for  the  Fourier 
series  fit  to  the  satellite  groundtrack  projection  and  the  f--function*. 
The  second  group  (viz.,  Forward  Transformation),  subroutines  D,13 


through  D.16,  generate  map  plane  x and  y coordinates  given  ({>  and  X, 

The  third  group  (viz.  Inverse  Transformation),  subroutine  D,18, 
generates  ({)  and  X given  map  plane  x and  y coordinate.  And  the  last 
group  (viz.  sensitivity  ananlysis)  , subroutine  D.17,  determines  length 
distortions  along  lines  of  constant  (Ji  and  A-  In  addition  to  the 
subroutines  cited,  the  final  three  groups  of  subroutines,  rely  upon 
several  of  the  first  group  to  perform  various  secondary  calculations, 

D.l  Subroutine  ROC  (T,  cons) 

This  subroutine  computes  the  radius  of  curvature  (eqn.  B.17 
Appendix  B)  and  the  first  and  second  derivatives  of  the  motion  along 
the  satellite  groundtrack  (eqns.  A. 16  and  A. 22,  Appendix  A).  Time  is 
input  es  T.  The  vector  C0NS(30)  contains  the  various  required  constants 
and  initial  conditions  defining  tlie  orbit  and  tlie  reference  ellipsoid. 
The  output  of  this  routine  is  passed  through  common  blocks  and  CONS, 

The  calculations  arc  based  on  the  eqviations  of  Appendix  A and  B and 
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TERM  = - 
P 

N(3)  = n 


DRDS(3)  = i 

A A 

CC(3)  = c = n X t 


DRDT(3)  = r:  - H 
DDRDT(3)  = r - H 
PHI  = 

g 

LAMDA  = X 

8 

DSDT  = |®R  I = V 
-g 

DDSDT  = V 
DPHI  = 

g 

DLAMDA  = X 

g 

XN(3)  = X,  Y,  Z 
DXN(3)  = X,  Y,  Z 

DDXl^(3)=  X,  Y,  Z 


_r  = satellite  position  vector 
H = height  above  the  surface  of  the  earth 


required  for  scan  vector  calculations 


Subroutine  ROC  has  external  references  to  ORBIT,  ROTATE,  CROSS,  EFRAME, 


PHIH,  DPHIDH,  and  VECPRD. 

D.2  Subroutine  ORBIT  (X,XF,TI ,TF,CONS) 

This  subroutine  uses  the  f,g,f,and  g solution  of  Appendix  A to 
compute  the  orbit  state  at  time  TF  for  all  species  of  elliptical  orbits 
(Including  circular).  X is  the  Initial  state  vector  at  time  TI  and 
XF  is  the  state  vector  at  time  TF.  CONS  is  the  vector  used  to  pass 
various  constants  and  the  output  of  subroutine  ROC,  The  position. 
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velocity  and  time  are  given  in  units  of  kra,  km/scc  and  sec  respectively. 
Subroutine  ORBIT  has  an  external  reference  to  NEIITON. 

D.3  Suhroutine  NEWTON  (CONS,  PHT , TI , TF) 

This  subroutine  uses  Newton's  method  to  iteratively  solve  Kepler's 
equation  for  change  in  eccentric  anomaly  E (eqn,  A. 8 Appendix  A), 
Subroutines  ORBIT  and  NEWTON  provide  the  means  to  compute  the 
satellite  position  at  arbitrary  given  times  as  a function  of  the  given 
initial  conditions. 

D.4  Subroutine  ROTATE  (M,TO,T,WE) 

This  subroutine  computes  the  earth's  rotation  direction  cosine 
matrix  M(3,3).  TO  is  the  initial  time.  T is  the  final  time.  And 
WE  is  the  earth's  rotation  rate.  6 is  passed  in  CONS. 

D.5  Subroutine  CROSS  (A,B,C) 

This  subroutine  computes  the  cross  product  of  two  vectors  (e.g. 

A X B = C)  . 

D.6  Subroutine  EFRAME  (A,B.N) 

This  subroutine  rotates  given  components  of  an  arbitrary  vector 
from  the  inertial  ref.  frame  to  the  instantaneous  earth-fixed,  equatorial 
frame.  A(3)  is  the  vector  to  be  operated  upon,  B(3,3)  is  the  rotation 
matrix,  and  N is  dimension  of  A.  The  output  overwrites  the  input  and 
is  returned  in  A(3). 

D.7  Subroutine  PHIH  (RXY , Z ,RN ,PHI ,H ,LAMDA , A ,B) 

This  subroutine  uses  a 2-dimensional  Newton's  method  to  find  phi 
(geodetic  lattitude)  and  H (the  height  of  the  satellite  above  the 
reference  ellipsoid  surface).  The  constants  and  corresponding  definition 
is : 
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Z = Z 

2 2 '' 

RXY  = (X  + y )'^ 

0 o ? !-- 
RN  = (x  + y + z ) ^ 

A = earth's  seroi-major  axis 

B = earth's  semi-minor  axis 

X,Y,Z  = Earth  fixed  components  of  the  satellite  position (x,y,z) 

See  equation  A, 22a  in  Appendix  A, 

D.8  Subroutine  DPHIDH (A ,B ■ H , HP ,HPP ,PHI ,DPHI ,DDPHI .X .RXY .R,MU, DRPT .DDRDT) 
This  subprogram  computes  the  first  and  second  derivatives  of  PHI 
(geodetic  lattitude)  and  H’(the  height  of  the  satellite  above  the 
ellipsoid  surface)  (see  Appendix  A eqn.  A, 22b  and  A, 22c),  X(3)  is  the 

vector  containing  the  earth  fixed  components  of  the  satellite  position. 
The  Fortran  definition  of  the  output  are: 

HP  = H 
HHP  = H 
DPHI  = Ip 
DDPHI  = (p 

D.9  Subroutine  VECPRD  (A.B,C.NU) 

This  subroutine  computes  the  inner  product  two  vectors.  If  NU 
/ 0 then  C = A*]5.  If  NU  = 0 then  C = A-A, 

D.IO  Subroutine  INTEG (Z , CONS , TO .TF) 

This  subroutine  integrates  to  determine  the  satellite  groundtrack 
map  plane  coordinates  as  well  as  the  coefficient  for  the  Fourier 
series  representation  of  these  coordinates. 

Z(40)  is  the  vector  containing  the  state  at  each  step  in  the 
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integration,  CONS(30)  is  the  vector  containing  initial  conditions  and 
constants.  TO  is  the  lower  limit  of  integration  and  TF  is  the  upper 
limit  of  integration.  The  components  of  the  Z vector  are  defined  as 
follows : 


7.^(t)  = 


X (t) 
g 


= / V(t) cos(f (t) )dT 


to 


t 

^ ^g(0)‘^  ^ V(T)sin(f  (T))dT 

Z (t)  = f(t)  = / dT 

to  p(t) 

t 

Z^(t)  = A^(t)  = -p  I v(t)cos(f  (T))dT 
to 

y t _ 

Z .^(t)  = A (t)  = -f,-  c [x  (t)  - X T]sin(^— )dT 
n+3  n g g ^ 


n=l,2-9,  10 
P 


Z^^(t)  = B^(t)  = - p / y^(T)d 


1 1 
2 P 


to 


"n+13 


(t)  = B^(t)  = -^  / y^(T)cos(^^)dT 


to 


= C (t)  = I / f(T)sin(^)dT 


n+23 


to 


In  this  program  one  preliminary  integration  is  performed  to  determine 

X and  y . Note  y„„=  B,  although, it  is  not  used  in  the  Fourier  series 
g og  i 

representation  of  y (t)  (so  that  v (t)  will  oscillate  about  the  x axis 
‘ g ■ g 

with  zero  mean)  . 

The  sole  output  of  these  subroutines  consists  of  the  Fourier 
fit  coefficients. 

Subroutine  INTEG  has  an  external  reference  to  RUNGE. 


60 


D.ll  Subroutine  RUNGE(TK,N, PELT, CONST. XK) 


This  subroutine  uses  a 4-cycle  Runge-Kutta  algorithm  to  approximate 
the  integration  of  the  state  XK  over  the  interval  [tK,TK+DELT] , as 
described  in  Appendix  C. 

TK  + BELT 

i.e.  }«(TK  + BELT)  = (TK)  +/  )aC(T)]dT, 

TK 

N is  the  dimension  of  the  current  state  vector,  Const(30)  is  the 
ubiquitious  vector  containing  the  necessary  constants! 

Subroutine  RUNGE  has  external  references  to  ROC  and  BERIV. 

B.12  Subroutine  BERIV(T.N,X. CONST, F) 

This  subroutine  computes  the  derivative  approximations  required 
in  subroutine  RUNGE.  T is  the  current  time.  N is  the  dimension  of 
the  state  vector  X.  CONST(30)  is  the  vector  of  necessary  constants. 

And  F(N)  contains,  as  output,  the  derivative  approximations. 

With  subroutine  BERIV,  we  complete  the  section  for  determining  the 
satellite  ground  m.ap  plane  coordinates. 

B.13  Subroutine  NS(CONS.PHI  ,LA.MBA,ALPHAO,XT.YT.TIMF,NU) 

This  subroutine  computes  the  x and  y map  plane  coordinates,  given 
ellipsoid  coordinates  PHI  and  LAMBA.  CONS  is  the  vector  containing 
necessary  constants.  ALPHAO  is  the  initial  equatorial  right  ascension 
displacement  angle  of  the  spacecraft  relative  to  the  inertial  X axis 
(usually  taken  along  the  intersection  of  the  initial  greenwich  meridian 
plane  and  the  equatorial  plane).  TIME  is  rime  t*  as  found  in  NS, 

If  NU^  100  seach  for  t*.  If  NU=100  bypass  the  t*  time  search,  XT  and 
YT  are  the  map  plane  coordinates. 
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Subroutine  NS  has  external  references  to  SETCON,  ROC,  ROTATE,  CROSS, 
EFRAME,  and  FNTS2. 

D.U  Subr outlne  SETCCN  (A. B, PIE, PERIOD. COFiS) 


This  subroutine  sets  parameter  values  not  passed  in  CONS,  A 
is  is  the  reference  ellipsoid  semi-major  axis,  and  B is  the  semi-minor 
axis  in  km.  PERIOD  is  the  satellite  orbital  period  in  seconds.  CONS  is  vector 
of  necessary  constants. 

D.15  Subroutine  FNTS2  (T , CONS , XCOEFFS . YCOEFFS , FCOEFFS.  F.  NX,  NY,  NZ) 

dx  dy 

^ ^ S dr 

This  subroutine  generates  values  for  x ,y  ,f,  -r-^,  and  by  the  Fourier 

g g dt  dt  dt 

series  representations  and  the  differential  eqns.  for  x ,y  and  f.  T 

§ § 

is  the  current  time.  CONS  is  the  vector  of  necessary  constants. 

XCOEFFS,  YCOEFFS  and  FCOEFFS  are  the  Fourier  series  coefficients, 

NX,  NY  and  NF  are  the  number  of  coefficients  for  x , v , and  f 

g ' g 

respectively.  As  output 
F(l)  = x„ 


F(2)  = y 


g 


F(3)  = yiti.c  ~ Franc  rotation  angle) 

F(4)  = f 
dx 

DX  = — ^ 
dt 


d V 

DY  = 

dt 


DF  = 


df 

dt 


D.  16 


Subroutine  FNTS2  has  external  references  to  SETCON  and  ANG , 
Subroutine  ANG(C,S,T.N) 

This  subroutine  computes  multiple  angle  sine  and  cosine  terras 


for 
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the  Fourier  series  for  x , y and  f.  C is  the  vector  containlnz  the 

g g 

cosine  (NT)  terms.  S is  the  vector  containing  the  sine  (NT)  terms. 

N is  the  number  of  terms  desired,  and  T is  the  time  which  has  been 
normalized  to  T = ir  (real  time) / (satellite  period). 

This  completes  the  section  for  determining  the  forward  transformation. 
D.17  Subroutine  ERR0RS(A,B,C0NS,FF,HP,K,U1,U2,  VI,  V2 , ANGl,  ANG2,  T) 


This  subroutine  computes  the  partial  derivatives 

3 b 
^x 
3 a' 


U1 

U2 


3 V 

VI  = 3b 
V2  = 

oX 

0 1 

and  the  length  distortion  factors  HP  and  K.  HP  is (r ) for  lines  of 

3s 

3 s ' 

contrast  A and  K is  (-' — ) for  lines  of  constant  b. 

os 

A = earth  semi-major  axis 
B = earth  semi-minor  axis 
FF  = flattening  factor 
ANGl  = lamda  of  current  point 
ANG2  = phi  of  current  point 

T = current  time  (as  found  in  NS  or  INVERSE) 

Subroutine  ERRORS  has  external  references  to  ROC,  ROTATE,  CROSS, 


VECPRD,  and  FNTS2. 

This  completes  the  section  on  sensitivity  analysis 
D.18  Subroutine  INVERSE  (TO,TI ,CONS , FF,X , Y , PHI ,LA>fDA) 

This  subroutine  computes  PHI  and  LAMDA  given  map  projection  plane 
X and  y coordinates.  TO  is  the  initial  time.  TI  is  the  current  time. 
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FF  is  the  flattening  factor. 


Subroutine  INVERSE  has  external  references  to  SETCON,  ROC,  ROTATE, 
CROSS,  EFRAME,  FNTS2,  NS,  and  ERRORS. 

This  completes  the  section  for  the  inverse  transformation. 

D.19  Subroutine  Scanner  (T,  R1 , R2 , R3 , XI,  X2,  X3,  Nl,  N2.  N3 , 

H.  DRl,  DR2,  DR3.  Cl,  02,  C3) 

This  program  computes  a linear  correction  to  time  when  the  scan 
speed  cannot  be  considered  infinite. 
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SPACE  OBLIQUE  MERCATOR 
PROGRAM  LISTING 

(Programmed  by  James  D.  Turner) 

The  subroutine  listings  follow  In  the  following  order: 
Subroutine 

SOM 

NS . . 

ERRORS  

SETUP 

INTEG  

ANG 

FNTST2 

SETCON  

INVERSE  

RUNGE  


ICS  . . 
DERIV  . 
ROC  . . 
ROTATE . 
PHTH.  . 
CROSS  . 
E FRA.ME 
ORBIT  . 
NEWTON 
DPHIDH 
VECPRD 


Page 


66 

69 

74 

79 

82 

34 

85 

86 
87 


91 


92 


93 

94 


97 


98 

99 

ino 

101 

105 

104 

105 


65 


1 pROr,«aM  s^*'  c i\put  .Output  . tapess  inp/T  . tapf  a=ou  t put. tape  io*TAPfi? 

1 , T ApF  m . TAPEl * . T APE3b  » 

c 

THIS  ^MOGRAy  acts  THE  OkIVER  FOR  THE  S.O.*^»  PAP  PROJECTION 
SUrtHOJTT'JFS,  THE  POOF  OF  OPERATION  IS  CONTRuLLLO  BY  THE  INPUT 

UARIATt  r flL>lPER. 

IF  Nj’AnFR  = 0 compute  THF  FuURTFR  COE  I F F I C I F N T S TO  THE  S.G.t* 

IF  N'J'^hFR  = I USE  thC  FOPWAftO  TRAfKFORHATTON  MODE 

TF  NJ'»«r«  = 2 USE  The  inverse  trAnsFOP**ATTOn  MOnE 

IF  Nj'*«fR  = 5 USE  the  sensitivity  ANALYSIS  MQOE 


15  C 

COMMON/OX YVAP/OY . OY 
CO*'*^ON/NSnATA/COROT  «DRPT  . C«  Hi  xx  .N 
COM*-ON/ROrnATA/A,R.XO.TO 

COMMON/XYrnFrc/XCOEFFs.YCnEFFs.FCflFFFs.Nx.NY.NF 
?0  COMmON/LE/R 

nf^rNSlDV  xCnEFFR(20)  , yCOEFFSc  ?0J,FC0EFFS(20).F('*) 
ni^FNSIDv  XOlFi .CONS( 3n  ».7(40>.GT(3).R(3i 
ni**FNSI0V  nOR'*T  ( * I iDROT  (3)»C(3).XXlf,>.N(3) 

HFAi  K 

25  real  NiUAmha 

CALL  SET JP( XO. TO.T.CONS ) 

CALL  S€T:of:<  A*8fPTEtP,COH^i 
CALL  MOCin..rrNS) 

no  70  isi . X 

30  70  cnNS(i  ♦ ?p»sTx(n  - h*n(t) 

FF=rONS( 7 ) 

ESO=?.*»*  - FF*FF 
AsCONSlB) 

RrCONSlR) 

35  WR I TF < 6. Sono J 

WRITF{6.SnnS)FF.A,B 

6000  FOR'*AT  { ///U9X  . • EARTH  RELATED  PARAMETERS*//) 

6005  FORMAT(m*x,*  Flattening  factoh=* .r^s, Si /44x. • semi-majoh=»»E22.5' 

)/44X«*  S-''T-«T\0R  AxIS=*«EI7.S//) 

40  RAPEGsPIE/IR*^. 

RrAnr5»4)MtiMprR 

write  (b.4)riuM”ER 

4 format  < I 9 ) 

IF  « flUMPE^  .F  Q.  0 ) GO  TO  loon 
45  REAni  10«  1 )M.^|Y 

REAn(l0.9)lXrCEFFS<I).I=l.NX) 

REAPi 12.1  )M.N7 

»FAn(12.9»(YrrEFFS(I).T=i.NT) 

REAP)  14.1  »M.  -jF 

50  RFAnii4,?)iFrrEFFR(r),i=i,NFi 

1 FORmATI?!*.) 

2 FORMAT  I 5* 1 R , t 0 ) 
jRTTF (6. SOI  0 ) 

6010  format ( ///4 . • rOEFF Ir IFNTS  fur  the  fOUPIER  FIT  XG.Y6  AND  F*//) 

55  .jH  r TF  ( 6 . SOI  S I 

6015  FORMATlXiy.*  X-rOEFFS . • . 90X . • Y-rOEFFS . • . 20X . • F * COEF FS • ♦ / / ) 
wR ITFl fa, SnpO)  f < XCOEFFSj  T ) , YCOFFFS  r I ) , FCOEFF^C I ) » . Ul .NX ) 


BtSl'AVAIUBlE  COPY 


or, 


bO 


65 


70 


75 


ao 


85 


90 


95 


100 


105 


110 


60?0  F0R“flT(?'4y.ria.7,t3Xiflfl.7.l3i.Clfl.7) 
jGC  TO  lOlO 
TF(*Ji;'1tiE:^.tO.’»'SO  TO  10?0 
TF  ( .f  0. 3 » GC  TO  1030 
1000  COMTIfjUE 
TF=T0  ♦ 3 

CALI  INTCr,  r 7 .CCN*^  • TO*  TF  ) 

Msuni 

write:  ( 10*  1 X 

wRTtF  ( 10.  ?M  xCOCFFSn  > * 1 = 1 .NX  ) 
writf ( i?.i  T 

wRlTril?,?)(YCOCFFSn).I  = l.NYj 
WRITf  ( 14. 1 )M.^F 

WRITf  ( 14.P  » (FCOEFFSn  > . 1 = 1 .NF  I 
WR  I TF  ( 6. 601  0 I 
WR ITF I 6. 6 01 5 » 

WRTTE(6,60?0»(f  xCOEFF5>(  I ) , YCOFFFSf  I ) ,FCOFFFS(  I » ) . Isl  .fjx  ) 

CALL  EXIT 
1010  CONTINUE 

PHl=.il55A6nn33El 
LA»'nA=.?57S7n*986El 
PHl=-.?a®A446'68F-2 
L A^nA=- , 1 ?957f 163? 

C test  toR^ARU  TRANSFORMATlOig 

WR  ITF { 6 . ) 

3333  FORwAT  ( Hi  ) 

WRITF(6.4^00) 

4500  FOR'IaT ( ///39y . • FORWA^O  transformation  TO  FIND  X.Y  GIVEN  Phi.La^DA 
1 •//) 

URlTF  16. 4ev05lFV-T  .lA’^OA 

4505  FOPMATf^Sx,*  PHl  =*.Fl?.6,/SSx,*  LAmOA  =*,FIP,6//J 
call  NS  ( :nriS  ,Rh  I.LAMpA.OfX.Y.T.Ol 
wRiTr(6.4sin)T 

4510  FORMAT(55X.#  T - S T aP  = * • F I ? . 6 // ) 
gRlTr(6«4SlPl 

4512  F0R«AT(49x,*  •'AP  projection  CooRO  T N A t^S  • / / > 

WRTTF (6.4615 ) X , Y 

4515  F0RMAT(44X.*  X =*,E16.7..  Y =*.E1A.7//J 
CALL  exit 
10?0  CONTINUE 

XTf5T=.10nqR66l76E5 
YTERTr- . 5R586A03’ 5E2 
YTFSTrO. TQl  36?62E  5 
XTr5T  = 0.37757167(46E5 
C INVEoSr  TRANSFORf^ATlON 

WR  I TF  ( 6. 47*»n  ) 

4790  format  ( 1 -^1  ) 

WRI TF ( b. 4600  » 

4800  F0RMAT(3Rx,«  T^JV^RSE  TRANSFORMATION  TO  FIND  PHI.  LAMOA  GIVEN  X.Y 
1 •// » 

wRITF(b.46ns»XTEST.YTFST 
4605  FnRMAT»4-4X  . • Xs».F18.7,*  y = *.F18,7//J 

CALI  INV-TmsF  ( TQ  , T ,CONS,  FF  ,XTE6T,YTEST,PHI,LA«0A) 

WRITr(6.-4aiO»FHl.LAMOA 

4610  format  ( /3RX . • ThF  (PhT.LAMOA)  uF  ThE  "^AP  PLANE  POINT  (X.Y) 
1«,//54X.*  PUT  = • . E 16. 7. /54X , • LAMDA  =*.ri6.7///) 
wR I TF ( b . 461 n ) T 


60  WMVABlt 


67 


115 


120 


125 


130 


135 


140 


145 


150 


155 


160 


165 


CAll  EXIT 
10.^0  rONTlNUE 

C <;rNSTTTVlTT  ANALYSIS 

URITr(6.5‘>nO» 

5500  POPMAT  M -11  ,47>  . • SEf^SlTlvITY  «NAlySI5  •//  ) 

WRITE  ( 6 • 5*.n5  ) 

5505  F0PMfiTl3SY.*  length  0l5T0RTl0.g5  FnR  PtS.  S Y T R I C ALL  Y PLACED  ON./ 
1 roth  stjcg  oe  the  satellite  ground  track  for  the  •/. 

2S5X,*  DI  SPt  ArEyfNT  INCRE»*FNT  MfLTA  = 55.66  km,.//) 

DEI P=fl. 

NDtLP=DE.P  *1.001 
OTrP/DEL^ 

Ts-OT 

no  300  JJrl .NEELP 

TrT+DT 

ETArT 

CALL  ROCTT.CONSl 
PHTsCONS ( 1 4 > 

LAPnA=CO'iS(  15  » 

PP=PHl/RAoFG 
AL=LA50A/pa0FG 
WRITE (6. 375 )PF  » AL  * T 

375  F0R'»aT(////35»,.  PHI.LaMHA  OF  THE  GRO  JND  TR  acK  . • 2F  1 5 . 5// • 

1 39X.*  TTwr  along  THE  sATFLLlTf  GROUND  TRACK  = •«FiO,5//> 

C GROUND  TRACK 

no  3,3  1 = 1,3 

-33  6T  ( T »sXX  I T ) - H*\  ( I > 

G-(PIE/1 50.  ) .A*?, 
aCLTS.sS/4  , 

no  P50  1 1 = 1 .q 

nSN  = - S I I T - 1 ) •DELTS 
XN  = n.SN*C  ( 11  ♦ GT  ( 1 > 

YN=nSN*CC?l  ♦ 6T(?> 

7N  = nSrj*:  ( 3 > ♦ GT  ( 3> 

LA»*nA  = AT  AN?  { Y^  « XN  » 

RXYsSORTi XM*X5  ♦ yN*YN| 

RNsSORT ( XM.XN  ♦ YN*YN  ♦ 2N*2N» 

CALL  PMl4<Rxy,2N.RN«PHT,H,LAMnA,A,B) 

CALL  NS  CoflS  .FHl  »lA'^0*  t 0 . . X « T , T . 0 1 

ANGi =LA5Da 

ANGPsPHI 

CALL  ERRORS < A .B»C0NS»FF . HP . K « 1 1 1 . U2 . V 1 . V2. ANG1  . ANG2.T 1 
IFf Tl.EO.SlGn  TO  5900 
UR  1 TF  ( 6*  495 )hP . K 
4?5  FORvAT( 37y , 9F90.6/I 
GO  TO  6050 
5900  CONTINUE 

wRlTr(6*4951  )»-P.K 

4251  FOP**AT  ( 37y  , 9F20.6.  * (SATELLITE  GROUND  TRACK!*  / 1 

6050  CONTINUE 
250  CONTINUE 
3n0  CONTINUE 
STOP 
END 


63 


\ 


«;|J»  rtnJT  T T ‘I*;  f C'"  .1  A n)A  • aI  PMAT  . X T , yT  , T I-wE  . nu  j 


c 


c*..* 

c 

T.-t;,  Si-rOWAM  '•-)»'PUTr«;  t Afju  Y f*AP  PLA-JF  CnOWOINATFS  GIVtN 

5 

c 

[ r tt 

1 A A . 

c 

TF  u.)  r 

'.«r  As  FOl'NJ  IN  SUF<.  IrjVtRSE. 

c 

f Pt.T 

10 

c 

9-^1 

r OroOFTIc  LaTTIT.iOC 

c 

LA-r 

A = L''  i&l  T“f  r 

c 

TT'  r 

= CuP^ff-T  tI'E  In  thf  program 

c 

«.■  T»  .-T 

15 

c 

Y T 

= Y.-a?  pj  ANF  COuhOINATF 

c 

c*»** 

Y T 

= Y-  -ap  p|  ariF  co.i«DP  atf 

c 

cn^‘T/.j/‘‘ry’  Ao/Dt»nY 


20 


25 


rnv«0  i/t.  ,r,/Ai  fh 

^n»'''Ol♦/''«;'“AT^/l:nQ"^T  .L'RrT«r»M«  x .n 
^n^•^Ol;/yY^r■F^«/YCoEPF<;,  YroFrP<.FcnrFF<;,fjy,^jy.NP 
rn*'**ot''/T  irPTA!  /x%  tDj<r.,rnY'' 
ro^ «-r;M/«r A'  /Tf'.'^FnT.^.rwni.OT.r 

I'f  ^ » .''f(  3)  . yf.j  tOXN*  A)  .r  n*i(  3 > .W  ( ^ > .Dfc.(  3 » t T<  3 ) .DT<  3) 
nT*-rt,«;lO';  dqt  ( \ , CC  ( 3 » , ( < » , OuC  T < 5 ) , F i 4 » , F e 4 > , ROT  ( 3 , 3 ) , R I ( 3 » 

ro*' -c  .l/■^3^nAT■''/a,.^1  xu.TO 

roA’'*rr./L*:/'^ 

Ffpsi-j  j r ( 3 » i TTROn  3 ) • n«rT  ( ^ > 

PT'^rfsI'i'j  rO*)«  ( l » tS  ( 3 » 3»  , X(»s)  . xr  (A  j ,rR<  31  . TPt  3) 

r»u),xri)EFF<;(;>n>.YCij^FFc(2ni,FcorFrSf20» 

RFAL  A 


C 


35 


40 


45 


50 


55 


TPTrT  1 

FAll  SF  T:ni.  ( A fj.  P !E  • F . f OfiS  » 

C=CO^.<F'-»t  ) 

CL=rC/';lLft**^A) 

SL  = <slj{-*.'*rA» 

•^sA  • A/S?^T  f A*A  tC  tr  «•  P*P*S*S» 

c 

C CC-rPJTF  TmF  'TA-iTH  FTxFO  CO«POrjrM«;  OF  THE  POINT  TO  RE  MAPPED 

c 

H H 1 I • : •rt 

RH  ) =M*:  .<;l 

c 

C r,r  T T-r  IiTtIal  GLFcs  ori  T i .-.E 

no  1 = 1.* 

35  R U » =H  I ( I ) 

Psrofps  i *'  > 

Tl=RJl»»>M>  ♦ ♦ K{3I*R(3) 

T3ri>  ( 1 ) . * ■ R < * » • ‘^  ( 2 I •TONS  ( » * R ( 3 I •CONS  < 25  ) 

T 3 = rf).|b  < 9 « t •rr ( ?3  » * CO'JS  { 2*.  » •Cor.S  < 94  ) * CON^  f 25  » •CONS  1 25  » 
ALPMAr«*:'?5  I T ?/SOR  T 1 T 1 • T 5 ) ) 

C 

TFfi  A-«'LA.rF.PTE/9.,  , L A^PA  . i f.3.*Prr/2.»ALPHas2.«PIF  - ALPHA 


69 


TF  r n ^ r .n  '"C  to  27 

T Tz  t'nA  *V'/  t P . •OTF  > 

fefl  >7  ro^  T I ,U' 

C 

TfOti-  T = r 
rof  «;  ( 1 ft  » zr . 

TOL'IZ  1 .f  -^o 

feb  C 

3 cMi  ‘^or ( T T . r'*K‘;  1 

c 

c 

rAii  KOTf.TMBrT.TO.ri.t.E:) 

70  C 

c 

no  1 z 1 , » 

so  ')r>  ( T I =^'  ( I ) - ( T ( I 1 - Mtr'  ( I > ) 

7‘>  no  i*'.!  Tzi  , ^ 

4S1  nw  « T » =t  M » 
vzrn»  SM 
nVPTzCO*  N f 1 7 » 

c 

80  C vr  = 3n«:rTTcN  PTHttal  s^^acf 

C rvf.  r w’FiOrTTY  I'J  ir.FRTTAL  ^^»ACF 

C nnyf.  z f.rrrj  that lOt.  If.  rffE:RTT«L  «;pacf 

c 

c n rs  T <r  vrcTOR  irj  thf  scA(.i  dipfctton 

fl*)  c 

r.Ck^K 

c 

c 

CALL  CHIsiCfifj.OCvj.OU) 

90  C 

c 

C riL  IS  TmF  CEIRtVATivr  OF  THf  SCAN  VECTOR  IN  THE  E-FRAmE  *JITH  CC'PONENTS 

C Tf  rf'SlTAL  SPACE 

C 

9S  nuM»zC.HT\  ♦ wC*u<?» 

nu(5)={jjf?)  • vr.ua ) 

nur  ^)=ujc 
C 
c 

100  c pctat-  tmf  co*'POfjEriTS  or  u..iu  to  tme  e-fhahf. 

c 

CALI  FF^ft^f  III. ROT. 3» 

CALI  EF®A»t  ml.  «RrT.3» 

C 

105  C T IS  TMf  iiMT  vector  along  IRACK  rITM  COMPONENTS  IN  THE  E-FRAme 

no  1 * I z 1 . \ 

TIM  sOFOT  I T 1 /V 

15  OT  I T j =L"l^l  T I M /■/  - UVHT  •ORPT  ( M / ( V*V  ) 

110  C r IS  T..F  Ilf  IT  VECTOR  In  TMf  CROSS  TRACK  OTRFCTION  WITH  COMPONENTS  IN 

C ThF  r , 

C 

CM  I CH0S<;i  N.T  ,r  \ 

OPi'TzCb  ISt  9k» 


70 


120 


125 


ISO 


135 


140 


145 


150 


155 


160 


165 


170 


OL/'"i.  Ot  S ( ?'»  > 

T'-.T  rC05  ( :or  «; » 1 H » I 
«;f.T  sSii.i  tf  4 > t 

rL^T-co** « ro  J' < 15»  > 

<;i  ( Ts«;i I rn.’j*;  i 15  » > 

c 

c nf  I**  T'lr  rcorvATivr  op  thf  nop^'al  v/f.ctor  in  thf  f-fra»'e. 

c 

nrjn  » '.f  *ci r ♦''PHI  - crir*SLt,i*Oi  A'ida 

p^  O > .<;i  r T •^PMI  ♦ CGT •CL(,  r •Ot 

HN  ( »,  ) iC  '7  P «T 
C 

TALI  < Il'u  T . *■  ) 

PAL  I CRr>>«-  ( T . F t 
C 

C nr  I?  tuF  '*FRI\/ATlVr  op  thf  C V/FCtOR  TN  the  E-FRAMF, 

no  I = 1 . 

7 9 nc  M )zl  ( M > F ( T ) 

c 

n = - «'Jh  OT  fl  ) *T  ( t » ♦ nRrT{?)*T|P)  ♦ DrOT  ( 3 J •Tt  3)  ) 

TP=nP ( 1 > *nT ( 1 1 ♦ nR(2)*nT(?)  ♦ OR(3‘«nT(3J 
T3  = np  ( 1 ) ♦ T f 1 j ♦ rR(2i*T(?j  ♦ ilK(3t*T(3> 

74  = - ( f,HnT  ( 1 ) .r  ( n ♦ L’RrT{P»»C»P)  ♦ DROT  ( S J *0  ( 3»  » 

T^rnK ( 1 ) .nr ( 7 I ♦ nR(2)*nC(2l  ♦ nR(SJ*nC(3> 

TA  = nK  ( 1 ) •r  f 1 ) ♦ rm?)*r(?>  ♦ Iir(3)*C(3I 

c 

c OR  Tt.r  FRnjFCTiOf  op  thf  oiRPLAcEMErjT  vector  onto  the  t.c  Plane. 

c 

no  3^4  7=7.3 

3.39  nR(7>=T'4HT>  ♦ TF,*C«r> 

c 

C nPLl^nr  TS  the  OERIVAtIVE  Op  thf  OT^iPUACET'TNT  vector  in  the  T»C  plane 

C 

no  49  1=7.3 

49  DEL  RET ( I ) = ( T1 4 73 7 ♦ T I 1 » ♦ T3*OT<I)  4 (T44T5)*C(T)  ♦ T6*DC(I» 

Tl=ni  (1»»T(7>  4 niJ(?)*T(?1  4 nu(37*T(3) 

T p = i I « 1 7 • Tt  ( 7 J 4 L(27»0T(?1  4 ij«37*r.T(3) 

T3  = '|  ( 1 7 ♦ 7 ( 1 j 4 'Jf37*T(3)  «.  U(3I*T(3> 

T4  = nii(  1 7 ip  n 1 4 nil(274C(?7  4 Imi(37*C{37 

T*>  = M f 1 > • Tr  ( 1 ) 4 i’(27«nf(P7  4 iii37»nC(3J 

T6  = 'I  n 7 ♦ *:  I 7 » ♦ 1'(274C(?7  4 U(3)«Cl.37 

C 

rFR'7=T34  ( T1  4T?  ) 4 T6#{T4*T57 
^77f‘'«iN'5Sn^  T ( T 3*  T ^ 4 Tb*TF,j 
Afl  = T3/.s  jvco 
R7J  = TF/SJ'^R0 
C 

C . !«:  r og  TmE  SCAN  VECTOR  as  PROJECTEP  ONTO  THE  T.C  PLANE 

C 

no  us  1=1.1 

4 5 - ( f 7 =Aa  •rut  * Rn«C  < I ) 
r 

r.As  ( n ♦ T 5 » / sM^so  - T3•Tr«^♦/^  SH.-iso^suMSOtsiiMiira  7 

ripr  7T*44T3t/S"*'^T  - Tb*7ERvi/(  Si  mSO  • RU'^SO  • S(  |MSO  7 


71 


175 


180 


185 


190 


195 


200 


205 


210 


215 


220 


225 


C ^r.tU  !«:  7‘'C  -:r.<IVATIVf  OF  !•••?;  SCAN  VLCTOW  As  PROJFCTCD  ONTO  TmE  T»C 

C nj  fti.r 

C 

no  *?!  1=1.^ 

‘'I  riv''T  ( I > = ■'A  • T ( T ) « AA*riT(ii  ♦ ♦ BR^ncrii 

c 

c 

CALI. 

cAi  I n-'C  T .r ) 

fAl.L  I 

c 

no  1 T =1  . ^ 

Ho  r ( T )=E  I T 1 4 r ( I ) 

Tni-;=(;H  ( 1 1 * TP  ( 1 ) + np(?)*Tp^2)  ♦ np«^)*TP<3) 

TFi-'vt  =00  ( 1 ) .r  ( 1 ) 4-  nh(?)*F(2)  * Dpi  3)  •EMI 

7FR'^?  = U'‘,.Pr‘T  n ) •TPU)  ♦ OFLP07  I2>*7P{P>  DC  I HOT  I 3 > • T p ( 3 » 

C 

PT  = Tf  R>’" 

pFf  T=Tf  ^^7  * 7£:R‘^2 

TF  f 'ji  .tO.inn^rc  TO  555 

c 

Ti=Ti  - -T/nrrr 

c 

355  rof  TT'iUP 

I ip  = nv  ( 1 ) ♦ r { 1 I ♦ ''.j<2i*C(p>  ♦ ri«(.3i*C(5) 

TF  I Ml,  L".  1 mil  GO  TO  20 
IF ( Af  5( tt-tOI  r ) ,lE. I .E-051  GO  10  ?0 
TOL  n=T I 
C 
C 

TCniir.  r = T :niiOT  + i 
TF(  irOU'lT.rQ.**  I GO  TO  37 
IF  ( TCOO'iT  .F  0 . ® J CONS(1f)=1 
TF  ( rO':i>  RETiiPrj 

GO  TO  b 

^7  •’RrTF<i^*innni 
lOl  O FOP-A  T ( 1 HI  I 

WRT  TF  { fa  I 1 lien 

lOl.l  FOF'-AT«*  FATIFT  to  CntVFHGt  In  HS  TT'*E  SFARCH.//) 

C 

GTOP 

20  cnr.TrjU” 

IF  { 'III.  £ 3.  ionise  T ) 4 11 

CAM  sCA'.r,rRrTi,Ri(l),Fli?)t«T«3).itll)iX(2I.ir(.3)«N<l)«N(2).N<3)« 

1 h,r}''(i>»rv(pi»Cj'"i5>ir{ii,C(?>«C{3)i 

51 1 rnr'T  inor 

TT-FrTl 

CALI.  fr.Tc;9,Ti,cn'G«XCOfFFS,yCi)7FFc.FcnCFr«;«F,NX*NY,r>;Fl 

C 

TTiGf.rtT  ( ? ( 1 > •«  ( 1 1 * H(2I*m<2»  * k(3)*M{3n 

nn  ^*'0  1 = 1 . \ 

3^0  nn I T I < I I - F •% { T I 

c 


72 


230 


235 


240 


245 


BHsr  available  copy 


C !«;  THr  r^AriUS  OF  riiRVAjuPF  IN  THE  CROSS  TRACK  niRECTION 

C 

RCTrf.H  t I ) ( 1 1 ♦ nR(2)*M(?)  + i)P  ( ( 3 > •A*A/9/P 

C 

C AS  OPTnTfjALLY  OEFirjEO  THE  kCT  VECTOR  IS  DOWN  THE  NEC.  N-AXIs. 

C 

RCTrAtjST  RFT  | 

ALT  HAC  = A r f.N?  (I  P.  RCT  ) 

TF<APS(A^PHAr).3T.(20.*PlF/l6n.I )rONS( 16)=1 
TF(rrjriS(iA>.FC.i  )REruRf. 

TFP«=TA*'(  ill  PH^C/?.  ♦ PIE/4.) 

TEP^?=W ( n * T t 1 T ♦ W(2)*T<P)  + w(3>*T(3) 

TMrTAS  = -AST^ttTE9vp) 

nsRCToALOOTTFR'^l/COSTTHETAS) 

XTrFtll  ♦ n*r''S('^l3>  ♦ THETAS) 

YTrFfP)  ♦ n*STNff'(^»  * TmFTAS) 

C 

RETURN 

FNC 


73 


1 


SUBROUT  I CPRQRS  t 'B*  CONS  .FF  .Hp.  K «U1  .U2  . VI  t y/2  . A%Gl  • AnG2  « TI^E  ) 


5 


10 


15 


20 


25 


30 


35 


HO 


H5 


50 


55 


C 

C«a*«*»«*««««***.*«««*»««*«***«««»**«****««»««**«*******»**«*«*****«**«**«**»« 

C THIS  RROGRAK  CO^BUTrS  THE  paHTTAlS  OEWIVATIVES  KEQUIHFD  TO  DETERMINE 

C LFNGTH  OTSTcRTIONS  TN  the  map  plant.  ThESF  pARTIALS  ARE 

C ALSO  Jsrn  9Y  SUBR.  tNVFRSC  to  fTIO  Phi  and  lAMGA  given  X AND  T. 

c 

C INPUT 

C A = SE'^I-MAJOR  Axis  OF  THE  EARTH 

C 3 = «E'*I-MIN0F  axis  of  THE  EARTH 

c PF  = Flattening  factOh 

C A\Gl  = LAMDA  OF  The  mapping  PT.  OFF  THE  GROUND  TRACk 

c ArG?  = PH  OF  THE  “APpInG  PT.  OFF  THE  GROUND  TRACK 

c TTMf  = THE  current  time  As  foiino  in  subroutine  ns. 

C OUTPJT 

C OARTTAlS 

C J1  = rx/C(PHl  ) 

C J9  = CX/0(LAHDA> 

c VI  = ry/n(PHi > 

C V?  = rY/C(LAMDA> 

C ITSTORTICNS  factors 

C HP  s LENGTH  distortions  ALONG  LINES  OF  CONST  LAMOA 

C •<  = LENGTH  distortions  AlONG  LINES  OF  CONST  PHI 

C«*««**«*«**«««***«*«««««*«««t*««««*«««**«**»»**»«««»*«*««*«***»«**»»*««»**««« 

C 

COMMON/NSnATA/CORDT.DRnT,CC.H,*,N 
COMMON/OF/DFT 
C0Mm0N/D*YMAP/U<CT .DYDT 
COMmoN/AVG/AlFhAC 
COMmon/OErviS/TEST (20) 

CO»»MON/XyrnrF«/XC0CFFs,YCnEFF|!;,FcnrFFS'NX.Nr,NF 
COMMON /sc  AN/ W.  CFC  T tWP.rii.PnT  tOTOT  . TP 
dimension  nTOTtJi .OTPOPHI ( 3 ) .niPDl  AM( X ) , TP  < 5 » 

niMFNSION  RI(3),  X 16 ) ,DR0t ( 3 » .DOROT J 3 » ,T( 3) .SS( 3» .CONS( 1 ) 

DIMENSION  nNOPhl ( 3 » .DNCLAMj  3).HUMmy(3».E(3).F(h>.G(3».N(3» 
dimension  wP  » * I . EwPOT  < 3)  . OwPOPriI  ( 3 ) « OwPDLAM  ( 3 1 , CC  ( 3 > 

DIMENSION  xCofFrs(20l.YCOFFF3<20).FCOrFFS«20»  .R(3».DR(3) 

DIMENSION  u < X ) .DCDPHI (x),nCOLAM(3).DRRMl<3l.nKL*M(3».ROT<3.3) 

REAl  M.N.l  AMOA 
REAL  K 
TIrTIME 

Call  ROctTi.rCNS) 

PHI=C0NS< 1 4 ) 

LAmOksCONSj IS) 

VsCCNS( l«) 

DVOTsCONS) I 7) 
nPMl=CONS<  ?6) 

DLAMOA=C3NS(?7) 

C r-FRA-r  COMPONENTS  cF  S.G.T. 

RXsX  ( I ) - H»N< 1 ) 

RYrXI?)  - f4*»i(2) 

R2=X  f 3 • • H*N»  3 ) 

Pir=2. •ASTNJ 1 . ) 

CL  = COS<A'jm  j 
SLsSINIAVr.l  ) 

CsCOS I ANG? ) 

S=SIN<ANS2» 


74 


60 


65 


70 


75 


60 


65 


90 


95 


100 


105 


110 


m=a*a/S9^t « a* ft ‘C^c 

C rcf'PCNENT';  or  the  pt.  to  be  Mftppro 

R ( 1 )z‘^*C*CL 

R(?»=»^»C*«;l 
R ( S » =8«B*y*B/ft /ft 

no  uB  1=1 . X 

•*9  RI  ( T )=K  ( I ) 

TO=n, 

C WE  IS  THf  rotation  RftTE  OP  THE  FftRTH 

CALL  rotate (RCT«T0. TI .wE ) 

C r-FRA'*r  C0'»PCNENTS  of  THE  PT.  TO  RE  i“APPEn 

no  X 1=1 . X 

3 DR  ( T )=«  ( t ) - f X { I » - H*'J(  r ) ) 

TANR?=2 . •AMG9 
C 

OMrPMl  = (A*A  - B«o  ) •>1*"***^*SIN(  TUNG2) /<  2.  •a*A»A*A  ) 

C 

no  45  1=1 . X 

45  nUP*ftY( I ) snPOTT I ) /V 

c 

C OnyPnTE  TmE  COaPONENTs  OF  THE  C VECTOR 

C 

CALL  CR39S  f N«DU“!“Y,CC  » 

C 

C T(I>«I=1«X  . . .OR t PHI .LAMDA I /OPWT . nRsRrPHT.LAWDAI  - R(I) 

C Rin.Irl.X  , . .DR  » PHT  ,Lft?*DA  « /OLAPOA  , DR  = R ( PM  I . L AmO  A 1 - R(T) 

T( 1 IsDMO^ht *C*CL  - M*S»CL 
T(2»=Of*nRMT*r*SL  - P*S«SL 

T 1 3 ) = (fa*B/ ( a*fl n ♦ (Omdphi •s  + h*C) 

c 

ss(  1 »=-''*r*SL 
SSI 9»=H*C*CL 
SSI 3»=0. 

c 

C F = DB.ORnT  , WHERE  (.)=  INNER  PRODUCT 

C 

C TP=DRDT/V  Is  the  U^4TT  VECTOH  IN  ThF  INSTANTANEOUS  VELOCITY  DIRECTION 

C 

C W=N  CROSS  UP  .WHERE  wP  IS  THE  SCAN  VECTOR  IN  THE  T.C  PLANE 

C 

call  VECPPDl TP, T.DRPDOTT , 1 ) 

CALI  VECPRniTP.SS.DRLDOTT.ll 
CALL  vEC»Hn<CC , T.ORPUOTC . I ) 

CALI  vLC^pnirr .ss .ORlOotc , i » 

CALL  vCCPRP  ( TP  . W . TDOTi.  . 1 > 

CAl L VECPROICC .W.COOTW. 1 ) 

C 

nFnpHi=DRpnoTT*TnoTw  * orpootc^cootw 

nFDI  AHSD^I nOTT.TOOTW  ♦ ORL DC Tc •COOT « 

C nF/D|PHT »=rFDPHl  » nF/niLAMuAlsDFnLAM 

c 

nTDPHi =- TFOPm! /OFDT 
DTOl  AMs-TFOLA^'/nFOT 


DO  ?7  1 = 1 , X 

OTPopHl 1 1 j snrr  T I T » •Otophi 


75 


lib 


120 


125 


130 


135 


l40 


1^*5 


150 


155 


160 


165 


27  niPOLdWt  T ] rOTT  r ( T > •OTPL  A*- 
C 

C RFLOW  PH!.LA''Pa  COOttCSPONU  TO  THF  OROUrjO  TRACK  COORDINATCS 

C 

CsCOSiPHi ) 

•irSiNJPHI  > 

CL=ros<LAMnA» 

«;L=<;iN(LAMnA) 

RC  = RX*C*Ci.  ♦ RY*C*SL  ♦ ( A • A/H/rt  ) •P?  ts 
ORCnPHi  = -KX  •«;*CL  - RY*«;*SI  ♦ « a* A/P/R  j »R2*C 
ORCnLAMr-PY  •r*SL  ♦ RY*r*rL 

c 

TFRMrDHCTPM  I *CPMT  ♦ DRCPLA^'*Oi  AMDA 

TrPwPsOROTn  » •C*'‘L  ♦ nRDT(2)*r*SL  ♦ A*  A*PROT  J 3 J ‘S/B/B 
C 

DRCnPMl s I TCRM  ♦ tERM2) •OTDPMl 
ORCDL AMs  I TFRv  ♦ T E RM2 ) • DTOL AM 

C 

ANGsPIE/4.  *■  AlppAC/2. 

TlsTAN  < A'JR  ) 

T?=1 ./COSf ANP)/CCS<ANG) 

TEP'^siJP  M » • rP  ( 1 ) ♦ WP(?)«TP(2l  ♦ WP(3»*TP<3» 

thftas=*a5;tn(  ter*'  ) 

0=Rr*AL03 ( Tl »/C0«(ThETaS> 

C 

C here  WF  RiiIlD  the  PAPrrALS  uF  the  c vector. 

c 

no  05  1 = 1 . X 

nRPHHl)=T(Ii  - PRDT f n •OTOPHt 
05  nPLAM  ( I » s^s";  ( T ) - OROT  ( I J •OTCLAM 

c 

TAUl =-S*Ct  •npH I - C«SL*0LAMDA 
TAU2=-S*Sl ♦np^ I ♦ C*CL*OLAMDA 
TAU3=C«^»HI 
C 

ONDPHI  ()  IsTAltltDTDPHl 
nfjFjPHl  ( r*  ) = TAt|2  •otophi 
DNPPHT  ( 3 ) sT  Ai  I’  *01 JPHI 
C 

nNOLAM(  l)sTAtll*OTOLAM 
ONOl  AM(?)sTAil2*OTOi.AM 
HNOl  AM|  5>sTAll?*nT0LAM 

C 

nvnPMi=nvnT*nTDPHl 

nvnLAM=ovnT*nTCLAM 

no  5 1=1 . X 

5 nuMMY { I ) SORPT  « I » /V 

CALI  CR0R<;  f DMCPHT  f DUMt^Y  . E » 

DO  10  1=1.3 

10  OtiMMY  I I ) snnRnT  ( I » •OTDPhI /V 
CALL  CR05<;iN.ru''MY*F) 
no  15  1=1.3 

15  nuMMY 1 n snROT f n •ovoPHT/v/v 
CAI  I CHnRC(N,ru'“*'Y.G» 

C 
C 
C 


170 


RFLOrt  v>r  cfi-pl)TE  the  nFRiVATiyr  or  the  c vector  w.r.t.  phi 


175 


180 


185 


190 


195 


200 


205 


210 


215 


220 


225 


C 


C 

c 

c 


c 

c 


c 


c 


c 


c 


c 


c 


c 

c 

c 


00  ?n  1=1.^ 

?0  OCOPhI  ( T »=E  J T > + Fd)  “ 6(1) 

no  1=1.5 
?5  nu*''^Y  ( I ) sr-ROT  ( 1 1 /V 

CALI  CROSS  (DrJCLA‘'.DU^*^Y,E:) 

no  50  1=1.5 

50  nu***'Y ( I ) sonROT  ( n •dtdlai«/v 
CALL  CKOSSJ  W.ru«*'Y  .F  » 

DO  55  1=1.5 

.55  nur*MY  ( I ) =nHnT  ( n •ovdlam/v/v 
CALL  CROSS{'i.ru''''Y.G) 

RFLOW  wF  CPNPuTE  THt  OFKlVflrlvr  OF  TMF  C VECTOR  k* . R . T . LAIQA 
no  40  1=1.5 

40  ocni  A«n  )=F(  T 1 ♦ F(I)  - 6(1) 
no  93  1=1.5 

P2  owPnpHi  ( nrOwFDTf  n •dtdphi 

DO  43  1=1,  5 

42  nuPnLAM(  T IsOWPCT  ( I ) •DTOLA'^ 


TER'^1=0R  ( 1 ) *60  ( n ♦ 0R(2)*CC(3I  + nR(5)*CC{3) 

TEfiM3  = OR3MT ( 1 ) tCC ( I ) + DRPHi ( 9 j *66 ( 3 ) ♦ OrPHT ( 5 ) ‘CC ( 3 ) 

TFRm5  = UR ( 1 ) •nCopHl ( 1 ) ♦ D« ( 2 ) •DCnpHK 9 ) ♦ OR ( 3 > •DCnpHl ( 3 ) 

nAcnpHi  = iHr*(YERw3  tfrn‘5)  - term  «[jRC0PHr ) /(rc*mc  ♦ termi^teh*'!  ) 

TEP«P  = DR:.Aw  f 1 ) tCC  ( 1 ) ♦ DRLAf*  ( 3 ) *cr  { ? » ♦ DRLAm  ( 5 ) *CC  ( 3 ) 

TEP'15  = 0R  ( 1 1 •nrcLAM(  I ) * dr  ( 2 » *tlCDl  A^C  3 ) + DR  ( 3 ) •DCDLAm  ( 3 ) 

nACnLA«=(wr*{7ER'“3  ♦ TFPV5)  - TERvi  •nPCDLA'5  ) /■  ( RC  •RC  + TER^' 1 • TE  R^' 1) 

P=CnNS(6) 

CALL  rr^r^pf rr«coA'S.)fCorFF<;, YCot ffc.fcoeffs.f w-^)f♦<'Jy.^FJ 

CF=C0S(P(4)  ) 

.sFssirn  p ( u ) ) 

ODXnT=OV)T*CF  - v*SF*nFT 
nnYnT  = DV')T*SF  ♦ V*CF*nFT 

nnxnpHi=DnYnT*cTrpHi 
nnxoL  AM=onxDT*oTnLA« 

DnYnpHi  = onYnT*cT''PHi 

nDYnLAMsOnYDT^CTCLAM 

CAi  L vec9p(0(nupnFHr . rP.o^pnoTT.  i » 

CALL  vEC3Rn(nk,pnLa«,TP,OwLOOTT.l  ) 

CALL  VtCPPn  f 'v’P  . OTPOPMI  .WDOTDTp.l) 

CAii  vEc=>Rn(wP.nTPnLA'*,rtnoTDTt  .d 

SU^S(J  = SORY  e 1 . - TFH'1*TrRM) 

nr,npMi  = (TYrT*roYrpHi  - DYOT«Du*t)pwn  / ( oxni^nxDT  ♦ dydt^dyoti 
1 -(OwPDOTT  ♦ WCOTOTP ) /SUMSO 

nr.Ol  A'liC  DxnT*''OYrLA’^  • OYOT*OuxDLA^*)/(DxnT*OXDT  ♦ nYOT*CYOT) 


best  AVAllAulE  COPY 


77 


J 


250 


1 -(nWLODTT  ♦ UCOTOTL )/SU«SQ 


235 


240 


245 


250 


255 


260 


265 


270 


C 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 


CMsCOSJThfTA*;) 


nDDPni  = ( :3rfrnpn*aLon(T  1 j ♦ KC*T2*nACUPHi/Ti/?.  \ /ch 
1 - «c*flLDr, ( Ti  1 « nwrnoTT  ♦ wuointp) /SU'^5Q/<  ch*cm  > 


DDOl  AM=»  ’)prnL'"'*ALOG{T1  > ♦ RC«T2*nACPLAW/Tl/?.  > /CH 

1 - RC  •Al.Dr.  t T1  » •Sh*  ( DwUrOTT  ♦ ^iiOTHTL  ) /Su^'^Q/ ( CH*CH  ) 
C = Cn5(F ( ^ J ♦ ThETAS  ) 

S=STh(FJ5)  ♦ ThETaS) 

OXDPHl  = TXnT  I 

nxnLAM  = P»nT  .nTcL  A'* 
nYPPHI=^r^T•^TCP^-T 
nYnLflM=nYrT*nTDLa« 

ro'^f'iTr  thE  partial  oerivatives 


rx/RiPHn.  nx/P(i«MnA).  r)Y/n(PHi>,  dy/diamua) 


LET  J = X 
y = Y 

1 = pmT 

2 = LA''0A 


ui=nxDPHi  ♦ nrcPHi*c 
UPsPXTLA'^  * nrjLAM«c 
VlsnYDPHi  ♦ nrcopt*S 
VP=nYOLA'^  ♦ nTiiLAH.s 

ESO=?.*p-  - rc*pp 
s=srfj(A»':oj 
WsSORTd.  - r«Q*<;«S) 
F1=M1 *01  ♦ VI *V1 

r?=n?»u?  ♦ v?*v2 


0*S*06fiHHl 
0*S*0GnLAM 
0«C»0G|)PMI 
0*C*OGf)l  AP 


COMP  /TF  LFAGTh  DISTORTIOMS  ALONG  LINES  OF  CONST  LAmDA 


HP  = U*U*  W*<;rj9T  I E 1 ) / I A*(1.-ES(JI» 

COMPUTF  length  niSTORTTONS  ALONG  L^NES  OF  CONST  PHI, 


KsW*SfiKT<F?)/(A*COS<  ANG2 ) ) 

RETURN 

ENP 


BEST  AVAIUBLE  COPY 


VECPRD 


105 


65 


5 


10 


15 


20 


25 


30 


35 


40 


45 


50 


55 


SUPfiOUTl'IF  SFTLP(XOtTO.T.CONSl 
niMFUSl^M  XO (1 » .CONS( 1 ) 

RE&L  41U 

REAL  '^UnfiYS.Mr 

real 

PIFz?. * AS  TH ( 1 . ) 

C PFRIOI  IN  «ECCN0S. 

P= ( 103.2^7 ) *60 . 0 
-«1U  = 396(>01.? 

C 

RO=fP*S5RTM1l)/2./PIL»**(?./3.I 
VO  = SCRT{ 

RADnFG=l . 745  32'9252E-02 

C POSlTTOr 

XOf1  )=0. 

XO(?»=RO*«ITM(«.*RAOOEG\ 
XO(3)=RO»ro^«?i*RADOLGl 
C VELOCITY 

XO ( 4 ) =VO 
XO(S)=0. 

XO(ft)=0. 

C TIME 

TO=ri . 

CONS (21 1 =RO 

c 

c orbital  elti-ems 

C * 

c 

C A SE'^T-PAjoR  AXIS  OF  ThL  ORBIT 

C F FCCFNTRTCITY  of  THE  ORBIT 

C ANGl  inclination  OF  THE  OhIBT 

C ANG0«  iONF^TTUFE  of  THE  ASCFnDING  NODE 

C ANGR  AHGii^'ENt  of  PFHIFOCUS 

C T TTY'F  OF  perifocal  PAsSAGF 

C 

TO  = n . 

4lJ=X.9fl50l  ?FO' 

SQ«l'  = 6.31T4f\7E02 

C POSiTTOt; 

x = xn ( 1 ) 

Y = xn ( 2 1 
7 = X 0 ( 3 ) 

C VELOCITY 

nxsxo ( 4 ) 

0Y=X0(5) 

nz=xo(bi 

C 

RO=SORT(x.x  ♦ Y*Y  ♦ Z*7) 

VO  = S(JRT  { D**OX  ♦ CY*OY  ♦ 07*0i) 
nr<  = x*DX  ♦ Y*nY  ♦ 2*02 
ATrP.O/RO  • Vf»VO/MU 
C 

C SFMI^A.IOR  Axis 

C 

A=1  .0/AT 

PfRTOO=2.*PlF*(A»*l.5l /SORT ( Mu i 
N = R0MU/  A/SCiRT  t A ) 


79 


60 


65 


70 


75 


60 


85 


90 


95 


100 


105 


110 


T1 =1  . n - U0/^ 

T?  = nO/SO'^U/SO»T  ( A J 

c 

c e:ccE‘jt«tctty 

c 

E=90RT(T1*T1  ♦ T?*T2» 

EO  = aTAN2(T?.Tl  ) 

»*o=ro  - 

TT  = TO  - '<n/N 

nnorwu* ( 1 . n/RC  • i.o/A) 

TFIE  .LE.  1 .r-lOJGO  TO  iQ 
T=i , 0/u/r 

PI=T*(DDO*t  - O0*0X) 

P?=T*<D03*Y  - C0*07» 

P^rT^lOOD*?  - CO*OZ) 

PsA*( 1.0  - E*E  ) 

HO=PO  - = 

DHO=DO/R3 
T = 1 ./50'^J/9QRT  (D) 

Ql=T*(0^0*ir  - hO*OX) 

Q2  = T*(D*J0»Y  - HO*OV) 

Q3=T»(OwO*Z  - hO*OZ) 

C 

ao  TO  12 
10  CONTINUE 

c 

C P ANO  Q VECTORS  FOR  thE  CIRCULAR  ORuIT  SPECIAL  CASE 

C 

Pl=X/RO 

P2TY/RO 

P,=7/RO 

Ql=nx/VO 
02=nY/V0 
Q3=nz/V0 
12  CONTINUE 
C 

W1=P?*03  - P^*Q2 
W?sP3*0l  - Pi *03 
w3sP1*QP  - p?*ai 
C 

ANGI=ACDSlw3) 

AN60M  = ATAfJ?{Wl  *-W?) 

ANG.rf  = ATA\|?(P3.(32» 
write(6. 1 nn ) 

100  FORMATdHI  ) 

WRITE(6.tn^) 

105  F0P«AT(45X , •TAPUT  PARAMETERS  OEFININR  THE  CURRENT  RUN, •.///» 
WRITE  < 6. 11 ni 

110  F0PMAT(47X,*  FCSTTION* . 20X. • VELOCITY*//) 

WRITE(6«11*^)XC(1),XO(4),XO(2J.XO(*.»*XO(J).XO(6) 

115  F0PMaTl45X.F1®.5.14X.Fl5.S) 

WRITE  <6,  Ipn  ) 

i?o  for^ati/z/sty.*  orbital  FLE^CnTS*.//) 

wfi  rTE(6.  I ) A 

125  FORMAT ( 3AX .FI  * ,5.5X . • SEMT-mAjjR  AXIS  OF  THE  0RR1T*J 
WRITE (6. 1 xn  )r 

130  FORMAT! ,Fi * ,5, 5X . * ECCENTRICITY  OF  THE  ORRIT*) 


BBI  AVAIWSIE  COPY 


r 


115 


120 


125 


WRITf  (6.  ns  M^GI 

135  roM«aT(3flv,m.5.5x,*  inclination  of  the  orbit*i 

WRITE<6.mn»ANG0'' 

l^O  FORHAT ( 3Ay ,F1 * .5,5X. • lONGITUt)t  OF  THE  ASCENOING  NODE*) 
yRITr(6. luS)A^6W 

145  FOR»^AT  ( 3ex  . Fl  5 ,5  t 5X  . • ARGu»*ENt  OF  PFRIFOCUS*) 

WRTTF(6. isn  )TT 

ISO  FORMAT  ( 3Ay  .FI  ' .5«5X  , • mt  OF  PERIFOCAL  PASSAGE*) 

WRTTF (6. 1 SS)PFr  icn 

155  format ( 3AX . Fi «, 5 , 5X . • SATFLLITF  ORBITAL  PERIOO*) 

RETURN 

ENP 


■31 


1 


5 


lO 


15 


20 


25 


30 


35 


40 


45 


50 


55 


SUFROUTI'JF  IMrfifZ.COfjS.TO.TF  » 

c 

C THIS  3RPGRA4  T'iTEGPflTrs  TO  i)ETFRMI»4E  ThE  FOURIER  COEFFICIENTS  FOR 

C THE  satellite  GROUND  TRACK  COORDINATES  AND  F. 

C**«««***»*«******«««*******«**«***««*««*«««**««*«**«*«******««*********«* 

c 

CO^MON/COlif  1 /'* 

rOf^'^ON/XYrOEFC/XCOEFFs.  YCOEFF<;.FCOEFFS«NXiNY  .NF 

oi^rr.oION  xco‘'FF^(2o»  .yCofffs(po»  .FcorFrs<2oi 

ni^ENSlON  7M»,C0NS(1> 
ni'^ENSlD'J  F f4) 

ILASTrO 

PIE  5?  . • ASTNM  , ) 
nELT  = CONJS<  1 0 1 
TsTD 
C 

CALL  ROC(n..CCNSl 
C 

TLASTsO 
T = TO 

CALI  1CS(7.C0^S•4  > 

DO  K 11  = 1. ^on 

IFJ{T  ♦ ori n .GE.  TF)  tLAST=l 
IF( TLAST.FO. 1 ) OELTsTF  - T 
C 

CALL  RUNGF ( T .4 .DFLT.COnS.7 ) 

C 

IF ( ILAST.Fo.1 > GO  TO  50 
5 continue 
50  COtJTINUE 
C 

C COMPUTf  oyCT  average 

C 

C r0NS(5»  = period  of  orbit 

CONS ( 30 ) :7 ( 1 )/C0NS(5) 

YAVr,rABS(7j4)  J/C0NS(5» 

0XAVG  = C0NS I 50 ) 

C 

C 

nELT  = CONS< 1 n » 

ILASTsO 

TzTO 

C 

C NTER'^S  = 5*(NU.  OF  COfFFICnNTS  TO  BE  FOUND)  ♦ 3 

C 

NTEP“S=35 

call  lCSf7,C0NS.NTER4S) 

7<2»=YAV3 

CALL  R0C<  n. .TONS) 

no  10  11=1.500 

iFtiT  ♦ DFiT)  ,GC.  TF)  IlAST=1 
IF ( tlast ,ro. 1 1 delTstf  - t 
C 

CAU  RUNSF) T.NTFRHS.OFLT.CONS.Z) 

C 

IF/ fLAST.ro. 1 > GO  TO  lOO 


BESrAVAIUBlE  COPY 


60 


65 


70 


75 


80 


85 


90 


95 


10  roNTlNUC 
inO  CONTINUr 

C TEHM  IC  THf  SC'ILC  FACTOR  Fqk  ThE  FoURIER  INTEGRALS 

TrR**=2  . /Cnr.s  ( 5 ) 

00  PO  I=‘*.NTFRHS 
?0  7(  I)=Zt  T ) •TER*' 

XCCrrFS(  1 »rOXO\/G 

YCCEFFSC 1 »=YAVG 

NPI»=(NTERMS  - 3)/3,  ♦ l.E-03 

no  ^0  1=1 .nfty 

XCnFFFS( 1*1 )=? ( I+3) 

YCOFFFSf 1*1  J =Z ( I ♦ 3 + NFIX) 

30  FCOrFFSE I )=7 ( I ♦ 3 ♦ NplX  + NplX) 

NXrMFiX 
NYrNFIX 
NF=NFIX 
WRITE (6. 9qq) 

999  FoR'"AT  (Hi  > 

wRITE(6.lnnOJ 

1000  F0RMflT(4?x.*  satellite  GROUND  TRACK  INTEGRATION*// ) 

WRITE (6. inn5lTC»  T 

1005  F0R'*AT('*8x,*  initial  THE  = *.Fin.4./48x.*  FINAL  T I ME  = * . F 1 2 . 4 / / > 
WRITE(6.lniO»CCN9(10) 

1010  FOPMAT(3?y,  « INTEGRATION  ST^P  SI7E  = (SATELLITE  PERIOD»/400.  =•. 

1 Fin.4»*  REC.*//) 

WRTTE(6« 1035) 

1035  FOR'^ATE  23x . • RESULTS  OF  GROii.gD  TRACK  INTEGRATION  FOR  THE  THE  IN 
1CRE'*FNT  3rLTA  s ( pER  lOn  ) / 1 00  . * /// » 

WRITE(6*iniS» 

I0l5  FOR'*AT(2lX.*  XG^il^X.*  YG*tlbx»*  F**  ;7X.*  PH  I • « 1 2 X « ^L  AMD  A • . 1 2X  i 
1 • TTME*/) 

DT=rONS(5)/inO. 

T = -nT 

DO  500  1=1. ini 
T=T  ♦ OT 

CALL  ROCd.COAS) 

CALL  FNTS9( T.CONR.XCOEFFs.YCOfFFS.FCOFFFS.F.NX.NY.NF) 
WRTTE(6.lnxrnF(H,F(2».F(4)  .CO.\)S(  14  ) .C0NS(  151  .T 
1030  FOP»*AT  ( 17X  , 5Fl  8 . 7 .F12.4  ) 

500  CONTINUE 

return 

END 
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20 


SURROUTI'JF  aNG<C«S*T.N) 

C 

c*««*«***«**»*««***««*«******«««***««**»«««»****««*«****»»*«t****«**»«*«»** 

C THIS  CO^*PUTFS  MULTiPLC  ANGLE  SINE  AND  COSINE  TERMS  FOR  THE 

C ARGUMEfjT  T s { P I E ) • < T T ME  ) / f HER  T 00  ) . 

C 

c time  r Current  time  in  the  program 

C PERIOD  = satellite  ORBITAL  PERIOD 

c«*******««**«***«««*«*««*«****«»**««****«*««»**«**«****«**««****«*«******* 

c 

oimensiom  r ( 1 > . s f 1) 

C ( 1 )=COS( T ) 

S(l»=SIN(Ti 

C (2)=C ( 1 » *0  I 1 » - S{ 1 ) «S( I I 
S(?)r2. *$<1  > *0 ( 1 ) 

TF(N.LE.?)RETLRN 

no  1 1=3. N 

SU)=2.*S{T  - l)*CU)  - - p) 

1 C(I)=2.*C(T  - 1)*CU)  - C(I  - P> 

RETURN 

END 
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55 


bO 


SUBROUTI^r  FNTs2rT.C0Ns*XC0EFPs,rr0EFFS*PC0EFFStPtNX,NY«NF» 


THIS  RROfiRiM  CO'^PUTFS  THE  SATElLITF  MAP  Pi  ANE  X.Y  COO«OINATES 
ANO  THF  F-FUNCTION  . INAODITION  ThE  1-ST  DERVATIVES  OF  XtYt  AND  F 

are  computed. 

A call  to  roc  for  the  time  T must  pRECEED  a call  to  this  program 

CONSMA)  = V(T>  AT  TIME  T Xu  THE  MAP  PLANE 
INPUT 

T s current  TIME  IN  PROGRAM 

x:ntFF«.NX=  X-FOURTER  fit  coefficient?  and  number  of  coefficients 

YCOFFFS.NYs  Y-FouRIER  fit  coefficients  and  number  of  coefficients 

FCOEFFs.NFs  F-FnuRIER  Fly  COEFFICIENT?  AND  NUMBER  OF  COEFFICIENTS 

OUTPUT 


C 

FM  » 

= 

X GROUND 

track 

c 

F ( ? ) 

; 

Y GROUND 

track 

c 

F<3> 

= 

tc-frahe 

rotation  angle  gamma 

c 

FJ4  > 

r 

F-INtEGRAL  oF  V(T)/fRA0lUS 

OF  CURVATURE) 

c 

Ox 

s 

DXG/OT 

, WHERE  XG  = 

X-GROUNO  track 

coordinate 

c 

nt 

= 

OYG/DT 

. where  yg  = 

Y-GROUND  track 

coordinate 

c 

DFJ 

= 

DF/DT 

□IMFNSION  XCOEFFsC 1 ) . YC0EFFS( l).FCCEFFS(n.Ftl) .C0NS( 1 ) 
DtMFNSION  C < » , ? ( 20  > 

COMMON/OF/ofT 
COMMON/nXY«AP/DX .OY 

Pir=?.*ASTN  n . » 

CALL  SETCOfJl  A.B.PIE.PFRIOD.COns) 

NFIX=NX 

IF( (NY/NX).GT.1»\FIX=Ny 

TIME  -UST  PE  NORMALTZFO  TO  T = P T E • T /PER  1 00  FOR  ANG  SUB  CALL 
TPrPlE*T/PFRiro 


CALL  ANGIC.S.TP.nfIX) 

coMPuTF  y-grouno  track  COOMUINATE 
Y = 0.0 

DO  1 1 = 2. NY 

1 Y=Y  ♦ TC0FFF?( I ) *0 ( I - 1» 

rOMPjTF  x-GROunD  track  COOhUINATE 
XsXCOEFFSf 1 j .T 
00  2 1=2. NX 

2 X = X ♦ XCOFFFSf  n •${ I-l » 

FT  = n. 

00  A 1 = 1 .NF 

6 FTrFT  ♦ FrOEFFS { I »«S(  I ) 


F <1  )=X 
F(?)=Y 
F(4)sFT 

TERMsPir-'PCRirC 

C COMPUTF  nx/OT  FOR  GROUND  TkACK 

nxscoNS ( 1 A j •ccs ( F m ) > 

C COMP'JTF  Oy/CT  for  ground  ThaCk 

OYsCONSJ 1 A> •?TN(F (4  M 


nF=n. 

no  A :=x.Nr 

8 OFrOF  ♦ I •rcoFFFS I n •C I I ) 
OFTrDF»T*'Rfc' 

F(3)=PIE/J>.  ♦ ATAf42lDY,DXI 

RETURN 

END 


BESTIVAIUBIE  COPY 
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25 


SUPROUTI'JF  srTCO\'(A.B.PIE.PE:RlOOtCONS) 

C 

c«*««****«*»**««««**«*««««**»**«*»«,««,*,*,»«, 

C THIS  sroGRAm  SCTS  VfiLUtS  OF  PARA'iCrERS  USED  IN  THE  CALCULATION  Of 

C THE  satellite  ground  TRACK 

C 

C NSTEpS  = Nl -RER  OP  STEPS  In  The  INTEGRATION  OF  S.G.T* 

C nELT  = TI-E  IncRE'^ENT  IN  EaCH  INTEGRATION  STEP 

c«»*««****«*»*«.*«««»*«««*««*«*««**,«, 

c 

DI^^ENSION  rONFl  I > 

F=1 ./29A,^ 

ESCsJ-'C"  . p,F 
A = f.37fl.  Ii5 
R=A*SQRT(i.  - ESC) 

Pir  = 2.*A3TNn  . ) 

PERTOO=10^.?ftT*60. 

NSTEPS=4nn 

OELT=PERTnr/N«TEPS 

CONS(6J=3rRIOC 

C0NS(7)=F 

CO^'S  I ) =A 

C0NS(9)=3 

CONS ( 10 ) rnELT 

CONS(l3)=n. 

C 

RETURN 

EmO 


1 


5 


10 


lb 


20 


?b 


30 


3b 


40 


4b 


bO 


bf> 


C 
C ' 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


c 


c 

c 

c 

c 

c 


c 

c 

c 

r 

c 

c 


<;ijf  .TOUT  T }f  I-!\;e:o«E  < TO  . T T .COnStJ-F  • X ♦ r . OHl  .LAVnA  ) 

thTS  c s A. 40  I AMOA  GTvEfJ  X AND  Y TN  THC  •lAP  PlANE. 

T'  PI  T 

X c A^SCIS’^A 

Y s D-«UIf.ATE 

ni'TP  IT 

n^T  s GCCOfTIC  LAyniPt'E 

(A^tia  s L0^^IT'ln^ 

TT  = TPF  DT<;pLACH“EM  vector  DFLTA  R LIES  ALONG  THE  SCArj 

VECTOR 


ro'  ''()  |/■|^"ATA/J'’«1T  .DPET  tCtH*  **  '^i 
rf'  -ir  i/Q-.V'/TE  ,rpoT  w .Ea-PT  .CT  . r 
•“O'  '*f’  i/"'- /I'ET 
rf'  •*O  j/''xx'  ap/Cx  .'■'Y 

ro^'^o  j/yfr<irr^/Y<*3EFf^,YroPEFL«FcorEF^»Ny»fiY«Nr 
CO'  "0  i/1  .ri-  T Al  / X'  «C»X  .«{  DXN 
ro*'vn  ./  A 'j-^/Ai  oj-AC 

Mf(»),''-4MJ«r^nT(3).rON«;ii)*ROT(3.x».xx(e,).-j(3’ 

IK  .n*.  f 3)  ♦ xf  M » .PX.^t  3)  .ODX'il  3).W(x».nwC3).TJ3).UTl3) 
nTi'rr;';lo'j  rcuJ.xrTEEF'scpoj,  YCtU  f F<;(;>n)  .ecoeffs<20J 
nt'r  'SliJ  r<Tl,T*?DT(3j,nRliT(f».E(3) 

ni'T*  •,  1 ov  TP  e ■* ) 

'/T I - f /T  tyE’ « J 

ri.  •’  • i.  "•’PA . * 

PEAl  L^'^OAO 
TOl 0=1 

TALL  -.t  TT -vf  { A . J .PTE  « f . rO'J«i  > 
r Ai  I .*uC  ( 0 . . r^\S  > 

r,r  T T^r  fiTTML  GuEcS  P'l  ll-'fc 

YrOL--<M>  .-.PPHr)!  Iv.aTFLi  the  average  velocity. 

TT  = x/  xL''"rr-^f'  ) 

I A’iTr  • 

TCf'-f  T = r 
no  Tn,i>;  TTrl  . t b 
p CAM  MU*  » T T .r*  ,*  ) 

rAii  »'*0''*Tri:j'“i,T0.Tl.*.E> 

VSTOKSi  . ■*  » 

TsC  J"?  M 7 > 

Y*  r 5r<-nT:-.  t4  I’irRTIAL  lmACF 
y = vM  :--TTY  I'j  IVfPTTAL  sPACF 
'"r  ,•  ; '.rr  t I r.<  AT  10’.  1'  INfHTl.lL  CPACE 

. *s  T..»  I-  T'<»  SCA,.  PIPfCTlOu 


c 


r\\  \ l.*-''-  «.<  »•  .w<‘  . 


I 


^7 


bO 


tib 


70 


75 


60 


66 


90 


95 


loo 


105 


no 


r 

TAIL  CiO':‘s  I v-i.  UOyfj^OU  ) 

C 

C 

C nu  TmJ:  '“CTtVATiVf  or  THF  SCAN  VFCTOR  IN  TH£  E-fRA^E  WITH  CO*»F'ONENTS 

C PI"«IT/V'.  «>PACC 

C 

OU<  1 S')  •'J(  2 » 

Oil  ci  )=D  I ( ? I - wr  •Lit  1 » 

01M  ^ ) =0  M ^ > 

C 

C 

C RCIAT:  THf  COvPONLMTS  qF  U.OU  TQ  ThE  E-FRA'IE- 

C 

CALL  lF  o A'-r  f ROT  . 3 J 
CAl  I Ef^Avp-  ( -iL  ,9rTi  3 ) 

C 

C T IS  TmF  iiMT  vector  along  TRACK  wiTH  COKPONEtJTS  TN  THE  E-FRAmE 

no  1 1:1.^ 

T » I > rOhOT I I ) /V 

1b  MT  I T ) zoOTPT  { T » /V  • OVOT^OrPTM  »/lV*Vl 

c r IS  THf  ll^IT  vector  T'J  THt  CROS*;  TRACK  DTRfCTIOfj  WITH  COMPONENTS  IN 

C THE  r-rOA***’. 

C 

CAf  f CRO«;c  < T ,C  » 

nPHT=CO'H;i 
DL/'  '*tJA=:  ( p7  \ 

cr.T  zC  US  Cof  s ( I ^ » 1 

SoTrSiNCjifiSe  !*♦  ) » 
n ( tzCGS  I rON^i  I 15)  ) 

<5Lr  TzS  ri  ( ro’jri !«)  ) 

c 

c ni.  Ic  T'-r  '*E"TvATlVf  or  THF  NORMAL  VECtOR  IN  THE  E-FRAMF 

C 

ONM  ) z-<;  ■ T *01  T •"PHI  - Co  r •SLf-T  *0)  AMOa 
••  T • '.t  " T *"PH  I ♦ Cpt^Cl-^  r*0i  amqa 
’>»  < T ) v-r *; T •npiJT 

c 

rn  ) c*‘''  = <^  < T •!.«■) 

r.H  i LM-'i;^  ( .P  » 

C 

C 'If  ic  T^r  ^lRTVATIVF  nr  THf  C VFCt^R  TN  THE  E-Fra»»F, 

C 

• ’ * 1 r 1 . ^ 

19  nc  I T ) =t  M ) ♦ ^ {}  \ 

C 

T1  SOI. , i j • r t 1 I ♦ rji?)*!!-**  • ini(5)»T(3) 

T'»r'(l)»*'Tl))  ♦l'l3)*rT|A) 

r^zMi)tr,i)  ♦ "j?)*!!?)  • riM«T»^i 

Tits''.  ♦ lUK^InCti) 

»‘(2l«rrFP»  ♦t.t3)*rc«<) 

TAs-.M)«*tT>  • U«2)»LIP)  • U(M*CM) 

C 

Tri.  I > * < T 1 . r 9 I . Ti,*J7**f5) 

••  * Tc»Tr^» 


r*  r 
JO 


115 


120 


125 


130 


135 


mo 


mb 


150 


155 


160 


165 


170 


AA  = T ^;/S 
PRsTf^/SJ'*«;0 
C 

c u 15  ThF  <^c■a^  vtCToR  AS  PHujEncn  ONTO  The  t,c  plane 

c 

00  iJb  1-1.3 

1*5  -J  ( I ) =AA«  T m ) ♦ pp  «C  ( U 
C 

OAr  ( T 1 + T?  > /S'l'^SO  - T3*TERM/(SlhtSQ*5i|V5;Q*SllMSO» 

PP=  I T.»  + Tt  » /S’1*'S0  - Tb*TEP'‘’/(Su.-1SQ*SUMSO*SUMSO) 

C 

C ni.'l'T  Tc  TH'’  ?rm\/ATlVF  OF  the  sCArj  VECTOR  W IN  THE  TC-FRAME. 

C 

no  si  I =1 , * 

61  n„rT ( 1 ) = ta* T ( T ) + aa*dt(T)  + oh^ccij  + bn*oc(i) 
c 

CAl  L CKOSc;  ( j.'a  , TO  j 

CALI  H.Ts?  ( TT  .CO^  vCOrFFs,  YCflEFFs.FcoEFFS.F.NX  .NY  ,NF  ) 

THFT ArF  ( 3 » -PTf /? . 

CTrrOSl  ’’HTTA) 

ST  = 5 t • U T ^r^  A » 

r)PYrDVi‘T*rn5l^  ( ^ ) ) - >'•  STN  ( F ( u n *nFT 
DnYs{)\/DT*<:T  JfF('J>)  * V*C05(F(4»mnFT 

nrv^F  TA=  ( Oy  •fiOY  - OY*ODx  » / { nx»nx  + nY*oYJ 

A1  rAA 
APrOH 
DAI  =t.A 
0A2-=D‘’ 

H1=-(A1»>T  f A2*rT) 

A1  ‘C  - A?  • 5T 
C 

c FT  = 1 nil..  T ft-P  1 . ( ri-CF05S-l«  > . FT  = 0.0  TI**E  = T-STAR, 

C 

FT=qi*{Y  - Ff1)»  * H2*<Y  - F(31» 

C 

* ^THFTA*A?1  ‘ST  ♦ »uA2  + HTHE  T A * A 1 ) •C  T ) 
np?T(DAl  - 1t6ETa*A2)*C7  - (Da?  * c TmLT A • a 1 1 • ST 
C 

nrFT=-LM  *Hv  - d2«nY  ♦ pPmix  - Fdii  ♦ Dn?*(Y  - F(2>) 

c 

IF  ( I AST  . ■ . . 1 r'''0  » 00  TO  20n 

c 

C FI  n'  ’IF'.  T^vr  rSTm/>TE  VIA.  NEwTOMS  *-FTHCn 

c 

TI  = TI  - -T/r)F'‘T 
C 

Tr(AMS(T-i||-TT).LT.l.r-OA)  LASI  = ir»P0 
rni  nr  r i 

0T  = T(ii : - r 1 

1 Ci'iO  rrr  t t I'f 

..rj  7 n ( ) 

3S2  F..P*'AT('  n.//»  failed  TO  COrJW^  rtbl  TN  iNVfRSE  T«STAP  SEAhCH*///! 

CAM  I A M 
^‘■0  rn'  T I I r 
T i-r  ?=I I 

C 

C MSI  Mi.i'.vja  .jf-  TMf  f.ROUNu  TRACK  AS  SiARTltjC  GulSSES  FOR  2-0  NErTqNS 


BESrAVAIlABLt 


Q 


C rMO")  fOR  T^r  PMl,  LAmHA  Uf-  POTNT  OFF  GROU^jO  THACk. 

c 

PMTOxCO*'^  M ^ » 

175  I A«nAO=C‘>f  «r  M “ ) 

np  tn  1=1.10 

f f;rri's  ( p H t n ) 

^ r *•  I ; M •■>  H T p t 
Csl  =COS<  ) 

180  RSI  r ■..•‘'HA'*  > 

k/l^Ti  SsA*//So®^<A«A*CF*CS  * B*H»SS*SS) 

c 

c r".<hr  iT  CftRTH  FixLn  cooROiijAir  csri'ioTrs 

c 

1 ^5  HI =pA"1 JS.rs.rsL 

H?=Pr  .)!  •'  .r  .«iL 

t«PAoT'JS«S*;/A/A 
C 

nwi  -p  1 - I » » ( 1 j - H*f.  n ) ) 

190  nRprP?  - cxxt^)  - 

HK^  = H ^ - J YX  ( ’ » - H*rj<  3 j ) 

’ tt  = ti  itr 

c 

C Al'JUST  TI'iF  FC«  SCAf,  RQTC 

195  C 

f.All  SCA4ir»MTI,Rl,R?.P5,tY(li,XXl?),yX(^),Nll),fg{2),N(3),H, 

1 r,u  ,0R2*rf-  5 .r  ( 1 ) .c  (?  > .c  r :')  I 

c 

C Rf  T I vaF  PLAfJt  X.Y  rSTI,.,aTFS 

200  C 

rn  t jSC  S.FmT''  .LAMRAO.  n.  « Xo.  YO.  lOfi  ) 

C 
C 

c ntr  .•’ApTrit  "rRiVATrvrs 

20b  C 

CAl  I nPFipr  I ft  , J .r  IfJS.rr  ,MP,K  .III  .UR.  V1  . VR.LAI^nAQiPHIO.  TI  ) 

on  Til  . / 1 Ml  • • VI ) 

c 

C rn  iPjTr  t^ri.  Rn  . LA^'OA  F ST  r .-lAir c , 

210  C 

PMlil'MlR  ♦ -If  1 T • ( ( X-XO  I *\/R  - rY-Yn»*U?) 

I A*  P/.SL  A ""lAO  ♦ OFL  r • I - ( X- YO  I *u  1 ♦ |Y-y0)*»H) 

IT  tArs{'’iI-PMTj).lF.  .l.F-OA.  Ar-in,  ARS  iLAMDA-lAMnAni.LC.l.E-OflJGOTOTO 
PMTnsi'HT 

2lb  I A'^'IA  .>slA'*i  a 

.•  0 rr.’  T r’ur 

PR  I n I u . 1 f 00 1 

1‘|  u rni.-^AriMl..  rAlLFO  CO^uFRul  T’''  I*KwKSF  PHl.LAMnA  SFAHCm*//) 

r . I I I X T r \\ 

220  fo  ro*  T I jUF  ! 

P/’T  f»-  1 !:: 

I 'if 


’l  > ' 

II  '' 


AO*AOd2  1A3 

unclassified 
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BEST’AVAIUBLE  COPY 


1 SUBROUTInjf  ICS(Z.C0N.S»MFIX) 

c 

c««** 

C THIS  PROGRfiF  SETS  ThE  INITIAL  CONDITIONS  FOR  THE  INTEGRATION  OF  THE 

5 C SATELLITE  GROUND  TRACK  ANO  THE  FOURIER  SERIES  COEFFICIENTS 

c*****«*««**«« 

c 

Dimension  z(i)tCONS(i) 

DIMENSION  XO(F) 

10  CrtMMON/ROCDATA/A.BiXOtTO 

c 

PIE  = ?.*ASTN( I . ) 

DO  •>  I = 1»NFIX 
5 Z(l)sO. 

15  C 

RETURN 

END 
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SUBROUTT'Jr  nrPIV(T.NfX.CONST«FI 

c 

c«*** •••*****«*••***«••«*•«•*•**••••****•***•*««*•*«*•«•«• ••••«***«••• 

C THIS  PROGRflf*  COMPUTES  THE  FUNCTION  APPROXIMATIONS  REQUIRED 

C TN  subroutine  RUNGE. 

€•«••«•*««•«•*•••*••*••«•****«****•**•**•«•••*••*•**«**••*•««•**••*••• 

C 

DtMENSION  xr  I ) .F(40) tCONSTC  t) 

OIMCNSION  C(9ft) »S(20) 

C CONST{i7>=COSnT 

C tONSTdftIrrsOT 

C C0NST{9r»s1 ./ROC  # WHFRE  ROCs  RADIUS  OP  CURVATuVE 

P=CONST(G) 

PIE=2t*ASTN(1 . ) 

F(  1 IsCONST  flftJ*CCS(X(3n 
Fl2  JsCONST  ( lfl)*SlN(X(  5n 
F(3J=CONST(lA»*CnNST<2A) 

C TP  N=4  wE  TNTPGRATF  TO  PINf)  DXAvG  AND  YAvG 

C 

IP(N.EQ.4)  F(4)=X(2) 

IF<N.E&.4>  GO  TO  lOO 
nXAVGsCONST(.BCI 
XTERMrXd  » - CXAVG*T 
YTERM=X(?I 
FTTRMsXIB) 

C M0RMALT7r  TI^'E 

TP=PIE*T/P 

C 

C NTERMs  ALlCwS  US  TO  VARY  ThF  NUMBER  OF  COEPPICIENTS  TO  BE  PIT 

C 

NTERMSs(V  - 3»/3.  ♦ l.F-n3 

NPIXsNTERMF  ♦ 3 

CAIL  ANG(C.S.TP,MERMS) 

DO  5 I = 4tfjFIX 

FH  )=XTER^'*Sf  I - 31 

F<I  ♦ NTEPMS)=YTERM*C ( I - 3» 

5 P<I  ♦ NTER»*S  4 NTERMSlaFTERM**;!  I -3) 
inO  CONTINUE 
RETURN 
END 
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subroutine  R0C(T»C0NS> 

C 

c«*»*««* ••«•««•*«**«••**«««»«*•••*••«•«••••••••••••«••*••••••••••••• ••*•••• 

C THIS  PROGRAH  computes  the  RAOlllS  OP  curvature  ruo  THE  first  and 

C SECOND  DERIVATIVES  OP  THE  MoTIqN  ALONG  THE  SATELLITE  GROUND  TRACK, 

c************* 

C 

DIMENSION  CC( ! ) . N < 3 » .ODROT < 3 ) , DORDS ( 3 » . ORDS ( 3 ) , DRDt ( 3 ) 

DIMENSION  V ( «0( 3) tUI 3) ,M( 3«3) tUVl 3) •WWR(S) , OhDT ( 3) .ODHOTt  3) 
DIMENSION  V0(F)«X(6)«C0NSC1) 
dimension  XN(* J ,nXN( 3) .DDXN t 3» 

COM*^ON/ROrnATA/A  .8*X0,T0 

COM»«ON/NSnATA/OOROT,ORnT,CC*H,x.N 

COM'^ON/INEPTAL/XNfDXNtOOxN 

real  LAMDft.fi,»'U«N 
MUsSRfiSOl .? 

C 

C 

c 

CALL  ORBlT(Xn«X.TO«TiCONS» 

C 

C XN  = IS  the  satellite  POSITION  IN  INLRTIAl  SPACE 

C DXN  s IS  the  satellite  VELOCITY  IN  INERTIAL  SPACE 

DO  13  Isl.3 
XN( T )sX( 7 ) 

13  OXNin  = X(T  + 3» 

DO  10  1=1,3 

10  V<l)=Xn*3) 

C 

CALL  ROTATE|MfTO*TfRATt » 

C 

C 

C RrRITl,  TsTIME  « RSQsR ( T > •R I T » 

C 

RSOsXCl)*xm  ♦ Xl2>*X(2>  ♦ X(3J*X(3I 
RrSORT(RSO) 

R3=RSO*S3RT(R«C» 

C 

C U TS  THE  ANGULAR  VELOCITY  OF  ThE  EARTH  g.R.T*  INERTIAL  SPACE 

C 

W(1  )=0, 

UC2)sO, 

W(3):RArE 

C 

C 

CALL  CR0SS(W.X,O) 

C 

C COMPUTE  ThF  E-FRAME  VELOCITY  WITH  COMPONENTS  IN  INERTIAL  SPACE 

C 

DO  A 1=1.^ 

6 DROT ( 1 ) =V( I ) • DC  n 
C 

C rOMPjTF  THE  E-FRAME  VELOCITY 

C 

CALL  EFPA*'rinRUT,M,3» 

C 

C WV  TS  THE  CORIOLIS  ACCLLERATIQN 
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C 

c 

CALL  CR0S^fU«V<UV) 

C 

C wyR  IS  THE  centripetal  acceleration 

c 

CALL  CROSSrW«C«UWR) 

C 

C ODROT  IS  THE  ACCELERATION  OF  THE  SATELLITE  IN  E-FRAME  WITH 

C COMPONENTS  IN  INERTIAL  SPACE 

C 

00  11  Isl.s 

11  OORnT< I )s>Mu*y ( 1 )/R3  - 2,«wv(l)  ♦ UUR<I) 

C nOXN  = TS  the  satellite  acceleration  In  InErTIAL  space 

DO  15  1=1. S 

15  nOXN(I)=-MU*XA(I)/R3 

c 

c COMPUTE  THE  E-FRAME  COMPONENTS  OF  THE  ACCELERATION 

C 

CALL  EFRAmE (nrRDT.M.3) 

C 

C COMPUTE  the  E-FRAME  COMPONENTS  OF  THE  POSITION  VECTOR 

C 

call  EFRAME<X.M.3» 

C 

OOWN=X(l ) «X ( 1 ) ♦ X(2)*X(?) 

LAMnAsATAN^(Y(2)tX<l) ) 

0LAM0A=( X ( 1 )*CRDT(2)  - X ( ? ) •OmdT ( 1 ) 1 /DOWN 

OOLAMDAsIvt 1 )*00R0T(2)  - x ( 2 > «HORnT M M /DOWN  - 2.*0LAM0A« 

1 (X<l)*nRDTni  ♦ X<2|*DkOT(2I  )/n0WN 
C 

RXYsS0RT(X<l)*X(t)  T X(2)«X(2n 
C 

CALL  PHIHfPXY.Xd)  tRtPHi  iHiLAhUAt  A.R) 

C 

CALL  OPHIOH( A»S*H.HP.HpP,pHl.UHHl.OOPMl.X.RXY.R.MU.OROT»DDROT> 

C 

C=COS(PHI > 

SsSTN(PHI ) 

CLsCOSILAMOA) 

SLsSlNJLAMnA) 

C 

C COMPUTE  E-CQHPONENTS  OF  DH/OT 

C 

OHDTriisHp*c*rL  - H*(s*CL*OPHr  ♦ c*sl*olamoa> 

0MDT«2»sHP«C*SL  - H#<S*SL*0PH1  • C*CL*OLAMOA) 

OHPTt 3) =HP»S  ♦ H*C*OPHI 
C 

c rOMPUTr  E-CCMPONENTS  OF  02m/0T? 

C 

OnHnTJl»sMpP*C*CL  - 2.»HP*(S*CL*0PHI  ♦ C*SL*0LAM0A) 

1 ♦H*r-C*CL*<DPHl*nPHl  4-  01  AMOA^OLAMOA  ) ♦ 2.*S*SL*DPhI 

2 •ola*»oa  - s*ri*DnpHi  - c*sL*nOLAMOA) 

nnHnT(2)sHPP«c*sL  - 2«*Hp*(s«sL«npHi  • c*cl*dlamoa) 

1 ♦H* ( .r tSL « (OPHT •OPHi  ♦ OL AMOA*OlAmOA I •2**S*CL*0PHI 

2 •DLA^'CA  ♦ C*CL*ODLAMUA  -S*SL*OOPMn 

n0HnTj3)sHPP*«  ♦ 2.*HP*C*0PHI  ♦ H*JC*0DPHT  - S*OPHI*DPHn 
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C 

C OROT  fq  THC  VCLOCITy  OF  THE  SATELLITE  6.T.  iN  THE  E*FRAHE. 

00  4 

9 OROT(nsOROTlT»  • OHOT 1 1) 

C 

C COMPUTE  THt  magnitude  OF  THE  VELOCITY  ALONG  THE  SPACE  CURVE 

C 

SUMsO. 

00  I?  l:l«3 

12  SUMsSUM  ♦ DROT(n*ORDT(I) 

OSDTsSQRT(SUM) 

C 

C nOROT  Ts  THE  ACCELERATION  OF  THE  SATELLITE  G.T.  IN  THE  E-FRAME. 

C 

00  16  1=1.^ 

16  ODRDTI  DsnOROT  ( I ) - ODhOTUI 
C 

C COMPUTE  THE  COMPONENTS  OF  THE  UNIT  VECTOR  AlONC  THE  SPACE  CURVE 

C 

DO  17  1=1.3 

17  OROSM  »sORnT<  I l/CSOT 
C 

C CONST  Ts  THE  DOT  PRODUCT  OF  THF  VELOCITY  AND  ACCELERATION 

C 

CONSTsDORDtII »*ORDTIU  ♦ ODROt t 2 1 •ORDT « 2 I ♦ DOROT C 5 ) •OROT ( 3 1 
vso=osoT*nsnT 

c 

C Nin.Isl.S  ARE  the  COMPONENTS  OF  THE  UNIT  VECTOR  FROM  THE 

C THE  SATELLITE  GPOUNp  TRACK  TO  THE  SATELLITE. 

C 

Nil  )=C*Cl 
NI2)=C*SL 
Nl3lsS 
C 

CALL  CROSSIN.CRDSiCC) 

C 

C R0SQ=1./(RCCI**2. 

C 

C here  JF  make  use  of  the  fact  the  CC  is  PERPENDICULAR  TO  OROT 

C 

TERMsIOTRDT  M >*ccn  ) ♦ D0r0T<P»*CCI2»  ♦ OOROTI3»*CC(3U/VSO 
ROSOsTERM^TER*' 

c 

SUMsO. 

00  19  1=1.3 

19  SUMsSUM  ♦ nRnT(I)*D0R0TIIl 
DOSnTsSUM/nsnT 

c 

CONSI 14 ) =PHT 
CONSI 151 =l AMDA 
CONSIlTlsnOSDT 
CONS  I 16  1 =n<J0T 
COf.'SI  191  sROSO 
C0NSI?61=nPHt 
CONSI ?71 sOL  AMCA 
C0NS(?8IsTFRM 


RETURN 

50  FOR’^AT  I //I  0 K . T5  1 
END 
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suBROumr  rotate(m,to»t.we» 

DIMFNSION 
REAL  M 
C 

C angular  rate  as  radians  PEr  mean  sioerial  day# 

c 

PIE=3.14159?A5 
WE=?.*PIE/A61 64. 09054 
THETAOsO. 

THETAsThETAO  ♦ WE*(T  - T0> 

C=C0S(THETA» 

S=SIN{THETA) 

DO  I I=lf^ 

DO  1 J=l»3 
1 M(1.J)S0. 

M<1.H=C 

M(1.?JSS 

M(2.l)=-S 

M(2.2)=C 

M(3.3)=l. 

RETURN 

END 
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SURROUTI>Jr  PMlH  lRXt»2.RN»PHl  iHtLAMOA.  A.6) 

C 

C THIS  3ROGRAR  COMPUTrS  PHI  Aft»0  H BY  A 2-0 1 ^•ENS I ONAL  TAYLOR  SERIES. 

C 

c RXY  = 1X*X  ♦ Y*Y)**0.5 

C 

R = A 

PHI=ASINC7/RN) 

HsRN  - R 
C 

DO  1 i=i « 1 n 
CsCOS(PHI » 

S=STN<PHI > 

TERM=A*t«C*C  ♦ 0*B*S*S 
TERM?rA*A/SQRT (TERM) 

q=A*A*(A*A  -R*Q ) •SIN(2.*PHI )/9./(TERM**1.5) 

Ell20*C  • (TFRP2  ♦ H)*S 
Fl?rC 

E?1=B*B*(0*S  ♦ TCRM2*C»/A/A  ♦ M*C 
E222S 

nELT=Ell*EPP  - E1P*E21 
OTFPI2RKY  - tTER*2  ♦ M)*C 
OTFPpsZ  - <R*B*TER«2/A/A  ♦ Hi*S 
PHTMsPHI  ♦ (nTFFl*E22  - 0 I FF 2*F 1 2 » /DEL T 
HN=H  ♦ (-nTFFl*E?l  ♦ 0IFF?*E11 )/DrLT 

TF(  ARS(PHTrj-P»-I  ) .lE.I.E-OA.  ANn.Afl«;iHN-H)  ,lE.  1 .E-O6I6O  TO  10 
PHlsPHiN 
1 HSHN 

WRITE(6«5t 

5 format* • ,//•  failed  TO  CONVERGE  7N  PHIH*» 

STOP 

10  PHlsPHlM 
HsHN 
RETURN 
END 
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SUBROUTiNjr  CRCSS(AiO»C) 

Dl^rNSlONJ  ) |R(1)*C(1) 

C 

C THIS  PROGRAM  COMPUTES  THE  CKOSS  PRODUCT  OF  TWO  VECTORS. 

C 

C<1  - A(3)*B(2) 

C(2»=A<3)*nM  ) - A(1)*R(3) 

C ( 3)SA { 1 1 tRl P)  - A(2)*P<1> 

RETURN 

ENO 
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SUPROUTl^Jr  EFRAME(A,B.N» 

THIS  program  rotates  COORDINATES  BETWEEN  DIFFERENT  FRAMES 

dimension  a ( 1 ) *8 ( 3» 3 ) 1 C ( 5) 

DO  I 1=1. M 

1 C<  I)=A( I ) 

DO  ? I=1»N 

2 A < I ) =B 1 1 » 1 ) *C ( 1 » ♦ B(I,2)*C(2»  ♦ B(I«3)*C<3) 

RETURN 

END 
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SUBROUTINir  ORPIT  ( X , XF  « T I ♦ TF  ♦ CONS  ) 


C 

C THIS  PROSRAM  CETCRMINES  THE  STATE  HISTORY  VIA  THE  CHANGE  IN 

c rccFfjTRic  anomaly  solution 
c 


c 

c 


c 

X ( I ) 

Is 

THE 

initial  state  vector 

c 

XF(  T ) 

IS 

THE 

FINAL  STATE  VECTOR 

C 

RO 

IS 

tmF 

INITIAL  RADIUS 

c 

VA 

IS 

ThE 

INITIAL  VELOCITY 

c 

RN 

Is 

THE 

Current  radius 

c 

A 

IS 

the 

SEII-NAjOR  AXIS 

c 

n 

IS 

THE 

DOT  PRODUCT  OF  POSITION  ANO  VELOC I T Y ( I Nl T I AL ) 

c 

c 


c 

dimension  X ( 1 ) .XF(  1 ) 

Or^rNSIO'4  CO'}St  3 ) 

real  mu 

CONS{20)=?.*ASIN{1. » 

CONS(l»  =3R®601.2 
MUsC0NS( I ) 

PIEsCONSCPO) 

C INITIAL  state  (INPUT) 

C 

RO=X(l)*Xin  ♦ X(2)*X(?)  ♦ X(3)*X(3» 

V0=X14)*X(U)  ♦ X(5)*X(«i)  ♦ X(^)»X{f>) 

C 

R0=SQRT ( TO) 

VO=SORT(Vn) 

APsP.O/RO  - VC*VO/MU 
IF(AP  .Li . n.''  ) GO  TO  500 
Azl .0/AP 

PPP= (2»0*PTE*fl*SORT(A)) /SORT { MU ) 

C 

D=X(l)*X(i4»  ♦ X(?)*X(5)  ♦ X(3)*X{A) 

C 

C0MS(2)  =1.0  - RO/A 

CONS( 3>  =n/sCRT(MU*A) 

C0NS(i4)  =SORT{''U)/(A*sQRT(A)  » 

CorjS(5)==»pr 

YYrrONS(4> 

T = TF 

nElTzT  -TT 
nrLT=AMOO(PEl T.PPP) 

5 TsDFLT 
C 

CALL  NE WTONI CCNS tPHI . TI . TD 
C 

SsSTNIPhi ) 

C=COS(PHI » 

F=1 . n - ( A/RO ) • ( 1 . 0 - c ) 
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G=(T-TI)-M.n/YY)*(PHI-S) 

SUMrO.O 

60  C POSITION  AY  TIME  TF 

00  140  Isi  .5 
IP3  = l4-3 

XF< I )=F*X ( I ) + G*X ( IP3) 

140  SUMzXF ( I )*XF( I ) ♦ SUM 

65  RN=SQRT ( SUM) 

FD=-(SQ^T(MU*A1/(R0*RN) )*S 
GOn  .0-(  A/RN)  • (1,0-0 
C VELOCITY  at  time  TF 

DO  145  1=1 .3 

70  lP3=I+3 

145  XF  ( lP3)=Pn*XM  ) ♦ GD*X(IP3) 

c 

GO  TO  lonn 
500  wRiTr{6,3nn) 

75  300  FORMAT ( H1 ) 

WRITE  ( 6.  ?®.n  ) 

250  FOftMAl ( ///52H  error:  HYPERBOLIC  EXCESS  SPEED  ACHIEVED  IN  ORBIT 
1 ) 

STOP 

aO  1000  CONTINUE 

RETURN 
END 
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SUPROUT  I 'JF  NFU  TON  ( CONS  « PHI . TI  t TF  ) 


C 

c*** •«*****«*••**«•*««•*«**•*•««**«««•«*•*«•««*«*»«******•«•«•*•••***• •**•** 

C THIS  PROGRAM  USES  NEWTONS  MfTHOO  TO  iTERATIVELr  SOLVE  FOR  THE  ANGLE 

C PHI  THAT  appears  IN  KEPLERS  EQUATION. 


C 

OIMENSIOV  CONS ( 1 ) 

YPI=C0NS(4>*ITF-TI» 
nELsAMOO(YPI.2,*CONS(20) ) 

C STARTING  VALUE  FOR  NEWTONS  METHOO 

PHirDEL 
C=CnS(PHt ) 

S=SIN(PHI ) 

EPSIL0Nsl.F-l2 

no  100  I = T . in 

U=C0NS(2» 

?=CONS( 3) 

YcCONStH) 

phtn=ph!-ipht-w*s+2*(i.o-C)-  del  ) / ( 1 . 0-W*C+Z*S ) 

IF  I Ans( PHTM-PH ) .LT.  EPSILON)  GO  TO  lOl 

PHlsPHiN 

CsCOS(PHI 1 

S=STN(PHI J 

ino  continue 

WRITE (6* in > 

10  FORMATdHi.*  Failed  to  convekge  in  subroutine  newton*//) 

STOP 

101  continue 

Pu*  tPHiN 


' r "JN 
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(>0 


65 


SUPROUTI>jr  DPhIOH(A.B.H»HP.HPP»PH!.DPHl.ODPHI.X*R)fY«R.MU»OROT* 

1 OORDT) 

niMFNSION  X i I ) tORDK 1 ) .DOROT (I > 

RCAl  LA'IOA.MU 

C THE  PRORRar  COMPUTES  THE  FIRST  AND  SECOND  DERIVATIVES  OF  PHI  AnO  H, 

C 

c RXY  s (X*X  ♦ Y*Y)**0,5 

c R = (X*X  ♦ Y*Y  ♦ Z*Z)**0.5 

c PRXY  r D(RVY)/DT 

C ODRXY  = n?(RXY)/0T2 

c 

ORXYsC  X(  n*ORCT  1 1 » ♦ X ( 2 ) •ORDt  » 2 ) J /RX  Y 
w=x f 1 ) •ORnT J ?)  - X(2)*0R0T(1> 

nORXY  = W*J/(RXY*RXY*RXr)  ♦ ( X ( 1 J •OOflOT C 1 I ♦ X ( 2 I •OOROT ( 2 ) > /RXY 
OZ^ORDT 15^ 

DOZsOORTTtSJ 

C 

CsCOSfPHt ) 

S=STN(PHI ) 
aNG22.*®HT 
C 

C COMPUTE  OPRI/OT  AND  OH/OT 

C 

TERM5A»A«r*C  •*  9*B*S*S 
TERt^PsA^A/SQRTiTERMl 
ETAsTER'*?  4 H 
rTA?5B*9*TFR»*?/A/A  4 H 
C/^N  = A*A  - R*R 

Qsa*a*C0V*9IN(ANr,  »/2./(TERM**1,5t 
All=0*C  - rTA*S 

ai?=c 

AP35B»B*3*S/A/A  4 ETA2*C 

A22  = S 

ClsDRXY 

C25DZ 

DELTsAll*a?2  - ai?*A21 
C 

0PHlslA25*Cl  • ai2*C2)/0ELT 
Hp5i-a2i4ri  4 aii«C2)/dElt 
c 
c 

C COMPUTE  DP(PHI)/OT2  AND  D2h/0T? 

nQ=A»A*CDN«(rrS*ANGI/(TERM**l.SJ  ♦ 3,4C0N*SIN(ANG1*SIN(ANGW4,/ 

1 i TCRM**2.«.  > ) •CPHI 
DrTarQ»DOHT  ♦ hP 
0CTa?=b*R*0*0FHl/A/A  4 HP 
Rn=0*C  - FTA^S 
B1?  = C 

B?1  28*0*‘5*«;/a/A  4 ETApmC 

fl?2=S  . , 

OlrODRXY  - nQ*CP^^I*C  ♦ 2 . tOE  T a •$  •PPH  I 4 ET  A •C  •OPH I ‘DPH I 

^ D?  = nOZ  - H*n*rG4''PHI*S/a/A  -2 . •DETA2«C»0PHI  4 rTA2*S»DpHl»0PMl 

c 

nELT?=Bl 1 •R?>  - Pl2*B2l 
C 

c 

OOPhT =t - B12*02 > /DELTZ 
HPP=I-B21«D1  4 Bl 1*02 » /0ELT2 
C 
C 

RETURN 

END 
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SURROUTINF  VFCPR0(A.B 
DIMFNSIO’^  an  ) »B{1) 
|F(NU.F>'3.niGO  TO  S 
THIS  RROGRAP  FORMS 
C = 0. 

00  1 1=1*^ 

1 C = C ♦ A( I ) •Bt  T ) 

RETURN 

5 C = 0. 

DO  > I=l»T 

2 C=C  ♦ A ( T ) .A  < I ) 

RETURN 

END 


C ‘NU) 

THE  INNEK  PROOUr.T  OF  TWO  VECTORS# 


BEST'AVAILABLE  COPY 


1 


5 


10 


15 


ao 


25 


30 


35 


<;UR»nuTI'jr  sr«N\rR(T»Ri«R?,R3.Xl#X2.x3«NliN2.N3«H#DRliDR2»DR5« 

1 Citr2.r3) 

c******************* 

c 

C THIS  3R0nRflH  COMPUTFS  A LINEAR  CORRECTION  TO  THE  TlwE  WHEN  THE 

C SCAN  ^ATF  CAN  NOT  BF  CONSIOFREO  INFINITE 

C 

C T 

C H 

C PI .R2.PX 

C Xl tX2»X^ 

C 
C 
C 

C««»**«********«**«« 

RFAi  i'll. N?. NX 
S1SR1  - *1 
S?=RP  - 
S3  = PX  - t^^ 

S=SnHT{Sl*^l  ♦ S9*S2  ♦ S3*S3> 

c 

‘?nOTHs-“»*  ( R1  •A  1 + S2*NP  ♦ S3*N.X) 

C 

nRnoTc=o^i *ci  ♦ pr2*c2  ♦ nR3*c' 
c 

TF  .FO.O  GO  TO  *> 

SIGNrnH'^DTr/APS<rROOTC  ) 

5 TFIORDOTC.Frj.n  »STGN=1  . 

TER'^sSOOTh/  < ) 

F=Ar0S(T'RM>*«lGN 
OTrO.OlSX*;*; 

FHAXsO.inOAA 
nT=nT*lE/F«A<i 
C 

T=T  4 DT 
RFTijRN 
FNP 


TIHE  As  FOUNO  in  NS  ASSUMING  INFINITE  SCAN  SPEED. 
height  of  SATELLITE  ABOVE  FARTHS  SURFACE 

mapping  pt. 

EARTH  fixed  RATFLLITE  POSITION 
NORMAL  VECTOR  TO  EARTHS  SURFACE 

OlSPlACEMENT  OF  MAPPING  PT.  FROM  SATELLITE  G.T.. 


^^SFAVAIlABircOPY 
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The  following  are  the  program  output  for  six  cases  and  four  sub-cases 


Case  // 
1.0 

*2.1 

*2.2 

*3.1 

*3.2 

*4.0 

**5.1 

**5 . 2 

**6.1 

**6.2 


* 

Assumes 

** 

Assumes 


Description 

Satellite  groundtrack  coordinates 
and  Fourier  series  coefficients: 
(x  , yg)“  plane  coordinates 
(I, A)  = ellipsoid  coordinates 

Forward  transformation: 

->  (x.y) 

for  (})  = 1.155860  radians 

A = 2,375706  radians 

Forward  transformation: 

for  (j)  = -.002888  radians 

A = -.122576  radians 

Inverse  transformation: 

(x.y)  ^ (^.A) 
for  x = 10099.660  km 
y = -59.587  km 

Inverse  transformation: 

for  X = 37787.170  km 
y = 791.363  km 

Distortion  error  analysis 

Forward  transformation: 

(4i,A)  = (1.155869,  2.375706) 

Forward  transformation: 

(4), A)  = (-.002888,  -.122576) 

Inverse  transformation: 

(x,y)  = 10099.59,  -59.58) 

Inverse  transformation: 

(x.y)  = 37787.67,  791.390) 

Instantaneous  scanner  sweep. 

constant  scanner  sweep  rate  of  314.9  deg/sec  (13.62  Hz). 
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SPACE  OBLIQUE  ‘CRCATOR 
SATELLITE  GROUND  TRACK  PFKMECTION 


CASE  i.O 


INPUT  parameters  nrPlNING  THE  CURRENT  RUN. 


tion 

velocity 

0 0 0 n 0 

7.3936 

601SU 

0.0000' 

47443 

0.0000' 

ORBITAL  ClEHENTS 


7291 .2m7l 
•00000 
1. 727Ae 
3,l41*i9 
1.^7080 

o.ooono 

6196.02000 


SEMI-major  axis  of  the  orbit 
eccentricity  op  the  orbit 
inclination  Of  the  Orbit 
LONGiTunr  OF  THE  ASCENDING  NOOE 

argument  of  peRifocus 

TIME  OF  PERTFOCAL  PASSAGE 
SATELLITE  ORttiTAL  PERIOD 


EARTH  RELATED  PARAMETERS 


flattening  factors 

SEMl-MAJORr 
semi-minor  AXISS 


.3S523E-02 

.637fl2E*04 

.6356aEf04 


BEST'AVMLE  COPY 


TABLE  D1 


SPACE  OBLIQUE  MERCATOR 
SATELLITE  GROUND  TRACK  PROJECTION 
CASE  1.0 


initial  times  0.0000 

FINAL  TIME:  6196.0?00 


results 

integration  step 

OF  GRCuNC  TRACK 

SIZE  s ISaTELLITE  PERlOOl/400.  = 

integration  for  the  time  Increment 

15.4901  SEC# 

delta  s iperiodi/ 

loo# 

Xb 

YG 

F 

PHI 

LAMOA 

TIME 

0. 

.9lF3S46t«03 

0. 

.lH14Cil9Ct0l 

.l570796E4ni 

0.0000 

.9145979E+03 

-.A935lfi3E-02 

, 140268?E*0l 

.1183A95E*R1 

61 .9602 

•0O95l98r^n^ 

.9O9i349E*03 

-.l7e3SJ6E-0l 

. 137H0pE*01 

.Ae24994E4nO 

123.9204 

.12l273?r*04 

.9n0l3fiflE*03 

••9A66900r*0l 

.132741)3f401 

,6732978E4nO 

lfi5.Aen6 

•l61SAB7r«04 

. 8A756AAE  + 03 

-.3S58993E-01 

.1276930E*01 

.529lAfe4E4riO 

247.8408 

•2020969E*04 

.ATISNOIE^OS 

-,4A979AIE-01 

.1222606C401 

.42S1114E400 

309.A010 

,2424959E*n4 

.AA?053AE+03 

-.S238A2PE-01 

.llbS9<i'iE401 

.3491731E+nO 

371.7612 

•2826A4Sr«04 

. A?«20«7E*03 

-.An594A3E-0l 

. 1107843E*<!1 

.28931 90E4«0 

433.7214 

•3232Ai?r4n4 

.8nT098PC*03 

••AAS6P49E-01 

,1048721E*01 

•2410a02E*00 

495.6616 

•363695nE«n4 

.7737973E+03 

-.7A2An37r-0l 

.988891tlE*00 

.90l0i92E+n0 

557.6418 

•403974Ar«n4 

.7414343E*03 

-.a:^657AAe-01 

.928539iE4'’0 

.1668938E400 

619.6020 

•4443l0nE*04 

.7n6l842E*03 

-.9fl729AnE-0l 

.8677883E400 

. 1 37lA98E*nO 

681.5622 

.4646302E«04 

1256E+03 

-•9743A20C-01 

.8067235E»1)0 

.11085l9E*n0 

743.5224 

,5299350E*04 

.ft?742A3£*03 

-.in37A25E+00 

.7454054E4I1I) 

.87l2l77E-ni 

605.4826 

.S652?49r^n4 

.5A4292SE+03 

-•in9A797C400 

.8838790E4l)0 

. A5443AIE-01 

867.4428 

•60549B9E*n4 

76flPE*03 

-.11 51A43E400 

.822l785E*00 

,453975lE-'^l 

929.4030 

,6497**67r*04 

.491 134AE*03 

••1201947E400 

.56033I1E400 

.268Gn74E-ni 

991,3632 

.666004AE«04 

.4416000E+03 

-.1247A07E+00 

.498359iE4')0 

.A979925E»R2 

1053.3234 

,7262S76r*04 

. 3902995E*05 

-»1?68144C+00 

.43628l2t»00 

-.78**7a54E-02 

1115.2836 

• 7664'>BAr404 

. 3374594E+03 

-•1 323A95E400 

.3741139E»I)0 

-.739flA21E-nl 

1177,2438 

•0OS6A94r*n4 

. >A3277r»E403 

- , 1 354O?0r400 

.3118720E*0t) 

•.395792GE-ni 

1239.2040 

•a4687iir«04 

.PP797A9E*03 

-.137AO97E400 

.249S892E*(10 

-,5473a44E»oI 

1301.1642 

•ea70AS4r«04 

.171777aE*03 

-.lA9RA2ftr*00 

.1872184E4°0 

-,A956799E»nl 

1363.1244 

.9272^40E*04 

,HMa973E  + 03 

-.1412A34E400 

. 1248319E4H0 

-.84i59i2E-nl 

1425.0646 

.96743eAr«04 

.57«62«2E*02 

-•1420961E+00 

.8242184E-01 

-,9859A52E-ni 

1467.0448 

• 1 007A2?r*OA 

-.P9Pi4?^E-0fl 

-•1423773E400 

.8353353E-13 

- . 1 l29A52E4nO 

1549,0050 

.l047A0«>r*05 

-.'S7*6?^2E*02 

-.14209AIE*00 

-.4242184E-I)1 

-.1273l52E+nO 

1610.9652 

,iO079e9r*nA 

-.n‘‘8’*73E  + 03 

-.»4i?s34E400 

-.124S319E*00 

-.1417A14E400 

1672,9254 

•112917Ar«0A 

-.171777SE+03 

-•1398^?8E*00 

-.1872184E*00 

-.1563425E4n0 

1734. A856 

• 1 166379r«0A 

-.2p797A2E*03 

- • 1 A78997c*00 

-.2495b92E*00 

-.1711720E4''0 

1796, A458 

. 1208974r*04 

-.?A»277nE403 

• • 1 ^54n90E  *00 

-.3118720E»00 

-,18b3312E4O0 

1658.8060 

, 12‘*878Sr»''^ 

- . 3374554E*05 

-.1 A23A95e*00 

-.374i139E»I)I> 

•.?0l9225E4«n 

1920.7662 

•i289n0Ar«0A 

3902^93E*03 

-• 1 28Al44r«00 

-.43428X2E4dO 

-,2180a36E4^0 

1982,7264 

,i329P39r*n4 

••44t6000£^03 

-.1 P47A07r*00 

-.498359ie*00 

'•.2348904E4n0 

2044.6866 

.1369489E*05 

•.491 i346£403 

-#120194 7E*00 

-.580331lE»00 

-.2525712E*00 

2106.6466 

,1409745E*ft5 

-.93P7682E+03 

-.1 151A43E*00 

-.6221789E*I>1) 

-.27l3nA0E4O0 

2168.6070 

.IHSOni^r^OA 

-.5A42529E405 

- • 1 n9A797E400 

..6838790E*I>0 

. -.?9l3941E4nO 

2230.5672 

,1490^0Ar*0A 

-.A974283E403 

- • 1 n37A25r*00 

-.7454054E*00 

“.3l303?5E*f'0 

2292.5274 

•l550Al3r*O9 

- . AaA 129AE403 

-,9743A?0r-0l 

-.8067239E*<>1> 

-.3367ai9E400 

2354.4876 

.l57093ir*OA 

-.7O€l842C^03 

-.9fi7?AA0r-0l 

-.8b77883E*00 

-.3G31O02E4Q0 

2416.4476 

• i6ii9S9r4>nA 

74l4943E*03 

-.A3697A8E-01 

-.928339ie»00 

-.3928o42E*00 

2478.4080 

•issiAiAr^OA 

-.77*7973E*03 

-.7A26n37r-oi 

-.9888910E*00 

-.4269>98E*00 

2540.36A2 

.lS9l98?r*0A 

-.8n«08A?E^05 

-.fcM5A949E-0l 

-.t04e721E*0l 

*.4669907E400 

2602. 32A4 

•l732359r*nA 

-.Ape20?7E*03 

-.fcn594A3E-0l 

-.110784SE»01 

-,5152295E4R0 

2664. 28A6 

,l772747r»0A 

-.AA2053AE403 

-•9238A22C-01 

-.116596'>E*01 

-.5750a36E*O0 

2726.24A8 

.U13147E«05 

-.8715«01E*03 

-•4A979felE-0l 

•.1222806E401 

•.8520pl9E*nO 

2788.2090 

.I853959r*09 

..8A7588AE4-03 

•.  ^S36993E-0I 

-.127fc930E*01 

-.75504A9C4nO 

2650.1692 

•l89397nr*0A 

-.900i368E*03 

-•PAGbROOE-Ol 

-.132740SE»I>1 

-. A992n83C400 

2912.1294 

.i93**A9ir^n4 

-,9o9l549C*03 

-.1783M6C-01 

-.137llO?E*'>> 

-.1106340E40I 

2974.0896 

,i974AXAr*nA 

-,91M5479E*03 

-.A935193C-02 

-.1402682E*01 

••1409a05E*01 

3036.0498 

.201924^r*0^ 

•.9163346E403 

•1 1939n3r-l2 

- , 1414fcl9E4l>l 

•,1796707E*01 

3098.0100 

.ao53A7ftr*o«. 

-.9145479E403 

.A9351%3e-02 

-.1402b8>E*l!' 

•,>183a08E401 

5159.9702 

•209S098t«^9 

•.9n9i349C*03 

•1 783A1SE-01 

- . 1371 102E*®1 

-.24d5n74C401 

3221 .9304 

.2l369iTr*n4 

-.9nni36AC403 

• PKGAAOOfOl 

-.132740SE»01 

•.2694205E*01 

3263.8906 

BESrAVMLE  C0?Y 

1C9 


.2l76^5?r*0^ 

• 22l7s«»nr*n^ 
,2257739c  + 0«i 
•2296l28t«09 
,2338905C*05 
,2378flb8r+09 
,2H192l8r*09 
,24S9953r*09 
,2499A74r*09 
.2540l78r+04 

• 256046ar*fi4 
.2620742r«n9 

• 286in09f:«09 
•270l24ftf ^OS 
,274l46iE*04 
,27ei70?r+os 

,282191 ^r*OS 
.26621  l4C*ri9 
•2902309r*09 
,2942497E*0«i 
,2982fc8?rl09 

• 3022A65f:«nf^ 

• S063n48(r«n9 
,3i03233r^n9 
,314342ic*.09 
,316S61<;f  4(19 
,3223Mi7r*09 
•3264n28r+09 
•3304P49r«n9 
,334448Pr+ns 
.3384728r*n9 
.34249eAr*09 

• 3465969r4>OA 
,3505459r404 
,3545A57F*09 
,3586l77e*OS 
,3626S12p4nA 
.3666A6PF*n9 
•3707P2AF*nA 
,3747a09f*0a 
,378799U4.nA 
«382839nt^nA 
,386879Art0A 
,39099l4F*OA 
.3949A35F*n9 
.3990nbnr4ns 
, 4030467^409 


• . 8A756AAF40S 
-.A715401E+03 
-.8520936C403 
-.A992097t*03 
-,8n!0862E*03 
-.7737973E403 
..74m543C*03 
-.7n6l84?E+03 
..6AAl29bE+03 
-.A9742A3E+03 
-.5A42529E+03 
-.93A768PE+03 
-.491154AE403 
•.4416000E«03 
-,3902993E*03 
-.3374534E+05 
-.9A3277nE403 
- . 29797A9F403 
-.171 777SE+03 
-.1148973E*03 
-.97? 6292E+02 
. lAA91*3E-08 
. 97«62«2E*02 
.tl489’r3E*03 
.1717775E+03 
.9P797A2E+03 
.?A*277nE403 
.33745*4^403 
.39n2993E*03 
.44l600nE403 
.491 154AE*0.3 
.93A76*PE403 
.5A42529E405 
.A9742fl3E403 
.AA*12«6E*03 
.7061842E*03 
. 741494 3e*03 
.7737973E+03 
.8fl30869E*03 
.A9®2057E*03 
.A9?0536E*03 
,A715401E*03 
, 8A758AAE+03 
.9n0l36AE403 
.9n9l 349E+03 
.9145479E403 
.9i63546E>03 


• 3S38993F-01 
,4.3979ftlE-0l 
.5238A22C.01 
.bn59463E-01 
•b8S6249t-0l 
,74»26037E-01 
,A.3657A8F-01 
,4n729A0F-0l 
.9743620e-01 
•in37A25E400 
,1 n96797t400 
, 1 1 51A4SE400 

• 1 P01947f*00 
,IP47507E400 
• 1P88144F400 

• 1 .323695F400 
,1.354n?0F*00 
, 1 .378997^  + 00 
•1398928C400 

• 1 412934F  + 00 
,142096lE+00 
,1423773^400 

• 1 420961F400 
,141P934f+00 
,1.398528^  + 00 
, 1 37A997E400 
, 1 354n?0E+00 
,l323A95r400 
,1288144^+00 
, 1 P47907p+00 
.1P01947F+00 
, 1 151A43F400 
, 1 n98797F+00 
.in37A?5F400 
.9743AP0F-01 
,9n72580E-01 

,A365788E-01 
,7K26n37E-0l 
,AM56949E-01 
,bn59483E-0l 
,S23AA2?E-01 
,4.397961E-01 
, 3938993F-01 
•P866900F-01 
,1 783918E-01 
.AS351S3E-02 
,2174229E-12 


-.1276930E+01 

-.122?606E«0l 

-,1165965E+01 

•,1107845E*01 

-,1048721E401 

-,98889iOE*00 

-.9265391E+00 

-,8677883E*00 

•.8067235E+00 

-.7454054E400 

•.6838790E+00 

-.6221785E+00 

-.5603311E+00 

-,498359lE400 

-,43628l2E^00 

-,374il39E*00 

-.3118720E+00 

-.2495692E*00 

-,1872184E*00 
-. 1248319E400 
-,6242184E-01 
-.7765209E-1P 

,6242184E-0l 
.1248319E+00 
. 1872184E^00 
.2495692E400 
. 3il8720E*00 

, 3741 139E+00 

.4362812E400 

.4983591E400 

.5603311E+00 

,622i78SE*00 

.6838790E+00 

.7454054E+00 

.8067235C+00 

,d677883E^00 

.9285391E400 

,9888910e+00 
.1048721E*01 
. 1107843E*01 
. 1165965E+01 
.122260bE«0l 
.127693nE40l 
. 1327403F401 
.157110PE+01 
. 140268PE+01 
, 1414619E+01 


-,2638367E4nl 
-,?94l392E4ni 
-.30l8330E^ni 
-.3078i84E4ni 
-.3l26423E*nl 
,31l670lE401 
,S08257bE+ni 
.3052A72E+ni 
,3026934E4ni 
. 3002A04E4ni 
.2981l26E*ni 
,P96lO80E*ni 
,294234 5E+«1 
.2924662E401 
,2907834E4ni 
,289l894E4ni 
,?876i03E401 
.2660944E4ni 
,28461 14E4nl 
,283lS23E4ni 
,26l7o87F  + r»l 
,2802727E4(M 
,2788367E401 

,2773931E+i1 
,2759340E4ni 
,2744510E4ni 
,?729^5lE4ni 
,27l3760E4ni 
,2697b20E4ni 
,2680792E4nl 
.26631 llE4nl 
,2644^74E4ni 
,2624328E*ni 
.2I»02a50E4O1 
.2578920E4ni 
,2552582E+A1 
,?522fl78E4ni 
,2468753E4ni 
,244889lE*nl 
.2400453E*01 
,2340s99E4ni 
,2263A60E4ni 
.2160f,35E*ni 
,20l6474E4ni 
,l807342E4ni 
.l505A77E4ni 
.lll8975E*ni 


3345.8508 

3407.8110 

3469.7712 

3531,7314 

3593.6916 

3655,6518 

3717,8120 

3779.5722 

3841.5324 
3903,4926 
3965.4528 
4027.4130 
4089.3732 
4151,3334 
4213.2936 
4275.2538 
4337,2140 
4399,1 742 
4461.1344 
4523.0946 
4585,0548 
4647.0150 
4708,9752 

4770,9354 

4832.6956 

4894.8556 

4956.8160 

5018.7782 

5080,7384 

5142.8986 

5204.6588 

5266.6170 

5328.5772 

5390,5374 

5452.4976 

5514,4578 

5576.4180 

5658.3782 

5700.3384 

5762.2986 

5824,2588 

5866,2190 

5948.1792 

6010,1394 

6072,0996 

6134.0598 

6196. nero 


coefficient^  fop  the  FOURIER  FIT  XG.YG  ANO  F 


X-COEFFS, 

y-COEFFS. 

F-COEFFS. 

6504981E*01 

,9163549E*03 

, 11O0839F-12 

408241 9F-08 

- ,2364052E-08 

-.1423582F*00 

6l8092nE-09 

.9164511E+03 

,2427713F-12 

1907859F-10 

,5234050F-08 

,3657487f-14 

973847qF*0l 

-.7500833E-10 

,2580546F*12 

315?1 80E-09 

,2233588E-08 

,l9i2630r-04 

1827289F-09 

.,9647972E-0l 

,29i6712f-12 

171354 3F-10 

,2n90698F-0e 

-,9321191f-15 

232nu99F-02 

-,3590699E-in 

,3577839F-12 

1 n 2u08r-09 

.2n3298lE-08 

.3275758F*07 

BESrAVAIlABLE  CC?Y 
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SPACE  OBLIQUE  HERCATOR 
FORWARD  TRANSFORMATION 
CASE  2.1 

FORWARD  transformation  tO  FIND  XiT  GIVEN  PHI.LAMOA 


PHI  = I.ISSSGO 
LAMOA  = P.375706 


x.T.z  earth-fixed  ellipsoidial  coordinates 

Xs-lBSS.tflSRXZ 
Y=  1787.339436 
Z=  5814.176682 


SCAN  VECTOR  IN  TC  PLANF  AT  TIME  T-STAR 

.405an20E*oo 

.9006060c*o0 

-.1556712E*00 


normal  vector  n at  time  t-star 

-.7870451E400 

.P872943F400 

.9138201E40n 


DIRPLACEMI  .«T  VECTOR  DELTA-R  AT  TIME  T-STAR 

-.P252466E402 

-.5021357F4I12 

.R445968r»01 


T-STAR=  5809.886615 


MAP  PhuJCCTION  cooroinates 


X 


. I/78774E*05  T 


.79i5939E»03 


BEsriv^iMBiE  copy 


SPACE  OBI.IQUE  MERCATOR 
FORWARD  TRANSFORMATION 
CASE  2.2 

forward  transformation  tO  find  x.t  given  PHI.LAMOA 


PHI  I -.0(l2B8a 

LAMOA  = -.1R2STG 


X.T.Z  EARTh-FIXED  ellipsoidial  coordinates 

X=  6SS0.26SOA8 
Y=  -719.851489 
Z=  -18.299651 


SCAN  VECTOR  IN  TC  PLANE  AT  TIME  T-STAR 

.lll6?66r»00 

.9813602F400 

-.I564344E*00 


normal  vector  n at  time  T-STAR 

.9935n73r*00 

-.n36884F*00 

-.42663CI5E-02 


displacement  vector  delta-r  at  time  T-STAR 

-.6470795F-01 

-.5472793r*02 

.87294llE*nl 


T-STARs  1553.239654 


MAP  phujection  coordinates 


X R .tO09966E*05  T = -.595e681E«02 


BESirAVAIWBLE  COPY 


SPACE  OBLIQUE  MERCATOR 
INVERSE  TRANSFORMATION 
CASE  3.1 


inverse  transformation  to  find  PHit  lamda  given  X«Y 


X=  .ln099P.6E*n5  Y=  -.59586e0E*02 


VECTOR  OElTA-R  map  PROJECTION  PLANE  COMPONENTS  AT  TIME  T-STaR 

X-COMPONENT  = -.4017802E+01 

Y-COmPONENT  = -.5565014E+02 


normal  VECTOR  N AT  TIME  T-STAR 

.9935n73E+00 
-.1 136d84E+00 
-.42b6305E-02 


vector  T-PRIME  map  PROJECTION  PLANE  COMPONENTS  AT  TIME  T-STAR 


X-COmPONENT  = 
Y-COmPONENT  s 


-.9974039E+00 

.7P01009E-01 


FINAL  PArtIALS  USFD  TN  THF  PHI,  LAMOA  INVERSION 


OX/0(PHI I 
DX/OlLAMDA) 
OY/0{PHI 1 
OY/OjlAMDA) 


-.63l2ei2E^04 
-.5418130E+03 
-.5379900E+03 
.fa3'54618E  + 0R 


THE  (PHI, LAMOA)  OF  THE 


wAP  plane  point  «X»Y) 


PHI 

LAMDA 


-.2A8fl447E-02 

-.1P25762E+00 


T-STARz  1553.259635 


bbovaiuble  copy 

I 
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SPACE  OBLIQUE  MERCATOR 
INVERSE  TRANSFOP*MATlOM 
CASE  3.2 


INVERSE  tRANSFOhmATTON  tO  PjND  PhI»  LAMqa  GIVEN  X.Y 
X=  .37787l7E+n5  Y=  .7913626E+03 

VECTOR  delta. R MAP  PROJECTION  PLANE  COMPONENTS  AT  TIME  T-STAR 

X-COmPONENT  = .1529013E+01 

Y-COmPONENT  s -.5565626E+02 

normal  vector  n at  time  T-STAR 

-.2871243F+00 

.2fl732l5r+00 

.9l37ft66E*00 


VECTOR  T-PRIME  MAP  PROJECTION  PLANE  COMPONENTS  AT  TIME  T-STAR 

X-COMPONENT  = -.9996228E+00 

Y-COmPONENT  = -.2746207E-01 


FINAL  PArTIALs  USED  IN  THE  PHI,  LAMOA  INVERSION 

DX/OJPHI)  = .6007216E+04 
nX/0(LAM0A)  = -.87fl5523E+03 
DY/DiPhU  = -,2i76227E+04 
DY/OilAMDA)  s -.2424541E+04 


THE  (PHI.LAMCAJ  of  the  "^AP  PLANE  POINT  (X.Y) 

PHI  = , ll5*.T77t*0l 

LAMOA  = .237f>793E  + 0l 


T-STARs  5809.798322 


BEST  AVAILABLE  COPY 


SPACE  OBLIQUE  MERCATOR 


DISTORTION  ERROR  ANALYSIS’* 
CASE  A.O 


LENGTH  DISTORTIONS  FOR  PTS . SYMMETRICALLY  PLACED  ON 

BOTH  SIDES  OF  THE  SATELLITE  GROUND  TRACK  FOR  THE 

DISPLACEMENT  INCREMENT  DELTA  =■  55.66  KM. 

PHI.LAMDA  OF  THE  GROUND  TRACK 

81.05172 

90.00000 

TIME  ALONG  THE  SATELLITE 

GROUND  TRACK 

= 0.00000 

— 1 

5s  'a 

3s'  1 

3s 

.998428 

1.000498 

.998941 

1.000219 

.999299 

1.000177 

.999605 

1.000114 

.999834 

1.000052 

(SATELLITE  GROUND  TRACK) 

.999938 

1.000031 

.999697 

1.000151 

.999320 

1.000340 

.998791 

1.000603 

PHI.IJWDA  OF  THE  GROUND  TRACK 

44.46685 

5.65505 

TIME  ALONG  THE  SATELLITE 

GROUND  TRACK 

- 774.50250 

1.000538 

.998860 

1.000299 

.999362 

1.000131 

.999718 

1.000032 

.999930 

1.000000 

l.nOOOOO  (SATELLITE  GROUND  TR/\CK) 

1.000031 

.999931 

1.000124 

.999725 

1.000275 

.999385 

1.000480 

.998915 

* 

The  first  (and  last)  point  at  *•*81. 05®  indicates  a slight  numerical 
difficulty  In  that  small  distortions  are  present  on  the  ground  track*, 
whereas  zero  distortion  Is  present  on  the  gro\»nd  track  elsewhere. 
Including  We  suspect  that  errors  In  the  Fourier  series 

approximation  (being  worse  at  either  end  of  the  fit)  are  the  culprit 
here.  It  h.as  also  been  observed  that  those  partlals  Indicate  small 
but  significant  loss  of  conformality  (which  is  counter  intuitive). 
This  Is  a point  which  should  receive  careful  attention  in  future  work 
on  this  problem. 
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BESI’AVAIUBIE  COPY 


PHt.L'I'OA  OF  THF  GHauNn  IrtUCK  -fa.t 

TI"E  ALOf'G  THE  GuTEILITE  GROUNn  TRACK  = 1550.70146 


1.000000  (SATELLITE  GROLlNo  THACKI 


PHT.LAvdA  of  the  GROUNO  track  -44.6987?  -18.6 

TT“E  ALONG  THE  SATEILITE  GROUNO  TRACK  = 2327.5970? 


l.OonOOO  (SATELLITE  GROUN()  TRACK) 


P-ir  ,LA»r)A  OF  THF  GRiiiiM  T“VK 


-Bl. 04555  -105.11219 


tt'E  ALfHHi  THE  4«TE(LirE  GBnU((n  TRACK  = 3103.78611 


l.oonoon  (satellite  grohNo  track) 


BEST'AVAIUBIE  COPY 
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PHT.L^'-OA  op  THC  ground  track  -4‘». 14007  I72.5e5b8 


TI'^E  along  the  SATEILITE 

GROUND  TRACK 

= 3878.27456 

1 , 0004H2 

.998914 

1,000275 

.999384 

1,000124 

.999724 

1,00003? 

.999931 

1,000000 

1,000000 

(SATELLITE  ground 

1,00003? 

,999950 

1.000131 

.999718 

1 ,000299 

.999361 

1,000539 

.998860 

Pht.lA^DA  of  the  ground  track  .23428  IfeO. 53058 

TI^'E  ALONG  THE  suTElLITE  GROUND  TRACK  s 4651.073fcl 


1 , 0005srt 
1 , 00031  *4 
I.OOOIhu 
1 . 00003n 
I • QOOOOU 
1 ,00003*1 
1 ,0001 ^7 
1 , 00C3l  •. 
I ,0005-*rt 


.998639 

.999346 

.999709 

.999927 

1.000000  (SATELLITE  GROUND  TRACKJ 
. 999927 
, 999709 
.999346 
. 998639 


BESrAVAIUBLE 


117 


PHT.LAHOA  of  THC  ground  track  44.56X29  14d, 42096 

TT'»E  along  the  f;ttTElLITE  GROUNn  TRACK  = 5423.16305 

l,n0053M  .998860 

l,n0029S  .999362 

1.000131  .999718 

1.0000A2  .999930 

l.OOOOOU  1.000000  (SATCLLlTe  QROUNq  TRACK) 

1.000031  .999931 

1.00012^  .999725 

1.00027**  .999385 

1.0004HU  .998915 

PHT.LA-^nA  OF  THE  GRuuND 
TI^'E  ALONG  THE  9«TEI 

.99fl4>'7 
.9«»89^l 
.99929-< 

.9996.1h 
.9998  ^*4 
, 9999k9 
.90999r, 

, 9S99n  3 
,999fa«*. 


BESr'AVyilMBLE  COPY 


TRACK  el*05l72  64.12310 

LITE  ground  TRACK  = 6195.99195 

1 .000498 
1 .000219 

1 .0001  77 

i.ocoim 

1.000052  (SATELLITE  GROUND  TRfirK) 
I .000010 

1.000001 
1 .000056 
1 .000121 


I 
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SPACE  MERCATOR 

POHWARIJ  TKANJsEORMATION 

CASE  S.l 

F0P<^A«0  tWANSFOwMATION  tO  FIND  X.Y  GIVEN  PHI.tAMDA 


HHI  s l.l'SSAGO 

LAMnA  s ?,375706 


XfV.2  TARTm-FIiFD  rLUlPSOIOTAL  COORDINATES 

Xs-lA5fl.4A3R^2 
Ys  l7a7.53943G 
2=  bfl1H.17fa6A2 


SCAN  vector  TN  TC  PLANF  at  time  t-star 

,i*o%Rn20r*no 

.9006n60F*n0 

*,l55fe7l2r*00 


normal  VrCTOR  f'*  AT  TIME  T-STAR 

-.?07O45ir»OO 

,2e7?9**5r»oO 

.9i36?0ir+00 


OISPlACEM^-,<T  VECTOR  DFUTA-H  AT  TImE  T-STAR 

-.225?4fc6r»02 

-.A021^^7F*02 

,fl44S9bAr*nl 


T-STAKs  %eOP. 075747 


i"AP  Projection  COorotnatcs 


I s .A77A7A7E*05  Y = 


,7913S01E*03 


BEST  AVAIWBLE  COPY 
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SPACE  OBLIQUE  MERCATOR 
FORWARD  TRANSFORMATION 


CASE  5.2 

forward  TRAW<;FOMMflTTON  tO  FIND  X.r  GIVEN  PHI.LAMDA 


PHI  : -.OnZASS 

LAMnA  s -.IPaSTG 


X.Y,2  FARTH-FIifED  FLLIPSOIOTAL  COORDINATES 

X=  6330.2B3C.^8 
Y=  -779.fl‘Sl4S9 
2=  -18.2996A1 


SCAN  VErrOR  IN  »C  PlANF  at  TIME  T-STAR 

.1  H6?66F  + nO 
.9813602F+00 
- . I564344r*00 


normal  vector  u at  TIWE  T-STAR 

.993Sn73r+00 
-.1 13Afl84F*00 
• .u26A308r-0? 


DISPUACf  Mf  i4T  VECTOR  OFLTA-R  AT  TIME  T-STAR 

-.8470795F*nl 

-.«>472793F*ft2 

.A729ullr*r»l 


T-STARs  I553.22a*;33 


MAP  PmiJJEPTION  coorotnates 


X = . 1 n0RR»i9E»0b 


.•>9«i7649E*02 


BESrAVAIUBlE  COPY 


SPACE  OBLIQUE  MERCATOR 
INVERSE  transformation 


CASE  6.1 

transfohmatton  to  find  Phi.  lamda  given  x.t 
. in.n99«i9E  + n5  Ys  -,5957649E  + 02 

vrrTCR  OElTA-R  map  PkOjFCTION  plane  components  at  TImE  T«STaR 

X-COMPONrNT  s -.4017801E+01 

Y-COmPONPNT  = -.5*>650llE  + 02 


normal  vector  N at  tImE  T-STAR 

.9935080E*n0 
-.1 156R55F+00 
-.‘♦243970r-02 


vrcToP  T-pRIME  Map  PRrMECTION  PLANE  COMPONENTS  AT  TIME  T-STAR 


X-COmPONFMT 

y-componfnt 


-.9974039E+00 

.7?01010E-01 


FINAL  partial^  USm  TN  THF  PHI.  LAMDA  INVERSION 


nx/0(PHi » 
nx/OfLAMnA ) 
OY/OlPHl ) 
DY/OILAMOA) 


- .b3l26l2E«04 
-.5U18128E+03 
-.5379899E+03 
,b354bl6E4-04 


THF  (PHl.LAMOAl  THE  MAP  PLANE  POINT  (X.Y» 


PHI  = 
I AMDA  S 


-.?8aA484E-0? 

-,l?2S7b2E+0n 


T-STARr  I5*t3,2l74bb 


BESr  ^VAIMSlECOry 
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SPACE  OBLIQUE  MERCATOR 
INVERSE  TRANSFORMATION 

CASE  6.2 

INVERSE.  TRaNSFOhMATTON  TO  FIND  PHl«  LAMQA  GIVEN  X»Y 


X=  ,37707A7E+n5  Y=  ,7913901E^03 

vrrTQR  nELTA-R  map  mkojfctton  planf  components  at  time  t-star 


X-COMPONFNT  = . l!S2fl720E*01 

Y-COmHONFNT  = -.5S656mE«»-02 

normal  VECTOR  N AT  TIME  T-STAR 


-.>a7n649r+no 
.2073fliir*no 
.91301 l7F*rO 


vrr r.jR  t-^^rime  map  projection  plane  components  at  time  t-star 


X.COmpONFNT 

Y-COmhONENT 


-.999b230E*00 

-.279S607E-O1 


PINAL  PAOTIALS  USEO  TN  THE  PHI, 


lamoa  inversion 


nx/o<HMi » 
nX/Oll AMOA ) 
nY/0|HMl » 
OY/DlLAMnA  » 


.6n07062E+04 

-.0705b53E*O3 

-.2l7660bE*O'4 

-.2424019E*04 


THE  CPHr,LA^nA>  UF  THE  MAP  PLANE  POINT  <XiY» 


mhI  s 
I 4MnA  s 


.115S660E*0l 

,?M7«;706E*01 


r-STARs  5SO9,0A4S9O 


best:avaiubie  copy 
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